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In considerable amounts of NO containing small 
amounts of NO. a continuous absorption occurs in the 
ultraviolet, obscuring both the absorption of NO and that 
portion of the absorption of NO» which lies below 2500A. 
The behavior of this absorption with respect to tempera- 
ture and the partial pressures of the constituents is rather 
convincing evidence that it is due to N2O;. When very 
small amounts of water are also present a group of bands 
occurs in the near ultraviolet, lying in the same region but 


not resembling the longer wave-length NO, absorption. 
These bands appear diffuse under low dispersion but 
possess an ordered arrangement. With increase of temper- 
ature the intensity of the bands decreases rapidly. They 
begin in the vicinity of 3850A, extending to shorter wave- 
lengths. The first members are broader and more diffuse 
than those that follow, indicating a predissociation process 
in the carrier, which is probably HONO. 





HEN the absorption spectrum of pure nitric 

oxide is photographed in moderate path 
lengths at atmospheric pressure and room tem- 
perature with quartz and ordinary unsensitized 
plates, the entire spectrum consists of the 0—0 
band of the y system with perhaps the 0—1 band 
also showing at the extreme short wave-length 
limit. If a very small amount of oxygen or of 
nitrogen dioxide be present, however, a con- 
tinuous absorption overlies this region of the NO 
bands, which is due to the molecule NOs. If, in 
addition, a very small amount of water vapor be 
present, a series of bands appears in the region 
4000-3000A which are not characteristic of the 
NO» spectrum. nor, of course, of that of NO. 
The authors have earlier reported! discovering 
these absorptions during work on the optical 
dissociation of N.O, and at that time erroneously 
attributed the bands also to the carrier N2O3. The 
presence of water vapor as well as both the other 
constituents is necessary for the production of 


' Melvin and Wulf, Phys. Rev. 38, 2294 (1931). 


the bands. A probable carrier of the bands would 
seem to be the molecule HONO, though a 
polymer of this molecule would also be a pos- 
sibility. In this paper they will be referred to 
as due to HONO, or, as is more often written, 
HNOsz. 

The absorption bands are shown as the two 
lower portions of Fig. 1, with the absorption of 
NO, above for comparison. Though the positions 
of such diffuse bands can hardly be given wave- 
length values with useful accuracy, these bands 
nevertheless suggest that they comprise two 
progressions, as is diagrammatically illustrated 
below the absorption in the figure. The bands 
set in remarkably close to the position of the 
first predissociation process in NOz. The first 
bands are very diffuse, sharpening noticeably 
as they proceed to shorter wave-lengths. 

If these bands are due to HONO, or a polymer 
of this, they should show a dependence upon the 
pressures of the three constituents NO, NO» and 
H,0, and probably upon the temperature of the 
mixture, depending upon the magnitude of the 
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heat of dissociation of HONO. The continuous 
absorption due to N2Os; should similarly show a 
dependence upon the pressures of NO and NOs, 
and also upon the temperature of the mixture, 
this latter following from the fact that the heat 
of dissociation of N2O; into NO and NO¢z has been 
found to be about 10,300 calories.? All of these 
conditions are fulfilled. 

In Fig. 2 the upper picture shows these ab- 
sorptions at three different partial pressures 
of NOs» in a 40-cm path length and constant 
partial pressure of H,O and NO, the first strip 
being the light source through the empty cell, 
the next lower being taken with H;,O at a pres- 
sure of 16 mm Hg and NO at roughly 544 mm 
Hg (this amount of NO is so great that the 
v’’=1-—>v'=0 band shows in absorption at the 
short wave-length end, even at room tempera- 
ture). The last three strips show the absorption 
for approximately 0.5 mm, 1.5 mm and 4.5 mm 
Hg pressure, respectively, of added NO». The 
lower picture in Fig. 2 is similar, showing the 
light source through the empty cell, through the 
cell containing 7.4 mm of NOs, and thirdly with 
5.2 mm of HO added, while in the three lower 
strips NO was added to the cell in three increas- 
ing amounts, roughly 35 mm, 215 mm, and 760 
mm Hg, respectively. 

In Fig. 3 the upper picture shows in a similar 
way the influence of water vapor, the first strip 
being the source through the empty cell, the 
second with the cell containing approximately 
7.4 mm of NO» and the third with 470 mm of 
added NO. The three following strips were taken 
after the addition of water vapor in amounts 
roughly 0.9 mm, 6.1 mm, and an amount which 
would have given 41.5 mm Hg had not some 
undoubtedly condensed out on the walls. In the 
lower picture of Fig. 3 these absorptions are 
shown at three different temperatures, the NO: 
pressure being 1.5 mm Hg, the NO pressure 380 
mm Hg, and the H,O pressure 20 mm Hg. 

It should be stated that the introduction of 
the gases into the absorption cell always in- 
volved condensing out in a small U-tube, con- 
tained in the circulating system which was in 
series with the absorption cell. This means that 


— and Daniels, J. Am. Chem. Soc. 53, 1250 
(1931). 
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the solid and liquid phases were always present, 
for a short time at least, prior to the above ex- 
posures. An all-glass magnetic piston-type pump 
effected the circulation and resultant mixing. 
The various amounts of gas added to the cell 
are given above as pressures they would have 
exerted had they been alone in the absorption 
tube system. 

There is a further observation which lends 
weight to the point of view that the bands are 
due to HNOsz. If, to a cell containing a small 
amount of NOs, water vapor be added in excess 
of saturation until a small amount of liquid 
water must exist as a film on the walls, the NO» 
may be reduced to a value so low that it cannot 
be seen. This is probably due to the formation of 
HNO; and NO, the NO making its appearance 
in absorption. If now NO be added in excess 
the HNO, bands become plain, the NO effecting 
an increase in the concentration of the HNO: 
which may be many fold, depending upon the 
degree of dissociation of the HNOsz prior to the 
addition. 

Certain of the bands give the impression of 
being related to the NO» absorption, because of 
the close correspondence of their positions with 
some of the WO, bands. However, in the case of 
such diffuse bands and such relatively low dis- 
persion this may well be fortuitous.* 


DISCUSSION 


At least two suggestions might be made as to 
the significance of the two groups of bands. 
They might be thought of as transitions in but 
one degree of oscillational freedom for which 
v’=0 and v’=1. This does not seem probable, 
however, because the relative intensities of the 
two progressions to be expected on such a basis 
do not seem to exist. With increasing tempera- 
ture the intensity of both progressions decreases 
rapidly as would be expected owing to the 
decrease in the concentration of the HONO 
through thermal dissociation. There is, however, 
no noticeable tendency for the weaker progres- 


3 If this were a real correspondence one might ascribe 
the bands to a selective enhancement of certain of the NO» 
bands by the mutual action of NO and H.O. It seems, 
however, that such an explanation would not differ greatly 
from attributing the bands to HNO; or a polymer thereof. 
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Fic. 1, Absorption bands occurring in mixtures of NO, NO. and H,O, probably due to HNO». The absorption of 
NO, is shown above. 


SOURCE 
NO+1L:0 
+NO» 
+NO:z 
+NO 


Fic, 2A. The effect of increasing nitrogen dioxide, showing at left the continuous absorption of 
N,Os, and at the right the absorption bands of HNO,. 


SOURCE 

| NO2 

» +110 
+NO 

, 4+NO 
+NO 


Fic. 2B. Similarly, the effect of nitric oxide. 


SOURCE 
NO: 
+NO 
+H:0O 
+H:2O 


oak d “ cae +-H:0 





Fic. 3A. The effect of increasing water, showing the decrease of the NO; absorption, and the 
increase of that of HNOsz. 


Fic. 3B. Effect of increasing temperature, showing the disappearance of NoO; and HNOs, and 
the appearance of NO and NOs. 
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sion of the bands to decrease less rapidly than 
the other progression, which on the interpreta- 
tion mentioned at the head of this paragraph 
would be brought about due to the partial off- 
setting of the decrease in the concentration of 
the carrier by the increase in the Boltzmann 
quota of molecules in the vibrational level 
v’=1. If the bands did represent two progres- 
sions for which v’’ = 0 and v”’ =1, then a minimum 
value for the initial separation of the levels in 
this fundamental frequency would be roughly 
750 cm~!. Now one may compare the change in 
the ratio of the numbers of molecules in the 
zeroth and first vibrational states at 298°K 
and 473°K with the relative intensity of the 
bands at these two temperatures. Assuming equal 
statistical weights, this frequency separation 
leads to a ratio of populations in the two states at 
room temperature of roughly 0.027, while at the 
higher temperature it is about 0.103. This 
change in relative intensity of only 4-fold might, 
however, be present in such a mixture of ab- 
sorptions without being noticeable. However, 
the ratio of intensities to be expected at room 
temperature, namely, about 40-fold, also hardly 
seems to be substantiated in fact, though in the 
absence of quantitative photometric measure- 
ments it cannot be ruled out. It is perhaps doubt- 
ful if even the latter would be very conclusive in 
this mixture of absorptions, and with the NO» 
absorption itself changing somewhat with the 
temperature. 


TABLE I. Absorption bands of HNO». 





WAVE NUMBER DIFFERENCES 
(cm~!) (cm~!) 


3843 26014 


WAVE-LENGTH 
(A) 





3681 27159 
3539 28248 
234 
3510 28482 
784 
3416 29266 
241 
3388 29507 
723 
3307 30230 
268 
3278 30498 
704 
3204 31202 
265 
3177 31467 
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Fic. 4. Diagrammatic potential energy curves. 


An alternative hypothesis seems in_ better 
accord with the observations, namely, that the 
transitions all originate from the same normal 
level and go to several levels in one mode of 
oscillation in the excited state and to two levels 
in a second mode of oscillation in the excited 
state, thus yielding the two partial progressions. 
From Fig. 1 it can be seen that, on the basis of 
the first hypothesis, there is the opportunity to 
obtain three estimates of the separation v’’=1 
to v’’=0, while on the basis of the second 
hypothesis there is the opportunity to obtain four 
estimates of the separation ve/=1 to v2’=0. 
Table I gives the data for these bands and these 
two sets of differences. In this comparison lies 
also a slight indication in favor of the second 
hypothesis, since in the case of the first estimate, 
instead of being constant the three observed 
values decrease progressively, being 784 cm", 
723 cm, and 704 cm™!, respectively, while in 
the second case the four observed differences are 
234 cm=', 241 cm™', 268 cm~! and 265 cm™!. The 
failure of the former to show constancy argues 
against the first hypothesis, and everything 
considered the second seems the more probable. 
The authors obtained similar results in a study 
of the effect of temperature upon the absorption 
bands of ozone.‘ In this case the bands are very 
numerous, and several levels of at least two 
modes of vibration in the excited state are evi- 
dent. In spite of this no bands were observed 
which came from levels of the normal state 
appreciably higher than the lowest. Since this 
work, Chalonge and Lefebvre,5 working with 


4 Wulf’and Melvin,’ Phys. Rev. 38, 330 (1931). 
*Chalonge and Lefebvre, Comptes rendus 197, 444 
(1933); ibid. 199, 456 (1934). 
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very large effective paths of ozone, have found 
bands at longer wave-lengths which come from 
higher vibrational levels in the normal state. 
The beginning of the HONO bands so nearly 
at the point of the first predissociation of NO» 
suggests that there may be a connection between 
these two phenomena. The exceptional diffuse- 
ness of the first bands indicates that at about this 
energy there exist discrete vibrational states of 
an electronic level in HONO from which there 
is a high probability of transition into another 
electronic level in the molecule, in which, for 
this energy, the molecule is dissociated, whereas 
for vibrational states higher than these the 
probability of transition into this second level 
grows less. Fig. 4 shows diagrammatically one 
arrangement of levels in the molecule that could 
lead to such circumstances,® pictured as if but 
one of the fundamental frequencies of this 
polyatomic molecule were excited, that is, as 
if one were dealing with a diatomic molecule. 
Franck, Sponer and Teller’ have discussed pre- 
dissociation in polyatomic molecules and have 
explained important differences from the be- 
havior of diatomic molecules that may in general 
arise. Gradual setting-in of predissociation is 
one of the important differences. In diatomic 
molecules it is very improbable, to be expected 
only where the potential curve of the upper level 
and the perturbing level intersect in an extremely 
acute angle. In polyatomic molecules, however, 
where the relative positions of the atoms are 
determined by three or more parameters, factors 
arise in the complicated motion of oscillation 
which can lead to a gradual setting-in of pre- 
dissociation, and here it is frequently en- 
countered. But the present case, in contrast, 
shows a gradual decrease- of diffuseness as one 
passes to higher v’ levels. This behavior is in 
general less to be expected in polyatomic mole- 
cules. It becomes possible, however, if the ex- 
citation is chiefly in one degree of freedom, and 
this, indeed, seems to be true in the present case. 
The experimental observation that the bands, 
starting very diffuse, become sharper with in- 


_*Franck and Sponer, Nachr. d. Gesell. d. Wiss. zu 
Gottingen, p. 241 (1928). 

‘Franck, Sponer and Teller, Zeits. f. physik. Chemie 
B18, 88 (1932). 
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creasing frequency seems to be quite certain. It 
is to be noted that both the increasing dispersion 
and, for the first few bands also the increasing 
intensity, act to minimize the appearance of 
this effect; it is, however, quite apparent in 
spite of these. This appears to be one of the 
simplest and most apparent cases of this phe- 
nomenon. How far this sharpening of the bands, 
after the first few, is continued, cannot be said, 
for they rapidly become faint and are lost in 
the continuous absorption of NOs. 

There are difficulties in the way of obtaining 
the bands more intensely. Temperature change 
is of little help, an increase decomposing the 
HONO, and a decrease reducing its concentra- 
tion through the reversion of NO2 to N2O, and 
the decrease in the saturation pressure of water. 
Increasing the absolute amount of NO: ob- 
scures the region by its own absorption. The 
factor which promises a decidedly favorable 
change is the pressure of NO. If this could be 
raised considerably above one atmosphere, an 
appreciable gain in the concentration of HONO 
molecules would be expected. 

It would be pertinent to raise the question as 
to the N2O; content of ‘‘pure NOs,” for in a 
mixture of as pure NO,»—N.O, as one can 
prepare there exists the equilibrium 2NOz2 
=N.03;+30:2. Using the free energy value given 
by Verhoek and Daniels? for the reaction 
NO.+NO=N,03, AFo9se= +441 calories, and 
that given in the International Critical Tables* 
for the reaction NOs= NO+40Oz, AFoo3° = +8374 
calories, one immediately sees that at room 
temperature this amount of N.2O; will be very 
small. For NOz pressures of a few millimeters the 
partial pressure of N.O; is of the order of 10~° 
that of NO». 

It is interesting that, in the path lengths we 
have used, the absorption in the red or orange 
that leads to the blue color of liquid N2O; has 
not yet been observed. It may be that longer 
paths will show it. It may also be that this 
absorption is due to a higher association product 
such as N,Os¢, appearing with appreciable 
strengths only at high densities as in the liquid. 


8 International Critical Tables, Vol. 7, p. 239. 





































DECEMBER, 1935 JOURNAL 


OF CHEMICAL 


PHYSICS VOLUME 3 


The Raman Spectra of the Isotopic Molecules H., HD, and D, 


Gorpon K. TEAL AND GEORGE E. MacWoon, Department of Chemistry, Columbia University 
(Received August 29, 1935) 


The Raman spectra of Hz, HD and Dz gases, at 3 
atmospheres pressure, were excited by Hg 2537 and photo- 
graphed with a high dispersion quartz prism spectrograph. 
A total of 42 lines were analyzed. The positions of the pure 
rotation lines for Hz agree with the measurements of 
Rasetti but the lines in the Q vibration band do not. 
New vibrational constants calculated from the Raman and 
emission spectrum data of Hz and HD give calculated 


values of the positions of the rotational-vibration Raman 
lines in good agreement with the values observed for the 
three isotopic molecules. Contrary to the theoretical 
predictions, the 0,0 line of the Q vibrational band of HD 
was observed. A source of error in Raman measurements, 
namely, the movement of the spectrum due to variations 
of the barometric pressure during the exposure of a plate, 
was detected and eliminated. 





INTRODUCTION 


HE large mass ratios and low moments of 
inertia of the three hydrogen molecules Hoe, 
HD, and D2 make them particularly suited for a 
demonstration of isotopic effects expected in the 
Raman spectra of diatomic molecules. Several 
years ago the Raman spectrum of the He 
molecule was investigated by Rasetti,' who suc- 
ceeded in photographing lines in the pure 
rotational and the 0-1 rotation-vibration bands. 
Observation of two rotational lines and the 
unresolved Q vibrational branch of Dez was 
recently reported in a note by Anderson and 
Yost.? In this paper we wish to report measure- 
ments of lines in the pure rotation and the 0-1 
rotation-vibration bands of the He, HD, and D» 
molecules. 

In the course of this work, a serious difficulty 
was encountered, namely, a shift in the spectrum 
due to fluctuations in the barometric pressure. 
Merton,* working in the visible with a high 
dispersion glass spectrograph, appears to be the 
only one to have reported this effect. Although 
we have found it a serious source of error in 
Raman measurements, especially for gases where 
the exposure times are longest, a search has 
revealed no previous mention of it in the litera- 
ture on the Raman effect. For this reason we 
report it in rather greater detail than might 
otherwise seem necessary. 


1F. Rasetti, Phys. Rev. 34, 369 (1929). 

2 T. F. Anderson and D. M. Yost, J. Chem. Phys. 3, 243 
(1935). 

3T. R. Merton, Proc. Roy. Soc. A113, 704 (1926). 


EXPERIMENTAL 


The first photographs of the Raman spectrum 
of hydrogen, taken in the region 2450 to 2850A 
with a Hilger E-1 quartz prism spectrograph, 
were entirely unsatisfactory for the measurement 
of the Raman displacements. The lines on two of 
the plates appeared to be incompletely resolved 
doublets and on others they appeared diffuse and 
the rotational structure of the vibrational branch 
was incompletely resolved. This loss of resolution 
was found to be due to fluctuations in the 
barometric pressure during the exposures. Such 
fluctuations are accompanied by changes in the 
refractive index of the air in contact with the 
prism and by consequent changes in the path of 
the light entering and emerging from it. By 
photographing the iron arc spectrum through 
adjacent holes of a Hartmann diaphragm, it was 
found that shifts of the order of 0.1A frequently 
occurred during the three hours between suc- 
cessive exposures. Measured values of the shift 
observed in this way agreed closely with those 
calculated from the change of the index of 
refraction of the air with pressure. The shifts 
were too large to be accounted for on the basis of 
a pressure effect on the frequency of the lines 
emitted by the iron arc. 

To eliminate the above difficulty, the spectro- 
graph was encased in a pressure-tight steel box 
having removable ends of aluminum. A crystal 
quartz window admitted the light to be analyzed. 
A rack and pinion device was provided for 
changing the position of the Hartmann dia- 
phragm without opening the box. The instrument 
was barostatted at 739 mm by use of a mercury 
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Fic. 1. Raman cell. 


manometer contact and a relay controlling a 
Hyvac oil pump connected to the box. The 
pressure was maintained constant to within less 
than 0.1 mm. The instrument was also thermo- 
statted at 32°C to within a few hundredths of a 
degree. 

Subsequent to the barostatting of the spectro- 
graph, additional exposures were made to de- 
termine more exactly the magnitude of the 
pressure effect. During two successive exposures 
of a plate to the mercury arc, the pressure in the 
spectrograph was maintained at 595 mm and 
then at 745 mm. Measurements of the displace- 
ment of corresponding lines in the superimposed 
spectra gave a value of approximately 2.6A/cm 
change in pressure as compared to the calculated 
value of 1.7A. The discrepancy between the 
observed and calculated values is probably to be 
attributed to the approximate character of the 
calculations and to strains set up in the spectro- 
graph mounting. 

The fused quartz cell in which the gas under 
study was placed was constructed as indicated in 
Fig. 1. It was approximately 60 cm long and had 
an inside diameter of 2.5 cm. The walls of the 
cell, which were approximately 2 mm thick, were 
calculated to withstand, with a safety factor of 
three, two atmospheres in excess of the atmos- 
pheric pressure. One end was drawn down 
symmetrically along the axis of the cell to the 1 
mm (I.D.) constriction A, to which the sub- 
sidiary light horn B was added. The Pyrex 
vacuum stopcock C was attached to the cell by 
means of a seven-step seal. The gas contained in 
the cell was illuminated by three low pressure 
Hanovia 2537 Mercury Arcs, the radiating 
elements of which were quartz tubes (1 cm in 
diameter and 133 cm long) bent into the shape of 


closed S’s. The arcs were held firmly against the 
Raman tube by a polished aluminum reflector 
which was bound around them. Each arc was 
supplied approximately 100 milliamperes by 
individual high voltage (0.5 kva) transformers. 
The arcs maintain a low temperature and there- 
fore were sufficiently cooled by the use of two 
fans. 

The illumination of the spectrograph slit by 
light coming directly from the walls of the quartz 
cell was minimized by a series of three platinum 
diaphragms D spun and spot-welded into a 
platinum sleeve and blackened by the elec- 
trodeposition of a thin film of platinum by the 
method of Lorch.4 The diaphragms allowed 
uniform illumination of a 3-mm length of the 
slit when the tube was placed about three 
centimeters from it. The cell was easily and 
accurately aligned along the optical axis of the 
spectrograph by sighting through a temporarily 
unblackened portion E of the subsidiary light 
horn and through the constriction A and the 
diaphragms D at light sent through the spectro- 
graph in the reverse direction. 

The mercury arcs were found to produce high 
concentrations of ozone resulting in the ab- 
sorption of \2537 and also the scattered light. In 
order to reduce this loss of light, a moderate 
stream of nitrogen was allowed to flow between 
the cell and the reflector. It is estimated that 
this method of eliminating the envelope of ozone 
around the cell decreased the exposure time by a 
half or two-thirds. The actual exposure times 
were from 100 to 200 hours. 

The exciting line that was scattered without 
change of wave-length was partially absorbed by 


4A. Lorch, Thesis, Columbia University, 1932. 
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vapor from a small quantity of mercury placed in 
the spectrograph. 

The deuterium content of the gas samples used 
for obtaining Plates I, II, and III was 99, 70, and 
60 percent, respectively. 

The iron arc was used as a reference spectrum. 


RESULTS 


In Figs. 2 and 3 are shown the microphotometer 
curves of the pure rotational spectra on Plates I 
and II. Plates I and III were taken after the 
spectrograph was barostatted. Due to the failure 
of the relay controlling the pump connected to 
the barostat, the pressure in the spectrograph 
was reduced from 739 mm to about 600 mm for 
approximately 10 to 12 hours during the exposure 
of Plate I. The resultant pressure shift of one of 
the mercury lines is indicated in the micro- 
photometer curve of Fig. 2 as an apparent change 
in wave-length Ad. Plate II is the only one 
satisfactory for measurement obtained previous 
to the barostatting of the spectrograph. Even for 
this plate, however, the line width was somewhat 
greater than obtains when the pressure is 
controlled, as may be seen by comparing the 
curves of Plates I and II in Figs. 2 and 3, 
respectively. 

While photographing the iron spectrum on 
Plate I an accidental jarring of the spectrograph 
caused a movement of the plate and resulted in a 
shift of the reference spectrum with respect to 
the Raman spectrum. In determining the fre- 
quencies of the Raman lines, allowance for this 
shift was made by determining its value for the 
Rayleigh scattered mercury lines. It was con- 
sidered best to obtain the apparent frequency of 
the irradiating mercury line from the Stokes and 








Fic. 3. 


anti-Stokes Raman lines on each plate, the values 
39,411.9 cm@!, 39,412.43 cm and 39,412.35 
cm being found for Plates I, II, and III, 
respectively. 

The error in measurement of the lines varies 
from 0.5 to 1.0 cm~ except for several very weak 
lines. It is estimated that the average values of 
the three plates are good to approximately 0.3 
cm}, 

In order to calculate the frequencies, the 
following expression for the energy® was used: 


E/hc=6(v+}) —x@(v+ 4)? 
+yaHe(v+3)> —20-(v+3)* 
+[B.—a-(v+3)+7e(v+ 3)? —6.(0+3)*? JJ(J +1) 
+[D.+B.(0+2) JJ?(J+1)P+ FF(J +1): 


The constants used are those of Jeppesen®: 7 as 
recalculated by Urey and Teal.* Below is given, 
in Tables I and II, a comparison of the calculated 
and observed values of the frequencies in wave 
numbers in vacuum. 

It will be seen in Table I that there is good 
agreement between the observed and calculated 
values of the pure rotational lines. Though a 
trend in the values for De indicates that the 
rotational levels for this molecule may be closer 
together than predicted by theory the variation 
from the theoretical values is within the experi- 
mental error for most lines. Very good agreement 
is had between the values for Hz in Table I and 
those found by Rasetti! and Bhagavantam.® For 


5 W. Jevons, Report on Band Spectra of Diatomic Mole- 
cules, Cambridge (1932). 

6 C. R. Jeppesen, Phys. Rev. 44, 165 (1933). 

7C. R. Jeppesen, Phys. Rev. 45, 480 (1934). 
(1935) C. Urey and G. K. Teal, Rev. Mod. Phys. 3, 243 

°S. Bhagavantam, Ind. J. Phys. 6, 319 (1931); 7, 105 
(1932). 
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TABLE I. 








DIFFER- 
ENCE 


PLATE 


PLATE 
II III AVERAGE CALC. 





586.3 
354.3 
354.5 
586.9 
815.3 


586.3 
354.3 
354.5 
586.9 
815.3 


587.0 —0.7 
354.3 
354.3 
586.9 
814.4 
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the Q vibrational branch of He, however: 
discrepancies appear between the values found in 
this work and the measurements of Rasetti. The 
possibility that the differences between these 
values and Rasetti’s might be due to the effect of 
pressure was considered since Rasetti’s observa- 
tions were made at approximately 12 atmospheres 
whereas the values in Tables I and II are for gas 
at 3 atmospheres. Bhagavantam did not observe 
an effect on the lines of hydrogen when the 
pressure was varied from 13 to 50 atmospheres, 
but no previous measurements at lower pressures, 
as used here, are available.* The cause of the 
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differences is not yet understood. The dis- 
crepancies are not cleared up by the measure- 
ments of Bhagavantam since these latter values 
agree in part with the measurements given in 
Table II and in part with those of Rasetti. 
Since agreement was had between the measured 
values obtained from three plates and the 
discrepancies between the observed and calcu- 
lated values, given in Table II, were consistent, 
it was concluded that the vibrational constants 
were incorrect. New vibrational constants were 
therefore calculated. Least-square solutions of 
the AG values for the first seven states were 
obtained for He from the data of Jeppesen,® 
Rasetti!® and this study and for HD from the 
data of Jeppesen,’ Mie™ and this study. Our 
experimental values for AG, are: 


AG;(He) =4160.2+0.2 cm~, 
AG,(HD) = 3631.4+0.2 cm“, 
AG,(Dz2) = 2993.5+0.2 cm. 


The emission spectrum data and Rasetti’s and 
the present Raman data were weighted 1, 3, and 
5, respectively. The best values of the constants 
thus obtained are listed in Table III. The 
relationship observed between the w,’s of Hz and 
HD agrees excellently with the theory and the 
agreement for the xw,’s, yw,’s, and 2w,’s, is about 
as good as could be expected from the accuracy 
of the data and from the number of levels used in 


TABLE II. 
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_"In this connection it is interesting to note that the 
nitrogen vibration line appeared on Plate III, and the 
measured value of its position was less than that recorded 
byRasetti for nitrogen at 12 atmospheres by approximately 


the same amount as observed for the hydrogen lines. 

10 Data of Rasetti as recalculated by R. T. Birge and 
C. R. Jeppesen, Nature 34, 463 (1930). 

11K. Mie, Zeits. f. Physik 91, 475 (1934). 















































TABLE III. 
pe’: He-HD 
He HD De Obs. Theory 
We 4405.30 3817.09 3117.05 0.86648 0.86616 
Xwe 125.325 94.958 63.034 0.7577 0.7502 
VWe 1.9473 1.4569 0.74526 0.748 0.650 
Zwe 0.11265 0.07665 0.031573 0.680 0.563 








calculating these constants.* The molecular 
constants of De given in the table are the average 
values calculated from those of He and HD, 
assuming the theoretical relationships between 
the reduced masses of He and D2 and of HD and 
Dz to hold. The value of AG; (De) calculated by 
use of these constants is 2993.2 cm@!, which 
differs from the experimental value of 2993.5 
cm by only 0.1 cm more than the estimated 
accuracy of the experimental value. The agree- 
ment of the observed positions of the rotational 
vibration lines of the three isotopic molecules 
and the positions calculated from the new 
constants is very good, as may be seen from 
Table II. 

The measurement of the heights of the in- 
tensity maxima of the microphotometer curves of 


* It seems improbable that the theoretical relationships 
between the constants of the isotopic molecules could hold 
exactly due to differences to be expected in the potential 
energy curves from the differences in the electronic zero- 
point energies of the hydrogen and deuterium atoms as 
well as the molecules. This effect, however, should be small 
and not detectable in the constants in Table III. 
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Plates I and II shows an alternation of intensities 
of the pure rotational lines of De, the lines 
involving transitions from even states being 
more intense. This, in agreement with previous 
work on emission spectra,’ allows one to 
conclude that the deuterium nucleus obeys the 
Einstein-Bose statistics and has a spin different 
from zero. The semi-quantitative measurements 
agree very well with those expected according to 
Placzek’s theory’ if Murphy and Johnston’s 
value of the spin equal to 1 is used. Though the 
theory predicts the absence of lines corresponding 
to V :0-1, J : 0-0 transitions a line for this 
transition for HD was observed on Plate III. 
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Any lattice in which the hydrogen atoms would be 
translationally identical (Bravais lattice) would have 
metallic properties. In the present paper the energy of a 
body-centered lattice of hydrogen is calculated as a 
function of the lattice constant. This energy is shown to 
assume its minimum value for a lattice constant which 
corresponds to a density many times higher than that of 
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the ordinary, molecular lattice of solid hydrogen. This 
minimum—though negative—is much higher than that of 
the molecular form. The body-centered modification of 
hydrogen cannot be obtained with the present pressures, 
nor can the other simple metallic lattices. The chances 
are better, perhaps, for intermediate, layer-like lattices. 
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METALLIC MODIFICATION OF 


(1) INTRODUCTION 


HE present calculation of the constants of 
metallic hydrogen has a threefold purpose. 
First, we wanted to explore the possibility of 
preparing metallic hydrogen under high pres- 
sures.' The hope for this, however, appears to be 
small, as the pressures under which metallic 
hydrogen would be stable seem to be unattainable 
with the present techniques. Second, it seemed to 
us important to check the free electron method of 
calculating wave functions and heats of vaporiza- 
tion? in a case where experiments give an upper 
. limit to the energy. The obvious instability of 
metallic hydrogen under ordinary pressure shows 
that its heat of formation is smaller than that of 
the molecular hydrogen. This condition proved 
to be only too well fulfilled. Finally, we hoped to 
obtain some information concerning the condi- 
tions for the metallic state and the line separating 
metals from metalloids in the periodic table.’ 
This purpose has been achieved by the present 
work to the extent expected, but this point, 
although perhaps apparent to the reader, will not 
be stressed in the present paper. 


(2) CALCULATION 


The calculation itself was a rather mechanical 
application of the principles outlined elsewhere.” 
First, the energy E(r,) was found for which the 
solution of the Schrodinger equation with the 
ordinary Coulombic potential —e?/r had zero 
derivative at distance r,. This was done by a 
numerical method and, in addition, checked by 
an analytic approximation,‘ which gave a direct 


‘It was J. D. Bernal who first put forward the view that 
all substances go over under very high pressures into 
metallic or valence lattices. : 

* (I) E. Wigner and F. Seitz, Phys. Rev. 43, 804 (1933); 
(II) 46, 509, 1934; (III) E. Wigner, ibid. 46, 1002 (1934); 
(IV) F. Seitz, ibid. 47, 400 (1935); (V) J. C. Slater, ibid. 
45, 794 (1934). E. L. Hill remarks in a ‘‘Note on the Sta- 
tistics of Electron Interaction’’ (Physik. Zeits. Sowjetunion 
7, 447 (1935)) that the expressions ‘‘parallel spin” and 
‘antiparallel spin” should be replaced by “parallel spin Z 
components” and “antiparallel spin Z components,” re- 
spectively. The former expressions are not quite equivalent 
with the latter ones according to the general formalism of 
the statistical interpretation of quantum mechanics and 
can be used only as an abbreviation for them. (This is 
done in the present paper also.) We believe that Dr. Hill’s 
further remarks concerning the papers of this reference are 
based on a misunderstanding. 

* Interesting quantum-mechanical considerations on this 
point have been put forward by F. Hund (International 
Conference on Physics, October 1934, London). 

* To be given in Section 3. 
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expression for the wave function needed for 
subsequent work. The plot (lower curve of Fig. 1) 
shows the results of the numerical method which 
gave the more accurate results. Its values lie 
about 0.02 Ry (Rydberg units) below those of 
the analytic formula. In this case E(r,) has no 
minimum, as in case of Na, but instead continues 
to fall to —* for r,—0. In Na the rising of 
E(r,) for small r, results from the energy neces- 
sary to press the valence electron into the inner 
shells, whereas no such shell exists for hydrogen. 

The upper curve of Fig. 1 is obtained by adding 
to E(r,) the Fermi energy for free electrons, 


(2.2099/r.22)Ry, 


where 7, is expressed in Bohr units h?/e’m. This 
corresponds to the calculation adopted in ref. 2 I 
and gives a minimum of E = —1.156 for r,=1.50. 
This would give a density of 0.80 (as compared to 
0.087 in the solid molecular form!) and a heat of 
vaporization of 48.7 cal. into atomic hydrogen, as 
compared with 52.5 cal., the heat of formation 
of the molecular form. We see that the molecular 
form is the more stable, even according to this 
very approximate calculation, but if this were 
all, it would probably still be possible to trans- 
form it under pressure to a metallic form. 

Unfortunately, a great number of corrections 
must be applied to the simple picture given in 
ref. 2 I which decreases the stability of the metallic 
form in this instance. These corrections are given 
for the flat wave functions in refs. 2 II and 2 III. 
They are: 


0 














Fic. 1. Lower curve—the energy of the electrons without 
momentum, the upper curve—total energy as a function of 
lattice constant. Both according to a naive free electron 
theory. 











































766 E. 





(a) The difference between the energy of the electron 
clouds (§2 in II) and the exchange energy (§3 in II) 
0.284/r,—0.017/r.2(Ex.—E), (in Rydberg units), (1) 


where Ex is a mean value of the energy of the higher s 
terms and has been estimated‘ to be about 20 Ry/r,’. 
The second term in (1) is, therefore, almost negligible. 
(b) The correlation energy, calculated in ref. 2 III, 
which can be represented for 7,>1 rather closely by 
—0.584/(r.+5.1). (2) 


(c) A correction for the Madelung number will be added 
in the next section. 


The energy after the corrections (a) and (b) have 
been added is plotted on a magnified scale on the 
lowest curve in Fig. 3. One sees that a great part 
of the binding energy is lost. 


(3) 


It is very probable that the wave functions 
have nearly the form e‘?-7/4, not only in Na, but 
also in the other alkalies excepting Li.’ In H this 
is not the case, however, and there will be further 
corrections introduced into the energy curve 
similar to those calculated for Li by F. Seitz.” !V 
Still the situation for H is rather different from 
that for Li. In H the wave function of the 
lowest energy electron (p=0) is strongly modu- 
lated but the modulation does not increase much 
for higher p functions. In Li the wave function 
for p=0 was quite flat and the modulation 
increased with increasing p. 

For the calculation of these effects an ap- 
proximate method has been used. The equation 
for the wave function of the electron with the 
momentum #, multiplied by e‘?:7/, is® 


H,xu= —Ax,—tu- grad Xu —2Xy/7 = EvXus (3) 
u=2hp/e*m. 


Here the e, does not include the kinetic energy of 
the flat wave functions, p?/2m=y?/4 Ry. The 
€.—€o will give us the change in the Fermi 
energy, and the wave functions x, the change in 
charge distribution, which are both caused by 
the fact that the electrons are not free in the 
lattice. The charge distribution is used in a 


5 Cf. F. Mott and C. Zener, Proc. Camb. Phil. Soc. 
30, 249 (1934); also J. Bardeen and E. Wigner, Phys. Rev. 
48, 84 (1935); for Li see reference 2 (IV) and J. Millman, 
ibid. 47, 286 (1935). 

° Cf. Eq. (11), reference 2 (II). The energy in (3) is in 
Rydberg units Ry and the length in Bohr units. 
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recalculation of the interaction of the electron 
clouds and the Fermi hole (exchange energy). 
We shall see that x, can be assumed to be the 
sum of a function of 7 alone (s part) and a 
function of r multiplied by (u-7) (p part). The 
boundary condition for the x, is periodicity in 
the lattice, which means zero derivative for the 
s part, and vanishing of the p part on the surface 
of the s sphere. 

Eq. (3) has been solved by a perturbation 
method. We shall need only those solutions of the 
unperturbed equation 


— Agni =niéni, (4) 


which are functions of r alone (/=0, or the s part 
of x), and those which are functions of 7 multi- 
plied by (u-r), (J=1 or the p part of x). The 
boundary conditions for these £ functions are the 
same as for the corresponding part of x. We 
have? ™ 


Evo = (3/42r,*)}, 
Eni = (34cn/uwn)(u grad) (sin war)/r, 


Eno =Cn(SiIN Wpr)/?, 


wn 
waa 








Noo= 0, n0 = ni = Wn”, 
1.4307 1.024 
with O1= ’ C1 = > 
Ts (2xr,)! 
T 1 
On = (n+3)-, CS 
Ys (2xr,)? 


when the last two expressions hold well for n>1 
only. 

Next we must calculate the matrix elements of 
H,. All terms give nonvanishing elements be- 
tween £,; and é,-; only, except for iu- grad which 
gives elements between £,; and £&n/14;. We have 


(77 ,) 00: 00 = —3/?.; 


(7,)10: 10= 20.2/1r,?—5.65/r,; 

(77 ,) no: no (Hy)ni; ni * ("+ 3)?4?/1,?; 
(77,)11; 11.= 20.2/r,? — 3.66/1,; (6) 
(Hy) n’0; n1= — (tuin/3*) barn 

(IZ,)00; 10= 1.36/r.; 

(H,)00; no=2(6)/(m+ 3) ar. 
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For »=0, this gives, by Rayleigh-Schrédinger 


formulas, 


(1.36)? « 


’ 
w= (n-+4)'4 


r, 20.2—2.65r, 


1.367, oo 


xo a 
20.2 —2.65r, 


2(6)'r, 
£10 ae 
n=2(n-+3)*a° 
These are the formulae to which we referred in 
Section 2. For 10, we diagonalized the three 
rows and columns of H, referring to £00, £19 and £11 
and used only for the remaining ones the 
Rayleigh-Schrédinger method. This was ad- 
visable, because (/7,)10; 10 and (/7,) 11; 11 were close 
together. This gave 
3 1.85 


-#, 20.2 —2.657,—6.72u2(20.2 —0.6r,)—! 
—0.0093. 





Eu 


(8) 


The mean value of ¢,—e€9 with the weight yp’, 
from »=0 to p=3.84/r, is the correction to the 
Fermi energy; it is plotted in the lowest curve of 
Fig. 2. It is negative, just as in Li, because there 
is no p level below the s level of the valence 
electron. The energy ¢,—€o is no longer very 
accurately proportional to yw’, because the effect 
of the next Brillouin zone is appreciable. 


(4) 

Next the interaction of the electron cloud and 
the more accurate exchange energy must be 
derived. The greater density near the nucleus 
makes it evident that the first will be greater 
than the free electron value assumed in the 
previous section. This is true also for the ex- 
change energy, because the Fermi hole is more 
effective in consequence of the uneven charge 
distribution. The two effects work against each 
other, but in contrast to Li, the former is greater 
here and decreases the binding. 

The calculation of these quantities is a little 
laborious and contains nothing of principal 
interest. It will be sketched only. 

The square of the calculated wave function, 
i.e., the electric density p, was approximated by 
an analytical expression’ containing a constant 


’ The reason for this is that only for this form were we 
able to perform the integrations. 




















Fic. 2. Bottom line, correction to the Fermi energy; 
middle line, correction for the Madelung number; top line, 
correction for the interaction energy. 


and an exponentially decreasing charge distri- 
bution around each ion: 


p=pitpe, pi=e(l1—c)/(4rr,3/3). (9) 


p2 itself contains two parts: first the exponential 
charge distribution around the ion, in the s 
sphere in which we are, and second the exponen- 
tial charge distribution around the other ions. 
The latter was averaged over all directions so 
that p2 also became spherically symmetric in each 
s sphere 


ecp* Sh(6r) 
ma ——| eo 44h | 
8r Br 


9 
k=e~*6rs(1 +-8r,+6?r,?/2)(Br,—1)-. 92) 


The two constants c and 8 were adjusted so as to 
fit the square of the calculated wave function (7) 
best. 

The potential corresponding to this charge 
distribution contains also two parts: #,, corre- 
sponding to the constant part p; and 9% corre- 
sponding to p2 which could be easily calculated, 
and hence also the interaction of the electron 
clouds. 
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For the calculation of the exchange energy the 
wave functions of the electrons with momentum 
p were assumed to be in the form e#?-*)/hy(r) 
where x(r) is the wave function of the lowest 
state. Under these conditions, the exchange 
energy is? 1 21 


e 


1 
—5 J etryete"ye(| nr" |) ara, 


r—r'| 


(10) 


where 





sin (r/d) —(r/d) cos (r/d) 7? 
g(r) -4 


(r/d)* 


r r\ 2 
~e-wid 141.6(-) 41.2(-) | (10a) 
d d 


with d=0.521r,. All the integrations could be 
carried out in elliptic coordinates. The result is 
evidently a quadratic function in c. The constant 
term is simply the exchange energy for flat wave 
functions (c=0), i., 0.916Ry/r,. The linear 
terms vanish as may be seen from the following. 
The integral of the first term in p,(r) times 
g(|r—r’|)/|r—r’| over r is independent of 7’. 
The product of this quantity and the terms 
containing c in p(r’), integrated over 7’, must 
vanish, since the integral of p(r’) is independent 
of c. By analogy the other linear term vanishes 
also and the sole remaining correction to the 
exchange energy of 0.916Ry/r, goes with c’. 
Hence this is a second order effect* of the 
deviation from flatness of the wave function (9), 
and smaller than the increase of the interaction 
of the charge clouds which is a first order effect. 

Next the change of the correlation energy, due 
to the uneven charge distribution, has to be 
considered. In the case of even distribution, the 
correlation energy is certainly smaller than 
0.292e?/r,, even for a very great lattice constant. 
For the actual uneven distribution in a greatly 
expanded lattice, however, this energy is as 
follows. Hartree’s equations, solved for many 
electrons in the field of many ions, give as wave 
functions of the electrons a sum of spherically 
symmetric maxima around every ion. If all ions 
are sufficiently far removed from each other, 


8 This differs from the case for Li. 
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these waves—let their form be u(r)—no longer 
overlap. The total interaction energy per electron 
is 


(11) 


fo | u(r) | *4arr?dr 


where ®(r) is the potential of the charge distri- 
bution | u(r) |*. Since no such interactions actually 
exist, Hartree’s equation gives an incorrect 
result by this amount, which will be compensated 
by the exchange and correlation energy. 

If Fock’s picture is used instead, the interaction 
energy between electrons with parallel spins 
vanishes in consequence of the enormously 
increased Fermi hole. There remains still the 
interaction between electrons with antiparallel 
spin, by which amount Fock’s picture is still in 
error. This remaining half of (11) must be 
nullified by the correlation energy. It follows 
that this energy must greatly exceed 0.292e?/r,, 
the largest value possible for flat functions. We 
see that the correlation energy is increased by the 
modulation of the waves by a greater factor 
than the exchange energy, and it becomes in the 
limiting case of an infinite lattice just one-half of 
(11), approximately }Ry per electron.® 

We have increased the correlation energy for 
flat waves, given in,?"™! by the same factor by 
which the exchange energy appeared to be 
increased. The total change of the electron 
interaction energy was positive, and plotted in 
Fig. 2 also. 

Finally, the correction for the Madelung 
number, as given in appendix 2, Reference*™ is 
given in the same plot. The sum of all three 
corrections has been added to the curve of Fig. 3 
and the upper full line of this figure obtained. 
This curve assumes its minimum of —1.05Ry 
(heat of vaporization about 16 Cal.) for r,=1.63 
Bohr units, corresponding to a density of 0.62. 
The possible error in these figures is quite large, 
because the correction quantities are so great. 


(5) 


The last quantity which had to be calculated 
for the thermal constants of metallic H was the 
zero-point energy corresponding to the _half- 


« * We emphasize this point since doubts have been raised 
whether the correlation energy actually exists. 
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Fic. 3. Energy of the lattice as a function of lattice 
constant. Lowest curve, for flat wave functions; second 
curve with all corrections, except zero-point energy of 
nuclei; the dotted lines contain the zero-point energy, the 
lower for heavy, the upper for ordinary H. 


quanta vibration of the nuclei. If the elastic 
constants of the medium were known, the elastic 
spectrum and the zero-point energy could be 
calculated following Blackman” or a paper of 
Barnes, Brattain, and Seitz." Since, however, 
only the compressibility can be derived from 
Fig. 3 (and this only very inaccurately, as it may 
be in error by as great a factor as two), it was 
performed by the older method of Debye.” 

From the curvature of the upper curve in Fig. 
3 one obtains for the reciprocal compressibility 
of a lattice of very heavy hydrogen, 


(C11 +2¢12) = 2.9K 10" dyne/cm? for r,=1.5 
= 1.510" dyne/cm? for r,= 1.6 
= 1.110" dyne/cm? for r,=1.7 
= 0.8 X10” dyne/cm? for r,=1.8. 


By analogy with other metals ¢1; = 2¢12=4c44 was 
assumed. The zero-point energy is then about 


_ 1° M. Blackman, Proc. Roy. Soc. London A148, 365, 384, 
(1935); 149, 117, 126 (1935). 
Barnes, Brattain and Seitz, Phys. Rev. 48, 582 (1935). 
" Cf. e.g., E. Schrédinger, Geiger-Scheel Handbuch der 
Physik, Vol. X. 
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0.0244(E’’)!Ry, where E”’ is the second derivative 
of the energy in Ry with regard to 7, in Bohr 
units. For the isotropic hydrogen with mass two 
the value is 0.0173(£”’)'Ry. The dotted lines of 
Fig. 3 give the final results for the lattice energy, 
the upper one is for light, the lower for heavy 
hydrogen. The density of the former is 0.59, its 
heat of formation out of atomic hydrogen 10 Cal. 
We think that this figure is somewhat too low 
because the correlation energy is greater than 
we assumed it to be. 


(6) CONCLUSION 


The great density of the hypothetical metallic 
form of H is not surprising. It seems to be the 
general rule that the density of the metallic and 
valence lattice of a substance be greater than 
that of the layer lattice, and this in turn has 
greater density than the molecular lattice. Table 
I illustrates the point. 

This table makes it natural to assume that all 
atoms are bound together by shared electron 
forces in metallic and valence lattices, and that 
this is the case also for the atoms within a layer in 
laver lattices. The different layers are bound 
together, however, more weakly, perhaps by 
van der Waals forces. This is well in accord with 
the general geometric considerations of A. Reis 
and K. Weissenberg. It seems sensible, there- 
fore, to speak about valence-like forces between 
two atoms in a lattice, if the density of the 
outermost electrons of both overlaps considerably 
at this distance." 

The heat of formation of the ordinary, molecu- 
lar He lattice is 52.5 Cal., practically one-half of 
the heat of dissociation. Although the figures of 
the preceding section may be in error by quite a 
few calories, it is evident that the molecular form 
is so much more stable (by about 40 Cal.) than 
the metallic, that the chances of obtaining the 
latter are extremely small. 

In the ordinary lattice, there are two ‘‘nearest 
distances.”’ The distance between two atoms of 
the same molecule is 0.75A=1.45a9, while the 
nearest distance between two H atoms of 


13 A. Reis, Zeits. f. Physik 1, 204 (1920); 2, 57 (1920); 
K. Weissenberg, Zeits. f. Krist. 62, 12, 52 (1925). 

14 This corresponds also to the ideas developed by L. 
Pauling, J. Am. Chem. Soc. 53, 1367 (1931); J. C. Slater, 
Phys. Rev. 38, 1109 (1931). 













































TABLE I. 

SUBSTANCE DENSITY DENSITY 
As, metallic 5.72 yellow 2.03 
Diamond 3.51 graphite 2.24 
Black phosphorus 2.70 yellow 1.83 
Se, metallic 4.82 red 4.47 
Sn, white 7.28 grey 5.76 








different Hz molecules is about 3.3A=6.3ao." If 
one compresses the ordinary lattice, the latter 
distance will decrease, the former increase, as 
follows already from the calculations of F. 
London." Thus, the increase of the pressure 
brings the intermolecular distances nearer to the 
intramolecular distances, and the lattice will 
become more similar to a simple lattice than it is 
at low pressures. It is possible that this goes on 
continually, and that one reaches a simple close 
packed lattice with metallic properties only in 
the limiting case of infinite pressures. 

During the increase of pressure, the total 
energy of the substance will increase by the work 
done (we assume very low temperatures). It is 
possible, thus, that the energy of the continu- 
ously compressed lattice becomes at a certain 
volume, assumed at a very high pressure, greater 
than that of another modification. In this case 
there will be a discontinuous change to the latter."” 
Whether or not such a transition will occur into 
the metallic modification, depends on the magni- 
tude of the compressibility of the molecular form 
at high pressures. The smaller this is, the greater 
will be the pressure necessary for a certain 
change in volume, and also the energy increase 
at this compression. If the extremely high 


16 Sy calculated nearest distance for the metallic form 
is 1.5A. 
1 F, London, Zeits. f. Elektrochemie 35, 552 (1929); 
H. Eyring, J. Am. Chem. Soc. 53, 2537 (1931). 
17 Such polymorphic transitions induced by pressure in 
ry ee by P. W. Bridgman, Rev. Mod. Phys. 
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compressibility at ordinary pressures,'* which is 
2.10-° cm?/dyne, would hold throughout, the 
molecular form would be stable for all volumes. 
This is not possible, however, and the com- 
pressibility certainly decreases with increasing 
pressure, but one cannot tell, at present, whether 
it decreases sufficiently to make the metallic 
form stable at any pressure. One calculates easily, 
that even under the assumption of the most 
advantageous compressibility at high pressures, 
the pressure necessary for the transformation is 
250,000 atmos., which is outside the scope of the 
present technique. 

The objection comes up naturally that we have 
calculated the energy of a body-centered metallic 
lattice only, and that another metallic lattice 
may be much more stable. We feel that this 
objection is justified. Of course it is not to be 
expected that another simple lattice, like the 
face-centered one, have a much lower energy,— 
the energy differences between these forms are 
always very small. It is possible, however, that a 
layer-like lattice has a much greater heat of 
formation, and is obtainable under high pres- 
sure. This is suggested by the fact that in most 
cases of Table I of allotropic modifications, one 
of the lattices is layer-like,!® the other either 
metallic, valence, or molecular. The difficulties 
for such an experiment will be greatly increased 
by the necessity of the formation of a nucleus 
for the new lattice. 


18 Measured by A. Eucken for the liquid form, Ber. 
deutschen Phys. Ges. 4 (1916). 

19 Diamond is a valence lattice, but graphite is a layer 
lattice (A. W. Hull, Phys. Rev. 10, 661 (1917)). Yellow 
arsenic and phosphorus are evidently molecular lattices, 
black phosphorus a layer lattice (R. Hultgren and B. E. 
Warren, Phys. Rev. 47, 808 (1935)), and metallic As is 
also only approximately a simple lattice. The red, mono- 
clinic selenium can be said to form a molecular lattice (F. 
Halla, F. Bosch, E. Mehl, Zeits. f. physik. Chemie B11, 455 
(1930)), while the metallic modification is a thread lattice 
(A. J. Bradley, Phil. Mag. 48, 477 (1924)). The situation 
seems to be most complicated with tin. Grey tin forms a 
diamond lattice, but shows otherwise no similarity to the 
valence lattice of diamond, while the metallic lattice has a 
rather complicated structure (H. Mark and M. Polanyi, 
Zeits. f. Physik 18, 75 (1923)). 
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The electron affinity of iodine atoms has been determined by a direct method in which only 
iodine atoms and electrons were involved. The value so obtained is 74.6 k-calories, which is in 
good agreement with values obtained by other methods. The concentration of the several ions of 
different masses (I~ and E~) have been calculated from their effect on space charge. It has been 
established that iodine has no effect on the thermionic emission of tungsten. 





I. INTRODUCTION 


ANY attempts have been made to de- 

termine directly the electron affinity of 
the elements and to this date comparatively few 
values have been obtained.* One of us‘ has 
described an empirical method of estimating the 
electron affinity of the light elements and it is 
our plan to check these empirical values by 
experiments in this laboratory. In the present 
paper we report an experimental determination 
of the electron affinity of iodine atoms obtained 
by a new method. We have determined the 
equilibrium of reaction 


X+E-—X-+AE 


at the surface of a hot tungsten filament and 
from the equilibrium constant we have succeeded 
in calculating the electron affinity: —AZo. As far 
as this fundamental principle is concerned it 
appears that Sutton and Mayer® have employed 
the same idea and for the same purpose. Since, 
however, the measurement of the important 
ratio of the ion to the electron current (I~ to E~) 
in our work differs radically from the method 
used by these authors, we feel that we should 
present our results. Iodine was used in our study 
for the reason that its electron affinity is well 
known and hence any new experimental method 
can be tested by the use of this element. 


1 Presented at the New York Meeting of the A. C. S. 
April, 1935. 

* This article is based upon a thesis presented to the 
Faculty of the Graduate School of the University of Min- 
nesota by Melvin Calvin in partial fulfilment of the require- 
ments for the degree of Doctor of Philosophy. 

*O. Oldenberg, Phys. Rev. 43, 534 (1933). This paper 
contains further references. 

* Geo. Glockler, Phys. Rev. 45, 111 (1934). 

* Sutton and Mayer, J. Chem. Phys. 3, 20 (1935). 


II. THEORETICAL CONSIDERATIONS 


The equilibrium constant 
In the case of the ionization of caesium on 
tungsten Langmuir® was able to measure the 
equilibrium 
Cs-Cst+E-—AU 54 (1) 


established at the surface of a hot filament, 
where AU), is the ionization energy derived from 
spectroscopic data. He obtained very good 
agreement with the theoretical value of the 
equilibrium constant calculated by the equation 


n(E-)Xn(Cst) w(E-)w(Cs*) 
n(Cs) _ w(Cs) 
(2am(E-)kT)? 
x as x 
1 


K(Cs, n) = 








eAU0/RT. (2) 


The n(x) are the concentration of the species 
“‘y’’ in numbers of particles per cc, the w(x) are 
the statistical weights of the normal states, m(x) 
are the weights in grams of the individual 
particles x and the other symbols have their 
usual meanings. If in our case the reaction 
occurring at the filament surface is 


I--I+ E-—AE (3) 
then 
K(I-, n) =n(I) Xn(E-)/n(I-). (4) 


Since the electron affinity is defined for the 
reverse reaction, we shall write 


R(I-, n) =n(I-)/n(1) Xn(E-) =1/K(I-, n) (5) 
and from Eq. (2) we have 
R(I-, n) =[w(E-) Xw(I)/w(I-) 
X (24m(E-)RT)3/h3 Xe-420/RT 1, (6) 


6 Langmuir, Phys. Rev. 43, 226 (1933). 
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If, then, k(I-, ~) can be determined experi- 
mentally, the electron affinity of iodine can be 
obtained from Eq. (6). The constant k(I-, 7) 
can be obtained from Eq. (5) if the concentrations 
of the interacting species can be determined. 
n(I) can be obtained from the vapor pressure of 
iodine assuming that every molecule striking the 
filament surface is dissociated into two atoms® 
and the ratio n(I-)/n(E-) can be obtained from 
the space-charge equation. Jf is in the determi- 
nation of this ratio that our method differs from 
the work of Sutton and Mayer.’ While they 
employed a magnetron set-up to differentiate 
between the ion and electron current, we 
employed the current-potential relations of a 
space-charge region to evaluate them. 


The space-charge equation 


The space-charge limited current between a 
pair of electrodes is given by the well-known 
Langmuir-Child equation’ 


log i=log C(e/m)!+# log V, (7) 


where C is a geometrical constant, e/m is the 
ratio of charge to mass of the carriers of the 
current and V is the potential difference between 
the electrodes. By taking current-voltage data 
in vacuum where electrons only carry the cur- 
rent and plotting log i against $ log V we can 
determine the geometrical constant C for any 
experimental conditions. By plotting similar data 
obtained with iodine present between the elec- 
trodes an effective value of e/m can be calculated 
which will be different from that of the electrons 
and will depend on the number of other types of 
charged ions present between the electrodes. 

In order to determine the relationship between 
the effective e/m of Eq. (7) and the number of 
heavy negative ions (I~) present we must re- 
consider the derivation of this equation. It is 
obtained from Poisson’s equation 


V?V=4rp, (8) 
and for cylindrical electrodes: 
2arp=1/u=i(m/2eV)}, (9) 


where p is the negative charge density and u is the 


70. W. Richardson, Emission of Electricity from Hot 
Bodies (Longmans, Green & Co., London, 1921). 
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velocity of the charged particle: ((2eV)/m)'. 
Eqs. (8) and (9) yield® 


rd2V /dr?+dV /dr=(2m/eV)%, (10) 


which gives on integration from r=0 to r=r, 
(dV /dr)=0 and V=0 at r=0: 
i= 2v2/9r(e/m)'V? 


or const X Vi=i(m/e)}. 


(11) 


In the special case we are considering we shall 
assume that there are present ions of two masses: 
m(E-) and m(I-). There will be then two 
currents: 7(E~) and i(I-) and 


p=p(E~)+p(I-) =(E-)/u(E-) +1(1-)/u(T-), 
or 
p=(1/v2){i(E~)(m(E-)/e)} 
+i(I-)(m(I-)/e)}} V4. (12) 


The form of this expression is the same 
as Eq. (9) with i(£-)(m(E-)/e)! replaced by 
{i(E-)(m(E-) /e)'+2(I-)(m(I-)/e)4}. Since these 
terms have the same functional relation to V and 
are subject to the same boundary conditions, the 
integration will be exactly analogous and we can 
write the solution directly 


const X V!=i(E-)(m(E-)/e)} 
+i(I-)(m(I-)/e)*. (13) 


Hence the effective value of m/e as determined 
from experiment is related to its constituents as 
follows: 


i(m/e)'=i(I-)(m(E-)/e)!+i(I-)(m(I-)/e)', (14) 

i=i(E-)+i(1) (15) 
i(I-) _m'—(m(E-))'__w(I-)Xe 
i(E-) (m(I-))'—mt (E>) Xe 


and 


F (16) 





and 


where u(I-) is the number of I- crossing unit 
area in unit time between filament and plate and 
u(E-) is the corresponding quanity for electrons. 


The equilbrium equation 


It is well known from kinetic theory that the 
concentration of a species is related to the number 
traversing unit area in unit time, by 


8’ Langmuir and Compton, Rev. Mod. Phys. 3, 191 
(1931); and Langmuir, Phys. Rev. 2, 450 (1913). 
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u(x) =n(x)(kT/2rm(x))?. (17) 


The mass of I~ is sensibly the mass of I atom and 
since the three particles (I, I~, E~) are formed at 
the filament their temperature is the temperature 
of the filament: 7,. The number of iodine atoms 
formed at the filament is twice the number of I, 
molecules striking it: 


u(I2) = P(I2)/(4rm(T)kTe)', (18) 


where P(I) is the pressure of iodine molecules in 
baryes and J¢ is the temperature of the iodine 
gas (about room temperature). The equilibrium 
constant (Eq. (5)) can now be written: 


3X10-" 
—- (19) 


k(I-, ) = 





Fk m(1)Te@} 
«| | x 
P(I,) L2m(E-) 


and combined with Eq. (6) permits the calcula- 
tion of the electron affinity. 


III. EXPERIMENTAL PROCEDURE 


The experimental set-up is shown in Figs. 1 
and 2 and consists of a diode (filament and plate) 
with an iodine reservoir. The filament is a 
tungsten wire 0.1 mm thick and about 42 mm 
long. The plate is made of tantalum (0.152 mm 
thick) and is in the form of a cylinder 45 mm 
high and 21.4 mm diameter. The glass tube is 
35.6 mm in diameter. The details of the electrical 
connections are shown in Fig. 1. The experi- 
mental set-up is essentially that used by Heinze® 
in measuring the work function by the cooling 
method. It is a Wheatstone bridge with audio- 
frequency current for balancing and direct 
current for heating the filament which forms one 
of the arms of the bridge. The various parts of the 
electrical arrangement are: 


Am, ammeter measuring approximately the current 
through the bridge. 

A,B, endcoils of Leeds and Northrup Kohlrausch bridge 
each having a resistance of 124.14 ohms. 
blocking condenser having a capacitance of about 
one mf. 
variable condenser having a range of 0.00003-0.04 
mf. It is necessary to balance out any capacitance 
that might be introduced into the bridge by the 
experimental tube or at other points. 
experimental tube. 


* Heinze, Ann. d Physik 16, 41 (1933). 





























Fic. 1. Reaction tube. 


parallel resonance circuit tuned to the frequency of 
the oscillator and thus providing an infinite im- 
pedance to the audiofrequency attempting to get 
through the plate circuit and preventing this 
auxiliary current from falsifying the bridge measure- 
ments. 

the filament of the experimental tube. 

manganin slide wire in mercury (10 ohms). 

slide wire of Leeds and Northrup Kohlrausch bridge 
(type 54161). Its resistance was checked on a 
standard bridge of the Physics Department and is 
27.72 ohms. 

lead storage “‘B”’ batteries. 

set of four two hundred-ampere hour lead storage 
cells. They could be arranged in any series or 
parallel combination so as to deliver up to 24 volts. 
It is very essential that these batteries be capable 
of maintaining a current of about 1.5 amperes 
steady to within 0.005 ampere over a period of 
one hour. 

manganin slide wire in mercury (one ohm). 


., Microammeter with shunts capable of measuring 


currents from 210-7—3X10- amp. It measures 
the plate current. 
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Fic. 2. Electrical connections. 


vacuum tube oscillator which supplies the audio- 
frequency current to activate the bridge. 
potentiometer (2700 ohms). 

potentiometer (10,000 ohms) used as a Wagner 
ground on the bridge. 

type ‘“‘K”’ Leeds and Northrup potentiometer across 
the standard resistance ‘‘S.” It is used to make the 
critical measurement of the current through the 
filament. 

the plate of the experimental tube. 

variable resistance (100 ohms) for controlling the 
current through the filament. 

variable resistance (25 ohms). 

standardized three-dial noninductive resistance 
(999 ohms). It is necessary in order to achieve a 
balance when the end coils are in. 

one-ohm Leeds and Northrup standard resistance. 
reversing switch for reversing the current through 
the bridge. 

arrangement for connecting the negative terminal 
of the ‘‘B”’ battery to either end of the filament. 
three-tube amplifier with telephone used as the 
balance detector of the bridge. 

voltmeter for plate voltage. 


The tube was outgassed by standard procedure 
and the work function of tungsten was de- 
termined by two methods. 
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The work function of tungsten 


The work function of the filament material 
(tungsten) was obtained from the temperature 
variation of the saturation current by the use of 
well-known relations. The object was to ascertain 
that the filament surface was clean. The value 
for the work function was found to be 4.59 
electron volts which is in very close agreement 
with the accepted value for pure tungsten. The 
second method employed was a measurement of 
the cooling effect on the filament due to the 
evaporation of the electrons. The value obtained 
is 4.60 electron volts in excellent agreement with 
the value obtained from the Richardson equa- 
tion. The temperature of the filament was 
obtained from the heating current and the 
Langmuir table." 

These measurements established the following 
facts: (1) The apparatus will actually measure 
heat values at the surface of the filament. (2) The 
methods used in the measurements of the various 
experimental quantities and the estimated cor- 
rections (for example end losses etc.) are valid. 
(3) The material of the filament is satisfactory, 
i.e., clean. (4) The tube is free from residual gases. 
Admission of iodine 

The tube having been prepared by outgassing 
and filament-aging was then filled with iodine 
vapor. After the iodine was admitted the filament 
was again allowed to age for at least thirty hours, 
always keeping liquid air on the tube containing 
the solid iodine supply. At this point the filament 
was subjected to the emission test and it was 
found to be unaltered. Moreover raising the 
iodine pressure even as high as 10-* mm Hg had 
no effect on the saturation emission charac- 
teristics of the filament. This seems entirely 
reasonable for before a molecule can affect the 
thermionic emission from a surface it must be 
adsorbed on the surface. At our high filament 
temperatures (above 2000°K) iodine atoms are 
not adsorbed. This finding is in agreement with 
the work of Van Pragh'! who studied the 
accommodation coefficient of helium on platinum 
which had been subjected to various treatments 
in iodine. He found that platinum treated in 


10 Langmuir and Jones, Gen. Elec. Rev. 30, 310—9 (1927). 
1! Van Pragh, J. Chem. Soc. 798 (1933). 
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iodine above 1650°K showed an accommodation 
coefficient toward helium corresponding to that 
of a clean surface. But when iodine was allowed 
to come into contact with the platinum at any 
lower temperature, the accommodation coeffi- 
cient changed immediately to that of a platinum 
surface covered with a layer of iodine. Since 
tungsten iodide is much less stable than platinum 
iodide it seems reasonable to say that at our high 
temperatures todine has no effect on the thermionic 
emission from tungsten. This is in agreement with 
statements made by Langmuir in private com- 
munications to others.” 

But it still remains true that the presence of the 
iodine vapor between the electrodes can and does 
influence the plate current. In the space-charge 
limited region and below the ionization potential 
of iodine (molecule or atom) the current is 
affected by the creation of heavy negative ions 
(I- etc.) which help to increase the negative 
space-charge density and thus decrease the 
electron current. As soon as the ionization 
potential is reached positive ions are formed, the 
negative space charge is neutralized and the plate 
current jumps above its value in a vacuum. If 


the number of positive ions is large enough the 
saturation current may even be larger than that 
in a vacuum. 


Current-voltage relation in iodine 


The data needed for the calculation of the 
electron affinity are simply a set of current- 
voltage tables taken in the space-charge limited 
region and at voltages less than the ionization 
potential of iodine. This information together 
with the temperature of the filament, the pressure 
of the iodine at which the measurement was 
taken and the gas temperature is all that is 
needed for the calculation of the electron affinity. 

The filament temperature is first set with 
liquid air on the iodine supply. The Wheatstone 
bridge in which the filament is located is balanced 
and the plate current is measured as a function of 
plate voltage. The iodine is then admitted to any 
desired pressure by warming the side tube in the 
proper bath. The filament will cool off because of 
the thermal conduction of the iodine vapor. 


"Van Pragh and Rideal, Proc. Roy. Soc. A134, 385 
(1931), and reference 5. 


IODINE 





P(I,), mm Hgs!0* 
Vac 

é 

26 


92 


LOGi +6 —— 


Temp.= 2332°K 











%LOGV — 


Fic. 3. Current-voltage relations in iodine. 


This cooling effect causes a drop in resistance of 
the filament which will become apparent on the 
Wheatstone bridge. By increasing the filament 
heating current so as to reestablish the filament 
resistance originally found in a vacuum, we assume 
that the filament temperature is returned to its 
original value. Again we find the plate-current as 
a function of plate voltage under the new 
conditions and repeat the procedure for several 
pressures of iodine. For each filament tempera- 
ture a new set of vacuum data must be taken. A 
sample run is shown in Table I. The filament 
temperature is 2332°K. The data are plotted in 
Fig. 3. 


IV. DiscuUSSION OF RESULTS 


The current ratio F 


The data given in Table I and shown in Fig. 3 
are used to calculate effective values of m which 
satisfy the Langmuir-Child equation under any 
given set of conditions and hence the value of F 
the ratio of the current carried by the iodide ions 
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TABLE I. Current-voltage relation in I>, filament temperature 
2332°K. 
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(I-) and the electrons (E~). From Eq. (7) we see 
that the intercept on the log i-axis of Fig. 3 
gives the constant C(e/m)!. The straight line 
referring to the conditions in a vacuum shows of 
course the relation when electrons only are 
carrying the current. Hence e/m is that of the elec- 
tron (4.774X10-" 6.06 x 107 x 1845 x 1.008—! 
=7.27X10°). For the case shown in Fig. 3 we 
find for a vacuum 


log C(e/m)!=0.66 — 6.00 =4.66 — 10, 
C(e/m)!=4.57 x 10-® 
hence C=4.57X10-°/7.27 X 10®=6.28 X10-" 


and 


in units where electric charge is in electrostatic 
units and masses in weights of atoms in grams. 
With iodine in the tube at a pressure of 9.2 x 10-* 
mm Hg we find the straight line log z vs. volts to 
intersect the log 7 axis at a point 


log C(e/m)'=0.33 —6=4.33—10 
or C(e/m)'= 2.14 10-. 
Using the value of C found above we obtain 
(e/m)'=2.14 K 10-°/6.28 x 10-5 = 3.41 X 108. 
Considering the carriers singly charged we use 


e'= (4.774 X10-")}=2.19 x 10-5 
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2.19X10-° 


——— = 6.43 X10-™. 
3.41 X 10° 


and e!/(e/m)'=m'= 


We need to know 


(m(E-))!=[1.008/(1845 x 6.06 X 10”) }! 
310-4 
and (m(I-))!=[127/(6.06 x 10?) }! 
=1.45x10-" 


in order to evaluate F from Eq. (16): 
6.43 X 10-4 — 3.0K 10-"4 
~ 1.45 10-"—6.43 x 10-4 





= 2.37 X 10-3. 


This procedure is repeated for each iodine 
pressure, a new geometrical constant being de- 
termined at each filament temperature. 


Calculation of electron affinity 


Eqs. (6) and (19) can be used to calculate the 
electron affinity. The following values of con- 
stants appearing in the equation were used: The 
statistical weight of the electron w(E-) is 2; 
the normal state of the iodine atom is ?P 3/2 and 
hence w(I) =4 and for the iodine ion the ground 
state is Sy giving w(I-)=1; k=1.371XK10-"; 
m(E-)=9X 00-8; h=6.55X10-?7; m(I) = 20.94 
xX 10-8; Tg =300°K; e=2.71828. The final equa- 
tion becomes 


— AE =4.575 (log F+2 log T, 


— log P(I2)+1.075). (20) 


The results of our experiments are summarized in 
Table II. The measurements extend over a 


TABLE II. Electron affinity of iodine. 
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average: 74.6 kcal. 
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temperature range of 380°K and involve a 166- 
fold change in iodine pressure. 


Concluding remarks 


In the development of the relation between 
the effective mass of the carriers of electricity 
(m) and the ratio (F) of ion to electron current we 
assumed that there were present only two kinds 
of ions: I- and E~. The possibility of the forma- 
tion of heavier negative ions I;~ and I;~ im- 
mediately arises. Should these ions exist in the 
tube in quantity, then the F values would be two 
large and hence result in high values of —AEp. 
Since the average value obtained for the electron 
affinity is apparently not too high, we conclude 
that our original assumption is correct. This is 
further supported by the work of Hogness and 
Harkness*® who have shown that in the low 
voltage arc in iodine, the only result of a primary 
process is I~ and that I~ is the result of the only 
primary process. They have shown that I,~ and 
I; exist only in very small quantity and that 
they arise by secondary processes after I- has 
been formed. 

Another assumption that was made in ob- 
taining the equilibrium equation, is that every 
iodine molecule that strikes the surface of the 
filament is dissociated. This is also justified by 
the results obtained. Langmuir in a private 
communication to others® has stated that such is 
the case. It should also be noted that the heat of 
dissociation of iodine is only 1.5 electron volts. 


13 Hogness and Harkness, Phys. Rev. 32, 784 (1928). 


OF IODINE 

The work of Stark and Bodenstein™ indicates 
that iodine vapor at low pressure would be almost 
completely dissociated at 1200°C. The filaments 
used were all above 2000°K and since it is 
reasonable to assume that the filaments used 
are not perfectly smooth so that each molecule 
has a very good chance to approach temperature 
equilibrium with the filament, it is to be expected 
that almost every molecule striking the surface 
will be dissociated. There is only one piece of 
evidence to the contrary. Van Pragh and 
Rideal’® in their experiments on the attack of 
platinum and tungsten by the halogen vapors 
come to the conclusion that the dissociation even 
at our high temperatures is given by 


u(I) =2u(I2) Xe~?/*?, 


where D = dissociation energy of iodine molecules. 
No explanation for this discrepancy has been 
suggested. In our work it may be assumed that 
the results indicate complete dissociation of 
every iodine molecule striking the hot filament. 

The possibility exists that the electrons suffer 
reflection at the filament surface. This situation 
has already been discussed by Sutton and 
Mayer.® Our experiments are not accurate enough 
to shed light on this point. From the results 
obtained it may be seen that the various as- 
sumptions made must be reasonably correct, 
since the value of electron affinity found checks 
the results of other methods satisfactorily. 


(21) 


14 —_ and Bodenstein, Zeits. f. Elecktrochemie 16, 961 
(1910). 

us i Pragh and Rideal, Proc. Roy. Soc. A134, 385 
(1931). 
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The experimental and theoretical absolute rates for the nonadiabatic decomposition of N,O 


are shown to be in good agreement. The amount of chemical inertia present in other reactions 
involving the singlet-triplet transition of oxygen is considered. A convenient method of con- 
structing potential functions for polyatomic molecules which fit the spectroscopic data, and 
which reduce in the proper way for. the various dissociation processes, is indicated, and is 


carried through for the N2O molecule. 








HE method of calculating absolute reaction 
rates by means of statistical mechanics 
as developed at the present time! involves in- 
formation obtainable from potential surfaces. 
For a large number of reactions interest is 
confined to the lowest surface. There are also 
cases of importance for which, when magnetic 
effects are neglected, two potential surfaces 
intersect. When these effects are considered the 
surfaces, instead of crossing, approach each 
other very closely. In such cases the reacting 
system may pass from one of the noncrossing 
surfaces to the other in the neighborhood of the 
point of closest approach instead of proceeding 
across the potential barrier, and thus may fail to 
react. The fraction of the systems reaching the 
activated state which react we take equal to x. 
The transmission coefficient, x, plays exactly the 
role of c in Eq. (10) of reference 1. 
Nonadiabatic reactions have been considered 
by a number of authors.? As an example of such 
a reaction the dissociation of nitrous oxide into 
a nitrogen molecule and an oxygen atom will be 
considered. That the N.O reaction proceeds 
abnormally slowly for a unimolecular reaction 
is well known from the work of Volmer and his 
collaborators.* Wigner and others‘ have remarked 
that the reaction might be expected to proceed 


1 Eyring, J. Chem. Phys. 3, 107 (1935). 

2a. Wigner, Zeits. f. physik. Chemie B19, 203 (1932). 
b. Zener, Proc. Roy. Soc. A137, 696 (1932); A140, 660 
(1933). c. Landau, Physik. Zeits. Sowjetunion 1, 88 (1932); 
2, 46 (1932). d. London, Zeits. f. Physik 74, 143 (1932). 
935" and Syrkin, Acta Physicochimica 2, 433 
1 ; 

3a. Volmer and Kummerow, Zeits. f. physik. Chemie B9, 
141 (1930). b. Nagasako and Volmer, ibid. B10, 414 (1930). 
c. See also Hunter, Proc. Roy. Soc. Al44, 386 (1934). 

4a. Wigner, Nachr. Ges. Wiss. Gottingen, p. 375 (1927). 
b. Pelzer and Wigner, Zeits. f. physik. Chemie B15, 445 
(1931-32). c. Eyring, Chem, Rev. 10, 103 (1932). d. G. 
Herzberg, Zeits. f. physik. Chemie B17, 68 (1932). 
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only about 0.001 as fast as the normally expected 
rate because of the change of multiplicity in the 
decomposition. In the ordinary thermal process 
a nonparamagnetic linear NNO molecule passes 
from a ‘2 state through an activated state 53 kg. 
cal. higher to form O(?P) and N2(!2). 

We next consider the properties of the mo- 
lecular potential energy. Among the conditions 
which we require for a satisfactory potential 
function for N2O in its normal state are: 


(1) Removal of an O atom to infinity must leave the 
function for a normal Ne molecule. 

(2) Removal of an N atom to infinity must reduce the 
potential to that for a 42 NO molecule. 

(3) N2O being a linear molecule with O on one end, the 
function must yield the same values for an O atom ap- 
proaching either N atom. 

(4) The minimum must correspond to the correct 
dissociation energy and moment of inertia (interatomic 
distances). 

(5) The potential must yield the known vibration fre- 
quencies. 


Certain less obvious conditions on the potential 
are next considered. When N:O decomposes 
adiabatically to give a 'Z Ne molecule the oxygen 
will of necessity be in a singlet state. Now O (‘D) 
lies 45.1 kg. cal. above O (*P),° whereas the next 
higher level, O (1S) lies 96.1 kg. cal. above the 
’P state. Because O (1S) is so much higher than 
O (‘D) we conclude that the state of oxygen in 
the normal N,O is such that if it dissociated 
without change of multiplicity into '2 Nez and 
therefore ‘singlet oxygen the latter would be 
in a 'D state. Now Ne ('Z) dissociates adia- 
batically into two nitrogen atoms each in a 
4S state. Thus one thinks of a !2 N.O as as- 
sembled without change of multiplicity from two 


5 Bacher and Goudsmit, Atomic Energy States (New 
York, McGraw-Hill, 1932). 
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Fic. 1. Potential energy surface for N2O for linear configurations. Contours, 10 kg. cal. intervals. 


‘S nitrogen atoms and one 'D oxygen atom. 
When a ‘4S nitrogen is removed from NO 
adiabatically it must leave a ‘2 NO molecule 
behind which is formed from a ‘4S nitrogen atom 
and a 'D oxygen atom. This corresponds to the 
NO (on‘x?, 4=-) of Mulliken.6 We use as his 
estimated values for the constants for the cor- 
responding Morse curve an equilibrium distance 
m=1.28A, a vibration frequency w=3000r,-* 
= 1430 cm™ and a heat of dissociation plus zero- 
point energy of 40 kg. cal. The zero-point energy 
of 2 kg. cal. leaves 38 kg. cal. for the heat of 
dissociation. 

We shall use a potential function V(r, re, y) 
for NO in which the coordinates and masses are 
as indicated in Fig. 1. 

The following definitions are introduced for 
convenience, where it is to be understood that 
the partial derivatives are to be evaluated at the 
point 7;=1.22; re=1.10; and y=180°. 


eV eV OnV eV 
Be ef anne as f ooo = = 
1 ( =") ’ ke ( =] ’ Ryo (2%) ’ ks (= ) ’ 


Hi = (mym_+-myms) (mm, +m2+ms)—, 


* Mulliken, Rev. Mod. Phys. 4, 51 (1932). (See Fig. 48.) 


ba = (mym2+ myms) (mi +m2+ms3), 

Bs =mm2(m,+m2+ms3), 

bss = P1772" (wie — bs’) (uit? + oto? + 2usrit2) 
and 

As = 276i, 
where c is the velocity of light and w; is the ith frequency 
of N2O in wave numbers. 


Now a straightforward solution of the me- 
chanical problem for three atoms (see Cross and 
Van Vleck)’ relates the frequencies to the force 
constants by the following equation. 


Mid? —ky 
Ms\? — Rie 


0 0 


usd? — Rye 0 
bod? — ko 0 
bss? — ks 


(1) 


The fact that N,O is linear is responsible for 
Eq. (1) reducing to a linear and a quadratic 
equation. The k’s are now determined to fit the 
observed frequency as given by Plyler and 
Barker.* For the two equal (degenerate) bending 
frequencies these authors give 589 cm-, so that 
k3=6.4X10-" ergs mol.-' rad.-*. They give 


7 Cross and Van Vleck, J. Chem. Phys. 1, 350 (1933). 
8 Plyler and Barker, Phys. Rev. 38, 1827 (1931). 
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w,;=1285.4 and we=2224.1 cm@ for the other 
two frequencies. These two values, then, fix the 
two roots, \”, of the quadratic. Since there are 
three k’s to determine however, and since in any 
case ky. must be comparatively small, we try 
taking it equal to zero. This leads to complex 
values for k; and ke which are impossible. Thus 
we cannot neglect the interaction between the 
two outer atoms. We eliminate one k by the 
requirement that k1/ke= 15.85 K 10°(22.2 K 10°)", 
which is the ratio of the force constants of a nor- 
mal*NO and a normal Ne molecule.® This gives 
ky=12.0X10°, ke =16.7X10° and kiz=0.9X10°. 
The units are in ergs mol.—! cm. 

For the equilibrium distances in the NO 
molecule we take the values proposed by Paul- 
ing,!° i.e., for the N-N distance r2=1.10A, and 
for the N-O distance 7;=1.22A. These values 
give a moment of inertia of 66.7X10-* as 
compared with 66.0 X10-*° found by Barker." 

We. now give a potential function for NO 
which satisfies all the conditions enumerated 
except the third. It is 


V=D,f(qi, ri, 61) +Dof(c2, re, de) 


+ED flas 11, 0) f(s, 72, di) 


+Dsf(as, 11, bs)f (Cs, 2, ds) sin’ y. (2) 
Here 
f(a, 71, 6) =exp—(2a(r1—6)) —2 exp—(a(ri—4)) 
and 
f(c, r2, d) =exp—(2c(r2—d)) —2 exp—(c(r2—d)). 


If to (2) we add all the additional terms obtained 
by replacing 7; by 73; and y by 4, where 6 meas- 
ures the angle between 72 and 73, then the third 
condition is also satisfied. The inclusion of these 
terms considerably complicates our potential 
without changing it appreciably in the region 
in which we are interested so that they are 
omitted. 

The‘constants used in Eq. (2) are collected in 
Table I. 


9 Mulliken, Rev. Mod. Phys. 4, 1 (1932). 
10 Pauling, Proc. Nat. Acad. Sci. 18, 293, 498 (1932). 
! Barker, Phys. Rev. 41, 369 (1932). 
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TABLE I. Constants used in potential function (2). 











D a b ¢ d 
1 40. 3.95 1.28 — _ 
2 172.82 — _— 3.10 1.09 
3 — 39.97 3.00 1.12 3.00 1.20 
4 — 57.89 3.00 1.12 3.00 1.10 
5 25.97 3.95 1.22 3.00 1.20 
6 74.27 3.00 1.22 4.00 1.10 
7 — 60.00 2.519 122 1.665 1.10 
8 46.83 3.00 l.2e 3.00 1.10 








The adjustment of the constants for a function 
of the type (2) is a simple matter but should 
be made in a certain sequence. 

(1) The first two terms are simply Morse 
functions for Mulliken’s NO (oe7‘z?, 42-) and 
for the normal nitrogen molecule, respectively. 

(2) The third term takes care of ki2. The con- 
stants a3 and c; are chosen as 3, approximately 
the value of a2 in nitrogen. b3; was chosen to con- 
tract the N—O bond while d; has the effect of 
slightly expanding the N—N bond. After these 
choices were made, D; is taken to make kj 
= 0.9 10° ergs mol. cm~. 

(3) The fourth term is chosen to make 0V/dr; 
=0 at 7;=1.22A and r2=1.10A. Now d, is taken 
equal to 1.10 so that f(cs, re, ds) = —1. Finally 
a, and b, are kept with the values of a; and ); 
for the same reason as before and Dy, is so 
chosen as to make dV /dr,=0. 

(4) The fifth term is arranged in a manner 
analogous to the fourth to make 0V/dr2.=0 for 
r,=1.22A and r2=1.10A. 

(5) In the sixth and seventh terms ky, ke and 
E are fitted. Here b, and 6; are set equal to 1.22, 
the equilibrium N-—O distance in N.O, and d, 
and d; equal to 1.10 the corresponding N-N 
distance. Thus, for 7;=1.22 and r2=1.10, for 
these two terms, 


f(a, ri, 6) =f(c, re, d) = —1 
Of(a, 71, 6) /dr1=Of(c, r2, d)/dr2=0 


df(a, 11, b) /or? = 2a*D, 
while 
O°f(c, 72, d)/de? =2cD. 


In this way the contribution to the energy is 
D;+D;. These are then chosen to give N2O the 
correct dissociation energy and at the same time 
to give values for the a’s and c’s in the neighbor- 





= | ~ ~~ 


~~ 


- me es om ee 


no 





|.09 
|.20 
| .10 
|.20 
|.10 
1.10 
1.10 


ion 
uld 


rse 
und 
aly. 
on- 
ely 


the 
time 
bor- 





hood of 3. as and a; are now chosen so that 
#V/dr;?=k,: and c, and c; so that 0?V/dre? = ko. 

(6) The final term, of course, adds nothing to 
energy, slope or curvature in the 7, and f2 
directions so long as the molecule is linear, i.e., 
y=0 or x. Ds is obtained by solving the equation 
b33\°=k3, and putting Ds=}k3. There is some 
arbitrariness particularly in our choice of as and 
“;. This could be removed by fitting the potential 
function to give more exactly the higher fre- 
quencies found by Plyler and Barker in the 
manner employed by Adel and Dennison for 
COz." A simple central force potential of the 
type they have used for CO, will not, in general, 
be valid for distances far from the equilibrium 
position, as they also point out. Our potential 
can hardly be seriously in error providing 
adiabatic decomposition gives the products we 
have supposed. 

The method of using polynomials in the f’s 
for the various bonds provides a way for con- 
structing potential functions to fit all the limit- 
ing conditions for even the most complicated 
polyatomic molecules. This should be a generally 
useful procedure. 

Fig. 1 gives the potential surface for linear N2O 
plotted from Eq. (2). Here 7; is plotted against 
ro/f where f=[m(m3+mye)/m2(m,+ms3) }', and 
the angle between the abscissae and ordinates is 
given by 


7/2—sin— [mom,/(m3+me)(mi+ms) }} 
= 58° 54’. 


The potential function (2) gives for the two 
nonbending normal coordinates 


—”; =. 0.4679(r; - Yo.) +0.8828(re a Yo) 
and Q2=0.7842(1r1 — 11) —0.6206(r2— 12). 


Here 7%; and 72 are the equilibrium distances. 
These give the directions of the major and minor 
axes of an equipotential line in the immediate 
neighborhood of the minimum if one plots it in 
the way used in Fig. 1. Because the potential in- 
volves terms higher than the second power in 
r, and 72, by the time one has reached an equi- 





1933) and Dennison, Phys. Rev. 44, 99 (1933); 43, 716 

3). 

assy and Polanyi, Zeits. f. physik. Chemie B12, 279 
t. 
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potential corresponding to the zero-point energy 
(6.66 kg. cal.) the major and minor axes of the 
near ellipse will have changed. However, even in 
the present unsymmetrical case, this change in 
the coefficients of the r’s is less than one percent. 
In the same way we expect our vibration fre- 
quencies to show a parallel small anharmonicity. 

We next consider how the potential surface 
for O @P) approaching Ne ('2) will cut the sur- 
face just drawn. At large distances the O attracts 
the near N in Ne due to the van der Waals 
forces and the Coulombic integrals, but it repels 
it by an amount equal to 3 all the exchange 
integrals between the two atoms.'* Now we take 
the attractive potential as 0.14 of the Morse 
curve value for NO in its normal state (0?r‘z, ?II 
in Mulliken’s notation) and the _ repulsive 
potential as 0.43 of this value, giving a net 
repulsion of 0.29 of the Morse curve value.'* 
For this Morse curve we take w=1906 cm, 
¥9=1.15A, a=3.05, heat of dissociation = 121.94 
kg. cal. and the zero-point energy = 2.71 kg. cal. 
In Fig. 2 the curve I is a section following the 
lowest part of the east west valley of Fig. 1, 
while curve II is —0.29 times the Morse curve 
values for normal NO. 


ABSOLUTE RATE OF REACTION 


Curves I and II intersect when 7;, the N-O 
distance, is 1.73A. This point of intersection of 
I and II is the activated state for the decomposi- 
tion of N,O. In the activated state the vibration 
in the r; direction becomes a translation, and our 
potential function gives us for the N to N 
vibration the value 2005 cm, while the de- 
generate bending vibration has become 361 cm. 
This gives us 3.9 kg. cal. for the zero-point 
energy in the activated state and an activation 
energy at the absolute zero of 50.0 kg. cal. 

Now not every system passing through the 
activated state will make the transition from 
I to II. If we indicate the chance of transition, 
averaged over the velocity distribution charac- 
teristic of the temperature, by x, we have for 
the specific reaction rate constant! 


14 Eyring, J. Am. Chem. Soc. 54, 3191 (1932). 
4a This method of estimating the potential energy be- 
tween molecules with these same percentages gives the cor- 
rect activation energy for the conversion of para- to ortho- 
hydrogen and very close to the right activation energy in a 
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3 4 
k’=xII (1—exp —hy,*/kT)—I*7—' II (1-—exp —hv;/kT) (exp —Eo/kT)RT/h (3) 


i=1 7j=1 


so that for the ordinary activation energy at temperature 7 we have 


d\n k’ 3 hy;* hv*\—! 
E=RT? -(= einen (1—exp — ) 
dT kT 











i=1 kR 


This gives for E—E, at 850°, 900°, 938° and 
1000°K the values 1.751, 1.799, 1.833, and 1.881 
kg. cal., respectively. Thus our calculated value 
for E is approximately 52 kg. cal., to be com- 
pared with the measured value of Nagasako and 
Volmer*® of 53 kg. cal. Hunter** agrees ap- 
proximately with Nagasako and Volmer for low 
pressures, but he gives values ranging from 50.5 
to 65 kg. cal. increasing with the pressure. He 
postulates three different mechanisms. We shall 
discuss his results later. 

We are now in a position to calculate k’ 
theoretically except for the factor «x in (3). This 
will also be estimated theoretically, but before 
doing this we determine the value to be expected 
from experiment. For this purpose we take 
Ey = 53,000 — 1.833 cal. The moment of inertia of 
the activated complex, I*, is 10.14X10-**, while 
I is 6.67X10-**. The appropriate frequencies 
have all been given. Now taking k’ = 0.00192, the 
value given by Nagasako and Volmer for high 
pressures at 938°K, we find «=1.882X10~. If, 
instead of using this EZ) derived from Nagasako 
and Volmer, we use our value of 50,000 cal. we 
find «x=3.5X10—. 

Now Zener?» (see also Landau”) has given the 
relationship 


Poe (x0 


P» is readily found by numerical integration, 
and if one takes c=47°e19?m!/(2kT)!h|si1—se| the 





wide variety of other cases. The very good agreement 
between the calculated and experimental values as given 
below was for this reason not unexpected. 


4r? €12°m 
_—~}) exp (—p'/2miTyap( f 
hp|si—Se| 0 


4 hy; 


hy;\— 
=F —e-y no(1—exp -—) +1) R+E. (4) 


j=1 





Py= 1 —exp|{ —4ne19"/hv| si —Se| } (5) 
for the chance that at the point of crossing of two 
potential curves, such as I and II in Fig. 2, 
the system will jump from one curve to the other. 
Here v is the velocity of the system, |s:—S2| is 
the absolute value of the difference of the slopes 
of I and II at the point of crossing. The inter- 
action energy is ¢12= { gif’ gedr where ¢; and ¢» 
are the eigenfunctions for the two states which 
cross, and H’ is an additional potential term 
(generally due to magnetic interactions) omitted 
from the Schrédinger equation in the approxi- 
mation which gave the y’s as a solution. We 
find, as a result of this new term in the energy, 
that instead of crossing the two lower branches 
of I and II unite to form a single new state and 
the two upper branches unite to form a second 
new state. These two new states, instead of 
touching, are separated by a distance of 2¢12 at 
the nearest point of approach, the activated 
state. The chief condition to be fulfilled for (5) 
to apply is that ¢;.<3mv*, where m is the re- 
duced mass of the system. In applying P12 to a 
chemical reaction we must average over-all 
velocities. Thus the quantity we require is 


ie) 


exp (—p'/2meT)dp) ; (6) 





values given in Table II are found for Pj: as a 
function of c. 

Systems in statistical equilibrium moving up 
along curve I in Fig. 2 have a chance Py, of 
passing to II per crossing in the forward direction. 
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If they fail they have a second chance, Py2, of 
passing to II upon their return. Thus, so long as 
Py. is small, the total chance of reaction per 
system passing through the activated state is 
2P:2 (more exactly 2Pi2—P,,?). We therefore 
set 2Pi2=x. Now |si—se| from Fig. 2 is 8.50 
X10-* ergs mol.-! cm. Using the value of « 
found when £, was taken equal to 53,000 — 1833 
cal. we find, using Table II, that €;2.=5.05 cal. 








TABLE II. 
é 1 0.1 0.01 0.001 0.0001 0.00001 
Pw 0.212 .276 0531 00774 .001027 .000128 


c 0.08 0.07 0.056 0.054 0.0548 0.0538 0.052 0.0°1 
Piz 07107 =.09942)  .01816 §.0'728 .01662 .05541 .05299 .05156 








per mole. Thus the amount, 2¢;.=10.1 small 
calories, by which the two upper branches of 
curves I and II fail to make contact with the 
corresponding lower branches is not visible on 
the energy scale we have chosen. Had we used 
our theoretical value E)=50,000 cal., we would 
have found for €12 7.15 cal. Either of these values 
for €y2 fulfill the condition of being small com- 
pared with the average kinetic energy RT/2 
=938 cal. In our averaging of Py. we have in- 
cluded systems having kinetic energies equal to 
and lower than €;2. This introduces some error 
into our calculation. However, the fraction of all 
particles having at the activated state a kinetic 
energy of €;2 and below is small, and even in this 
region Py» is approximately correct. Thus treat- 
ing this region more carefully could hardly 
change our results by more than a few percent. 
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THE NATURE OF €19 


Besides the kinetic energy of the electrons 
and the electrostatic potential energy of electrons 
and nuclei which are ordinarily considered in 
obtaining the energy levels which cross and which 
do not contribute to €12, we have the following 
types of interactions which do contribute: 


(a) The relativity corrections, 

(b) A perturbation associated with the kinetic energy 
of the nuclei, 

(c) The interaction of the spin of an electron with its 
orbital magnetic moment, 

(d) Spin-orbit interaction between two different elec- 
trons, 

(e) Spin-spin interactions of the electrons inside the 
oxygen atom, 

(f) The magnetic interaction of the spins on the nitrogen 
molecule with those on the oxygen, 

(g) Terms responsible for the hyperfine structure. 


Now all the effects except (f) involve only a 
single atom so their magnitudes can be estimated 
from atomic spectra. Of these interactions only 
(c) and (e) seem as large as the experimental €;». 
Condon” finds from experimental data that 
the spin-orbital interaction between O ('D2) and 
(®P2) is 106 cm~, or 302 small calories per mole 
for the isolated atom. If this value for €;2 still 
held at the activated state, Table II shows us that 
x would be around 0.1. Actually, however, what 
must happen is that from both the 'D and *P 
states of oxygen molecular orbitals are formed 
with no angular momentum along the molecular 
axis, i.e., with \=0, while the other components 
of the angular momentum precess so rapidly 
about the molecular axis that coupling with the 
spin is partly destroyed. We still have spin-spin 
interaction inside the oxygen atom and between 
the nitrogen molecule and the oxygen. Of these 
two effects we should expect the latter to be 
negligible compared with the former simply on 
the basis of the distances between the electrons. 
The rough estimation of this effect of the nitrogen 
molecule from the point of view used by Wigner" 
in discussing ortho-para conversion indicates that 
we can neglect it. 

The operator for spin-spin interaction is 


D> (4u*((S:S,)ri? —3(Siti;)(Sirs;))}/ri®, (7) 


i<j 


46 Condon, Astrophys. J. 79, 217 (1934). 


6 Wigner, Zeits. f. physik. Chemie B23, 28 (1933). 
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where 7 and j are summed over the four # elec- 
trons, giving six terms in the sum. §; and §; are 
the vectors corresponding to the angular mo- 
menta of the spins of electrons 7 and j, respec- 
tively, and r;; is the vectorial distance between 
the electrons while 1;; is this distance. uw is one 
Bohr magneton, eh/4amc. Now the order of 
magnitude of (7) is 6y?r1.-*. If the average of rie 
is taken to be 1 Angstrom, we obtain 7.3 cal. per 
mole, the order of magnitude given by our calcu- 
lations for €,2. The theoretical treatment by 
Breit of spin-spin interaction in He and Li', 
taken with the experimental work of Houston 
and Schiiler, has been presented in convenient 
form by Bethe!’ and indicates that the above 
estimate for oxygen is right as regards order of 
magnitude except that there is still a factor de- 
pending on the symmetry of the eigenfunctions. 
This factor for a singlet-triplet interaction is 
zero. Thus only spin-orbit interaction contributes 
to €12. A part of the smallness of x is of course 
due to tunneling, which requires a corresponding 
revision upward of é1e. 





i*g* (2am*kT)! 827(82°A*B*C*)(RT)! 2» 





k! = k- 
1122129 h’ 


-@— Folk? — 


STEARN AND H. 


EYRING 


k FOR OTHER REACTIONS INVOLVING A CHANGE 
OF MULTIPLICITY FOR OXYGEN 


If we can calculate the c defined in connection 
with Table II, we can obtain x = 2P;2 immediately 
for any reaction. The four quantities required are 
the reduced mass of the oxygen against the 
rest of the molecule being formed or decomposed, 
the difference of the two slopes at the point of 
crossing, the temperature, and €;2. Now €12 will 
be nearly the same as in N;O so that only 
|s3—Se| offers any difficulty. 

As a further example we consider the reaction 
H.+O=H.0 at room temperature. Bear and 
Eyring'® in their Fig. 3 give, for the value of 
|s3—S2|, 8X10-* ergs mol.-' cm. m=2.934 
X10-*4 g mol.-!, €1.=5.05 cal. mole. We take 
T =300. This gives c=1.5X10~. 

Thus in general for a reaction involving a 
change in multiplicity, such as CO+O =COsz, we 
would expect « to be equal to 10-* within a 
factor of about 10. 

The specific reaction rate constant! for the 
reaction 


H:+O=H,0 (8) 


’ 


II (1—exp (—hv,*/kT))-! 
2h? int 


(1—exp (—hve/kT)) 





h hs 
By neglecting the frequencies which are high 
this gives 


k= (x/gi)(m*/m me) 1h(A*B*C*)*I.1¢e—20/kT 


= «1.83 KX 10~-"e—#0/#T cm? mol.— sec.“ (10) 


(11) 


The subscript 1 refers to 0 and 2 to He and the 
activated complex is indicated by an asterisk. 
m and y are the masses and frequencies, re- 
spectively, while 7 and g refer, respectively, to the 
statistical weights due to the nuclear spin and to 
the electronic states. g; is readily found to be 
6.76 at 300°K; the product of the three moments 
of inertia, (A*B*C*), is 1.9110-"* according 


17 Bethe, Handbuch der Physik, Vol. 24, No. 1 (Berlin, 
Julius Springer, 1933). 


= «1.11 X10'%e—#0/*T cm’ mole sec. 


—_— (2ame2kT )} 


Sr2lekT\—! 
———) i (9) 


h’ 2h? 





to the dimensions given by Bear and Eyring," 
while J2=4.56X10-. Taking E,)=8 kg. cal. ac- 
cording to these last authors, and x=1.5X10™, 
we find that k’, from Eq. (11), is 269, a very 
slow atomic reaction. Thus oxygen atoms in the 
presence of hydrogen undoubtedly disappear by 
some mechanism other than (8), as for example 
H.+O=OH+H, in which the multiplicity 
change can be avoided. 


PRESSURE EFFECTS IN THE N2O REACTION 


It is interesting to consider briefly Hunter's 
results* on the effect of pressure on the rate of 
reaction. Since we may write for the specific 
reaction rate!® 


18 Bear and Eyring, J. Am. Chem. Soc. 56, 2020 (1934). 
19 Wynne-Jonesand Eyring, J. Chem. Phys. 3, 492 (1935). 
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kT 
k'= K*—-=exp — 
h RT 


(12) 





[AE*+A(pV)*] ae 
exp — 


h 


R 


we may, of course, apply the usual thermodynamic considerations. Evans and Polanyi?® have dis- 
cussed the effect of pressure on reaction rates from a similar point of view. The process of taking 


partial derivatives gives 


(—=") (=) (—= (~*) (~) 
wee OF; ap ), av/ \aT/, 


Because of a misprint in their article they have a 
positive sign before the last term in (13). This 
does not affect their final results since, of course 
(dp/dV)r is a negative quantity so that they 
find 


((d log k)/dT)y=1.12((0 log k)/AT), (14) 


for the reaction 


C;sH;N+C.H;I =CHWN’ 
\ 
CH; 
in acetone. Now Hunter measured ((0 log k)/ 
dT)», so that we cannot invoke (12) to explain 
his high activation energies as we might had he 
measured ((0 log k)/dT)y. 
There is still another pressure effect also 
considered by Evans and Polanyi. From (12) we 
obtain at once 








d log k’ AH*+RT AE*+RT pAv* 
So 
Pp 


dT RT? RT? = RT? 


We assume that the NO molecule excludes other 
molecules from a prolate spheroid formed by the 


935) and Polanyi, Trans.’ Faraday Soc. 31, 875 


(13) 





rotation of an ellipse about the molecule’s 
major axis a=b+r;3 where b is the minor axis 
and 73, as before, is the distance between the O 
and the far nitrogen atom. This volume is then 
(x/6)(b+173)b?. Now the molal volume of liquid 
NO is 36 cc, and 73;=2.32A in the normal 
molecule so that we have }=4.16A. Thus by 
increasing 73; to 2.83A, as in the activated state, 
the molal volume is increased by only 2.8 cc. 
Consequently at a pressure of 40 atmospheres 
the quantity pAV increases the activation energy 
by only 2.7 small calories. We must therefore 
conclude that if Hunter’s interpretation of his 
experimental results is correct, then at high 
pressures the decomposition goes by some new 
mechanism, as he suggests. 

Although the picture of the nonadiabatic 
unimolecular decomposition of N2O at moderate 
pressures with an activation energy of approxi- 
mately 53,000 calories seems very satisfactory, 
an explanation of Hunter’s high pressure results 
is not given. 

We wish to take this opportunity of expressing 
our appreciation to Dr. Wendell Jackson for 
joint consideration of some of the initial phases 
of this problem, and also to Dr. G. E. Kimball 
for helpful discussions. 
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The absolute rate of the recombination of three hydrogen 
atoms is calculated entirely theoretically. The manner in 
which rotation determines the dimensions of the activated 
complex in cases having little or no activation energy is 
discussed. The theoretical data are in good agreement 
with the experimental rates of Steiner and of Amdur. An 
immediate consequence of the theory is that energy trans- 


fer occurs most effectively among particles which can 
react with each other, free atoms being more efficient than 
molecules. A qualitative application of potential surfaces 
to the problem of energy transfer as met in velocity of 
sound experiments and in experiments on maintenance of 
high pressure rates of unimolecular reactions is made. 





INTRODUCTION 


N previous papers? there has been developed 

and applied a theory of the rates of homo- 
geneous reactions which takes explicitly into 
account the behavior of all the degrees of freedom 
of the reacting system. In order to evaluate 
exactly the expressions that are obtained it is 
necessary to construct a potential energy surface 
for the system. In the applications that have 
already been made? ©) this was not feasible 
and it was necessary to obtain approximate 
information about the surface from the data of 
spectroscopy. For diatomic and triatomic sys- 
tems, however, it is possible to construct a 
surface which is in agreement with experimental 
data.* In this paper we shall show that from such 
surfaces one can calculate rates of reaction 
agreeing with experiment. From the surfaces 
considerable information is obtained concerning 
the mechanism by which energy is transferred 
during the three-body collisions by which atomic 
recombinations take place. 

We shall consider in detail the reaction 
H+H+X-—H.2+X. Before proceeding to this 
case, we first consider atomic reactions of the 
type A+B—AB. Although this reaction is very 
improbable, it is of great interest because from 
an investigation of the properties of the potential 
surface one learns how to take into account the 
relative angular momentum of two bodies which 
are going to react and also what is meant by the 


! Parker Traveling Fellow of Harvard University, 1934-5. 

2? (a) Eyring, J. Chem. Phys. 3, 107 (1935). (b) Gershino- 
witz and Eyring, J. Am. Chem. Soc. 57, 985 (1935). (c) 
Wynne-Jones and Eyring, }. Chem. Phys. 3, 492 (1935). 

$ 7) ie and Polanyi, Zeits. f. physik. Chemie B12, 279 
(1931). 
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collision diameter for such a reaction. This 
information can be carried over to associations 
of diatomic or polyatomic molecules in which 
there is no activation energy. 

These results are of interest also in connection 
with the important work of Pelzer and Wigner‘ 
who took less explicit account of the fact 
that molecules usually approach with relative 
angular momentum. 


THE REACTION AB—A+B 


The general equation for the rate of any 
reaction is 


, 


k= «(F.'/F,)(RT/h)e-®" 27. (1) 


This equation gives the concentration of the 
systems at a saddle point on the surface times 
the rate at which they are crossing the barrier 
“times the probability («) that once they cross 
the barrier they do not recross it. The usual 
potential energy surface for two atoms (Morse 
curve) is two-dimensional and it has no saddle 
point (maximum) except at infinity. It would 
seem as if Eq. (1), therefore, could not be 
applied. We must, however, remember that the 
Morse curve gives the energy of a system of 
two atoms as a function of the distance between 
them only for the case in which the atoms are 
approaching along the same straight line. In 
general this situation will not be encountered ; 


4(a) Pelzer and Wigner Zeits. f. physik. Chemie B15, 
445 (1932). (b) Pelzer, Zeits. f. Electrochemie 39, 608 
(1933). (c) Wigner, Zeits. f. physik. Chemie B19, 203 
(1932). 

5 Reference 2(a) Eq. (10). In accordance with the nota- 
tion used in reference 2(c), x has been substituted for c, the 
transmission coefficient, in order to avoid confusion with 
concentration. 
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Fic. 1. The distance between the atoms in the activated 
state in the bimolecular association of Hs and Ds. 


the atoms will have a relative angular mo- 
mentum, which may be considered as a rotation 
in two degrees of freedom. The kinetic energy of 
rotation for each quantum state and for each 
interatomic distance must be added to the 
potential energy arising from electronic inter- 
action, the latter being given by the Morse 
curve. The sum gives the energy that must be 
used in the exponential of the partition function, 
if we take as our zero of energy, for both the 
activated and normal states, the energy of the 
nonrotating molecule in its lowest vibrational 
state. For the case of two atoms we then have 
the activation energy Ey equal to D, the energy 
of dissociation. Now the addition of the energy 





of rotation to the rest of the energy of the 
molecule will cause a shift in the position of the 
activated state. For the dissociation of a dia- 
tomic molecule with no rotation there is no 
maximum, so that the activated state is at 
infinity. When we have a rotating molecule we 
do get a maximum and as the energy of rotation 
increases the distance 7; at which the maximum 
occurs becomes smaller. (See Fig. 1. These curves 
have meaning only for the integral values of J, 
the rotational quantum number.) 

This maximum locates the position of the 
activated state for the corresponding value of 
the rotational quantum number /. Knowing this 
we may calculate the activation energy for each 
value of J and by summing over all J’s we obtain 
the rotational partition function. Because of the 
Boltzmann factor, rotational states of very high 
energy are weighted out, so that at each temper- 
ature, for the association of two atoms there is 
an apparent collision diameter not very different 
from that for the collision of two diatomic 
molecules of the same kind. By differentiating 
with respect to 7 the total energy E= Eyin t+ Erot 
for each value of J and setting this derivative 
equal to zero, we find 7;, the separation between 
the atoms in the activated state for the corre- 
sponding rotational quantum number. Fig. 1 
gives the plot of 7; against J for He and Dz at 
300°K. . 

The activated state for the dissociation of a 
diatomic molecule thus lies at the top of a 
different maximum for each value of J. For the 
specific rate k, of the reaction 


AB—A---B—A+B 


we then have: 


Zatak(24(my+me)kT)*h 0,7 > (2J+1) exp {—(J(J+1)h?(822ur s*)-! 
J=0 


k= 


—(D—Ev))(RT)“}kTh- exp { —D(kT)-} 





(2) 


Snin(20(my-+ma)kT)h-4(1 —exp { —hv(kT)“} on) (2 +1) exp { —I(T+1)2(82ure2kT)-4} 
J=0 


=LM-' where L and M are the numerator and 
denominator, respectively, of the preceding frac- 
tion. For the reverse reaction 


ko=L{(2arm kT )*h-3(2amekT )th gy nim}. (3) 








xk is the transmission coefficient, i.e., the 
fraction of systems which, having crossed the 
maximum, do not return, m, and me are the 
masses of the atoms in the molecule; u.=mm2/ 
(m,+m2); ¢, and o, are the symmetry numbers 
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of the activated and normal states, respectively. 
v is the vibration frequency of the molecule. 
Za, 4, and i are the statistical weights due to 
nuclear spins and remain equal since ortho and 
para forms are always present in equilibrium 
amounts; however, for the statistical weights of 
electronic states we have g,=g,=1 and g,=4. 

The zero of energy is chosen to correspond to 
the nonrotating molecule in its lowest vibrational 
state ; hence D is the energy of dissociation from 
this level to the asymptote of the Morse curve. 
Ey is the energy of vibration measured from the 
lowest vibrational state as taken from the Morse 
curve at the activated point r; for each J. The 
rotational term in the denominator of (2) 
involves the constant 7, the equilibrium separa- 
tion of the atoms in the normal molecule, so that 
this moment of inertia is not a function of J. 
(We neglect the small interaction of vibration 
and rotation.) Thus at reasonably high temper- 
atures we can replace this summation by 
82°uro’kT(c,h?). The rotational partition func- 
tions for the normal and activated states of He 
and D2 at 500°K are given with considerable 
accuracy by the areas under the appropriate 
curves of Fig. 2. (The true curves are steplike.) 
The maxima of these curves give the most 
probable J’s for this temperature. With these 
most probable J values we can calculate the most 
probable values of 7; (the collision diameter). 
A more useful average 7; can be obtained from 
the equation 


8x2urPkTh-2= >> (2I+1) 


exp — {(J(J+1)h?(82?ur;??)-! 
—(V—Ey))(kRT)*}. (4) 


We find that for hydrogen 7;=4.4A and 7;=4.8A; 
for deuterium 7;=4.1A and 7;=4.9A. For Cle 
and Bre one gets 7;’s of the same order of magni- 
tude. In general, if the association of two 
particles, atoms or radicals, requires no activa- 
tion energy the apparent collision diameter, 
corresponding to the dimensions of the activated 
state, will be determined by the relative angular 
momentum, and at room temperature will be 
from 4 to 5A. We have applied this generalization 
in our treatment of trimolecular reactions.» 
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Fic. 2. Weighted number of states for each value of J. In 
normal state D=E,. (Bimolecular association.) 


From Eggs. (2a) and (2b) we can now calculate 
the rates of bimolecular association and uni- 
molecular decomposition of Hz and De. If we 
assume that « is the same for the reaction of D» 
and Hg for He we find 


ky=«3.6 X10!4e—?/2T sec.—!, 
ko = «2.78 X 10'* cm? moles sec.—!, 


and for De 


ky = «3.61 X10!4e-?/2T sec.—!, 
ko = «1.99 KX 10"4 cm? moles™ sec.~!. 


The transmission coefficient x is something like 
10-14 instead of being of the order unity as it 
frequently is when no non-adiabatic process is 
involved. ® 

We see that Hz and D2 decompose unimolecu- 
larly at the same rate except for the difference 
in zero-point energy. This is because the mass 
factors cancel out in the general equation for the 
rate and in addition the ratio F*,o¢/Frot is the 
same for Hz and De» although the individual 
rotational partition functions for the activated 
and normal states (F*,o¢ and Frot) for He are 
quite different from those for De. We also note 
that the rate ke for the unimolecular association 
of hydrogen atoms is nearly (2)! times the rate 
of association of deuterium. This is to some 
extent fortuitous, since it arises from the can- 
celation of a number of terms. This ratio of 2' 
for the relative rates strikingly exemplifies the 


6 For an excited diatomic oscillator with a reasonable 
sized dipole the interval before radiation is frequently 
around 10~? sec. In the case of hydrogen, where there is no 
dipole but only a quadrupole moment, we should expect a 
mean life of the order of seconds, so that x may well be as 
small as 1/y~10-"" 
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fact that one is not justified in postulating simple 
mechanisms or rate expressions on the basis of a 
single numerical agreement. 

The reason that the bimolecular association of 
the atoms plays a negligible role in the actual 
process of recombination is that x is so very 
small. The usual assumption will be shown to be 
correct that at atmospheric pressure about every 
thousandth collision among atoms is a three body 
one and forms a stable molecule whereas only 
about every 10~'* diatomic collision is stabilized 
by radiating energy. 


THe Reaction A+B+C—-AB+C 


We first consider the case in which A, B and 
C are all hydrogen atoms. A system of three 
atoms has nine degrees of freedom. In the normal 
state all these are free translations. In the 
activated state we have a triatomic quasimole- 
cule. For three hydrogen atoms for all except 
very large distances the configuration of lowest 
potential energy is that in which all three atoms 
lie on the same straight line. This important 
type we discuss first. The other important type 
to be treated later is the isosceles triangle 
configuration for the activated complex. The 
nine normal coordinates of the linear complex 
are three of translation of the system as a whole, 
two of rotation of the system, two for bending 
vibrations and two for motions along the straight 
line. We make the usual assumption that the 
two rotations and three translations are separable 
from the four internal coordinates. Of the four 
internal coordinates there are two transverse 
vibrations of equal frequency. Upon writing out 
in full the potential energy for the system the 
two transverse vibrations are seen to be separable 
from the two motions along a straight line up to 
terms involving the fourth power of the bending 
angle. Because this interaction is so small it is 
principally the transfer of energy between the 
two remaining degrees of freedom which is 
responsible for the stabilization (or dissociation) 
of the diatomic molecule by the third body. 
These last two degrees of freedom then are those 
which we now investigate. 

Eyring and Polanyi* in collaboration with 
Wigner have shown how to construct a potential 
energy surface for three collinear atoms such 
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that the classical motion of a mass point on this 
surface will describe the behavior of the system. 
Such a surface for 3H atoms is shown in Fig. 3. 
r,; represents the distance between atoms 1 and 
2 and 72 the distance between atoms 2 and 3. 
The contour lines are equipotentials spaced at 
intervals of 5 kg.cal. Both Figs. 3 and 4 were 
constructed in the usual way* with 14 percent 
Coulombic energy and the following constants for 
the Morse curve for hydrogen: heat of dissocia- 
tion = 102.4 kg.cal. ; vibration frequency wp = 4375 
cm~!; a=1.94A-!; r>5=0.74A. This surface has 
many features in common with other surfaces 
which represent the interactions of three atoms, 
any two of which can form a diatomic molecule. 
The angle that the axes make with each other 
is always 60° for three like atoms; otherwise it is 
a function of the relative masses.* For any set of 
three atoms the surface. will have the same 
general appearance. In particular, for the regions 
of high potential energy, in which we shall be 
principally interested, the surfaces are divided 
by some line into two regions, each of which is 
developable, i.e., the energy in each may be 
expressed as the sum of two potentials each 
depending on a single coordinate.’ It is this 
feature which enables us to discuss simply the 
transfer of energy between the two degrees of 
freedom associated with the two distances 7; and 
tr, and to determine the probability that a 
system with a given amount of energy will react. 
These two degrees of freedom of the linear 
molecule fix the relative translation of the three 
atoms. In what follows, for the sake of definite- 
ness, we shall refer to the hydrogen surface 
shown in Fig. 3 in which the dividing line is 
1:=1. The behavior of other systems of three 
atoms which can form diatomic molecules is 
entirely analogous. The cases in which one of 
the atoms is an inert gas or in which the cata- 
lyzing particle is a molecule will be discussed 
later. 

Although the surface is somewhat approximate 
in the region of low potential energy it is much 
more accurate for energies greater than 45 kg.cal. 


7 See for example Kimball and Eyring, J. Am. Chem. 
Soc. 54, 3381, etc. (1932) Figs. 1, 2, and 3 in which this 
feature is to be observed. In these particular cases, as well 
as for the three hydrogen atoms, the line makes equal 
angles with the coordinates 7; and re. 
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per mole. We shall see that the greater part of 
the systems that associate (or dissociate) cross 
the line 7:;=7r2 in the region of high potential 
energy, so that the part of the surface that is 
inexact does not enter into our calculations. The 
calculations in the section immediately following 
depend only on that portion of the surface 
corresponding to potential energies greater than 
45 kg.cal. 

The section of the surface to the right of the 
line 7; = 72 is developable in terms of the Cartesian 
coordinates x; and y,. The part to the left and 
above ”=72 is developable in terms of the 
Cartesian coordinates x. and ye, the latter of 
which is taken perpendicular to the axis 72. The 
angle between y, and yo and between x; and x2 
is then 120° (see Fig. 3). The energy of the mass 
point representing the system when it is in the 
first region is 


3u(a?+y:*) + V(x, 1), (5) 
dx (m,+m2)ms3 2 
where ¢=—etc.;5. 4«=——————__=-m 
dt M,+Mo+m3 3 


62 ero nas ¢ pron) 
10 


Angstroms— 


Fic. 3. Energy contour map of linear H; molecule. 


for our hydrogen atoms.® Similarly, in the upper 
part of the surface the energy is 


2U(k2?+ Yo”) + V(x2, yo). 
Along the line 7;=72 we have the relation 


This relation expresses the physical fact that 
suddenly as 7; becomes equal to 72 motion within 
the system which was pure translation becomes 
redistributed between vibration and translation. 
This redistribution also takes place in the vibra- 
tional energy. This abrupt change in develop- 
ability in one set of coordinates to developability 
in another probably often holds in ordinary 
molecules where there is no internal translation. 
Where this is found to be so it will provide a 
very convenient means of discussing many prob- 
lems. Let us-now suppose that in the x, y: valley 

8 See Eyring and Polanyi, reference 3. Fig. 16 and ac- 
companying text. Their coordinates —b and g correspond 
to our x; and 4». We should also like to mention that in 
their Eq. (25) on p. 307 sin g should be replaced by — sin ¢. 
Accordingly, the angle between 6 and c in Fig. 17 should be 
taken as 60°, as we have done in our Fig. 3, where 7; and r2 


correspond to b and c, respectively. We are indebted to 
Dr. O. K. Rice for kindly pointing out this error in sign. 
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there is a system (mass point) whose total energy 
is greater than 108.6+R(J) kg.cal., where R(J/) 
is the rotational energy corresponding to the 
angular momentum of the entire system. (This 
R(J) creates a ridge along the edge of the valleys 
converting the central plateau into a shallow 
basin, in the same manner as it forms the 
maximum of the potential energy curve for the 
diatomic molecules, which case we have discussed 
above.) If the point representing our system is 
to remain in the valley (in spite of having the 
energy 108.6+R(j), which is enough to dissociate 
it), this energy must be distributed in the x; 
and y; directions in such a way that the part of 
it in the y, direction, }my.?+ V(x, 1), is less 
than 108.6+R(J) kg.cal. Since the surface is 
developable in this region we can have no 
transfer of energy between the x: and y; direc- 
tions, so that the system cannot dissociate 
before it reaches the median line. When the mass 





b d f l 
=i f dxf ay dp. { dp, exp 
a c e k 


where /~ is the normalization factor. We deter- 
mine the limits a, b, c, etc., in the following 
manner. We are interested in the fraction of 
systems for which 


3u(ar?+y:?) + V(x, yi) >108.6+R (8) 
and at the same time 
3Uyo’ + V(x, v1) >108.6+R, (9a) 
which, by (6), becomes 
4u(0.5%1+0.86641)?+ V(x, y1) >108.6+R. (9b) 
There are the additional conditions 
0<ij,< © (10) 
and 
0< 3uy?+ V(x, y1)<108.6+R, (11) 


which restrict the systems considered to the 
valley. These four conditions taken together 
restrict us to the fraction of systems originally 
composed of a molecule and an atom in the x, 7; 
valley which dissociate into three atoms after 
crossing into the x2, ye valley. These conditions 
give the upper and lower limits of the probability 
integral H. The two strongest conditions are 
(9b) and the second inequality in (11) which we 
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point reaches this line it will pass into the other 
valley, since along the part of this line where 
V(n, 72) =45 kg.cal. (which is the region we are 
considering) no activation energy is needed, as 
the surface shows. As soon as it crosses this line 
the energy will be developable in terms of the 
coordinates x2, ye, 42 and we where # is given by 
Eq. (6). We see, therefore, that there is a definite 
probability that enough energy will be trans- 
ferred to the ye direction to make the system 
dissociate, i.e., cross over the ridge of the valley 
onto the plateau and remain there. By the 
principle of microscopic reversibility, this will 
also give us the probability that a system of 
three atoms will react to give a molecule and 
an atom. This is the reverse reaction in which 
we are, of course, also interested. 

There remains only the formal problem of 
setting up the required probability integral. We 
first consider 


—{(gu(tr+yr?+ V(x, y))(kT)“}, (7) 





will call (11b). (11b) gives the upper limit for 1.4% 
Wi = (2(108.6+R—V))ut=g. 


From (9b) we can obtain the lower limit for 7; 
for, solving, we get 


Ui = 2g —1.732x, = 2g — 3g. 


But this is the lower limit only until y= ~—g, 
otherwise condition (11b) would not be fulfilled. 
Letting 7: = —g and solving for #1, we get 4 = 33g. 
There is also a lower limit on %, since y cannot 
be greater than g, i.e., 4:=g3~*. We can, there- 
fore, divide our space into two regions, in one of 


which 


2g¢—3'm=Hx~g and g3t=%,=3'g; (I) 
and in the other 
—gs=yn=g and 3ig=ij=-. (If) 


If we replace the coordinates x and y by the 
coordinates measured along the line 1.=72 and 
n perpendicular to this direction, our potential 
at the activated state can be written as a function 
of & alone, to a good approximation. The rate 
with which the reaction proceeds is given by the 
density in phase space of the points representing 
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the reacting system per unit of length normal to 
the line 7;=72 multiplied by the velocity of these 
points in the 7 direction. Since we can neglect 
the small dependence of the potential upon 7 the 
integral over unit length /_,}dy simply yields 
unity. The limits of the integral over & as well 
as the potential along this line 7,=72 are de- 
pendent on the amount of rotation and therefore 
upon the quantum number J for this linear 
symmetrical activated complex. For J =0 corre- 
sponding to no rotation an interesting situation 
arises. To get the number of molecules formed 
we should integrate ¢ from 7;=72=0 to 1=fr=l 
where / is of the order of a mean free path, since 
at this distance the attractive potential will have 
a maximum due to external perturbations. 
Actually if all such collisions resulted finally in 
the formation of a stable molecule we should 
get a very large rate of combination. We must, 
however, consider the subsequent fate of the 
diatomic molecule formed with such a very large 
amplitude of vibration, i.e., the corresponding x. 
A subsequent collision normal to this axis will 
almost certainly decompose the molecule whereas 
a collision in the direction of the axis is apt to 
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stabilize it. Thus the chance that the molecule 
will persist after the next collision is approxi- 
mately d/l where d is a normal atomic collision 
diameter and / is the length of the molecular 
axis (the mean free path). We should, therefore, 
(when J =0 as for other cases) use approximately 
d rather than the mean free path / as the upper 
limit in the phase integral (7) (if x is to be taken 
as approximately unity) since this integral is 
very nearly proportional to the upper limit. 

Now instead of considering the surface for 
each value of J separately we use the surface for 
J=0 and integrate out to 7,=7r2=4.4A which is 
the position of the rotation ridge separating the 
plateau from the valley for the average J at this 
temperature (300°A). This procedure introduces 
but little error since the potential surfaces change 
only slightly with J and even choosing the 
position of the rotation ridge 1A too far out 
would only multiply the true rate by a factor of 
approximately 1.25. Underestimating the dis- 
tance to the ridge by a similar amount under- 
estimates the rate by a similar factor. 

Introducing the limits of integration into (7) 
gives the equality 


4.434 339 g , 
Hai [exp =| vi eryjaq ff exp — tae ayer) Paedy 
° . img te * =29g—233 


0 


+ f { exp — [4ulit-+ 2) (RT) tay | 
z=g3t =—g 









We have evaluated this integral numerically. The region fo which V(é) is less than 45 kg.cal. 
makes a negligible contribution. Even keeping the original assumption that the surface is ruled, 
the contribution along the line 7;=72 near the minimum is unimportant. Now the fact that the 
minimum is higher than corresponds to such an assumption decreases still further the contribution 
to H from this region since it introduces an activation energy for particles crossing 7;=72 near the 
minimum. We estimate the error made by neglecting the region where V(é)<45 kg.cal. to be not 
more than 10 percent. If we take for the zero the energy of three separated hydrogen atoms we 
find that H=2yukT-1.6410-°/h?. It can be shown that if all systems crossing the line 7;=72 with 
enough energy actually reacted, the corresponding partition function F1’ would be equal to 2xukT- 4.4 
X10-*/h?. The ratio H/H’ gives us a measure of the transmission coefficient x. We find that approxi- 
mately } of the systems that cross the barrier with enough energy react. In this case there is only 
one degree of freedom which can take away the excess energy. Where there are more, «x will be 
larger, and for polyatomic molecules it will approach unity. 
The rate of the association reaction k; is therefore given by 


2 
| 8x1*RT (oth?) IT (1 exp —{hvs(kT)—})-\(2aem*kT) h-411 0 
Lala i=1 
k3= ’ 
Znln (2 amk T) 9/24,-9 








(12) 
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where IJ*, o*, 6=(kT/2rm*)' and »; are, re- 
spectively, the moment of inertia, symmetry 
number, average velocity across the barrier and 
vibration frequencies of the activated complex; 
m is the mass of one of the atoms. gata(gnin) =} 
exactly as for the bimolecular association. In 
this expression the most uncertain quantity is 
the moment of inertia of the activated complex, 
since the atoms may be from two to four Ang- 
stroms apart when they cross the line 7;=72, 
and the distance / is also a function of J. We take 
as an average, as well as a most probable value, 
r,=3A. This makes our maximum error a factor 
of two, although the probable error is much less. 
The one other uncertain quantity is the magni- 
tude of the frequencies of the transverse vibra- 
tions. The absolute rate is not very sensitive to 
changes in this quantity, but without it we 
cannot predict the isotope effect with Hz and De. 
The evidence available seems to indicate that 
these frequencies should be somewhat lower than 
those given by surface IV. If we assume the 
frequencies are very low the vibration partition 
functions take the classical values kT /hy; 
=kT2rm'/hf? and we get a factor of m from 
the two frequencies of the activated state. This 
m in the numerator with the m? in the denomi- 
nator of (9) gives an m= as the final dependence 
of the rate on the mass. This agrees with the 
results of Farkas® and with those of Amdur.!° 
This agreement gives us some information about 
the values of the vibration frequencies that we 
must insert into Eq. (12). To give the classical 
mass dependence these frequencies must be of 
the order of magnitude of 100 cm~. Using this 
frequency we get 7 as the numerical value of the 
partition function for the two transverse vibra- 
tions. Although the bending frequencies are low 
the one associated with ¢ is sufficiently large to 
make the total zero-point energy approximately 
that of the normal state in agreement with the 
results of Farkas and of Cook and Bates.® 
Substituting in (12) the numerical values, we 
get k;=3.1 10" cm® moles sec.—!. 

The process just considered involves the three 
atoms 1, 2 and 3, where 2 lies in between. 
Atom 3 comes up to the unstable configuration 
1—2, atom 1 leaves and 2 and 3 form a stable 


® Private communication. 
10 Amdur, J. Am. Chem. Soc. 57, 856 (1935). 





molecule. This mechanism accounts for those 
cases in which the atoms are moving along a 
straight line, or at such angles as are included 
by the amplitudes of the transverse vibrations 
of our activated state. When the amplitudes of 
these vibrations are very large, i.e., for atoms 
approaching at very sharp angles, it is no longer 
legitimate to make the assumption that the 
energy is still separable in terms of the coordi- 
nates that were good approximations to normal 
coordinates for the linear case. This being so it 
may be there is some second region of our many 
dimensional energy surface which is as effective 
as that just considered. Such a new region would 
make possible an additional mechanism for the 
reaction. We find that this situation actually 
arises for configurations in which one of the three 
atoms is moving perpendicular to the line forming 
the other two. The pair moving along the line 
reacts and the third carries off (or for the 
reverse process furnishes) the excess energy. 


SECOND TYPE OF COLLISION FOR 3 H AToms 


Fig. 4 is the potential energy surface for two 
hydrogen atoms moving along a straight line, 
the distance between them being 7 with the 
third hydrogen atom moving perpendicular to 
and toward (or away from) the center of this 
line. The distance from the third atom to the 
center of gravity of the first two is 7. The three 
atoms thus always form an isosceles triangle. 
In Fig. 4 only a little more than half the surface 
is drawn as it is symmetrical about X =0. Here 
again we see that a line AB divides the surface 
into two regions in which different sets of 
coordinates are appropriate. In the region to the 
southeast of AB the most convenient coordinates 
are the parabolic cylindrical set," u, v and w, 
where 


x=c(v—u); y =2c(uv)}; 2=W. 


In this coordinate system our equipotential lines 
are nearly the parabolas v=const. In the region 
to the northwest plane polar coordinates p and @ 
with origin at the point c may be used and in 
this part of the surface the equipotentials are 
given to a good approximation by p=const. 


1! Adams, Smithsonian Mathematical Formulae (1922), 
p. 105. 
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Fic. 4. Energy contours at 5 kg.cal. intervals for 3 H 
atoms forming isosceles triangle. 


In these regions the surfaces are not develop- 
able, i.e., there is transfer of energy between the 
u and v directions and between the p and @ 
directions. However, because of the general 
shape of the surface in the neighborhood of the 
line AB we can treat it as consisting of two 
valleys intersecting at an angle of 130° by the 
same method as was used for the two valleys 
meeting at 60° in Fig. 3. From this point of 
view this type of collision is seen to be less 
efficient than the linear one first because the 
angle differs more from the most favorable 
angle of 90° and secondly because the height of 
65 kg.cal. for the minimum in the pass has the 
effect of excluding many more colliding systems 
than did the former 15 kg.cal. height. Therefore, 
this type of collision cannot as much as double 
our calculated rate. 

We have, then, that, within a factor of two or 
three, the calculated rate of recombination of 
hydrogen atoms is 3.1X10'5, as compared with 
the experimental results of Amdur!® and of 
Steiner” of 2X10! and 1X10"*, respectively. 


CASES WHEN THE THIRD PARTICLE IS NOT AN 
H Atom 


For definiteness we discuss the case when the 
third atom is helium and take the distance from 


12 Steiner, Trans. Faraday Soc. 31, 623 (1935). 
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He to the nearest H as 72. A new map will not 
be drawn. Instead the differences will be indi- 
cated by reference to Fig. 3. The surface instead 
of having an angle of 60° between the 7 and 12 
direction now has 50° 46’. The effective mass of 
the point representing the system instead of 
being half the mass is now 4/3 the mass of a 
hydrogen atom and distances in the 72 direction 
must be multiplied by (3)! before plotting.* Be- 
cause unexcited He cannot form a valence 
compound with the H atom the lower valley of 
Fig. 3 disappears. Instead we have only a very 
shallow van der Waals valley at a considerably 
larger value of re (around 3.75A). For smaller 
values of re the surface rises rapidly. Van der 
Waals attraction is of course present in the case 
of 3 H atoms but is unimportant compared with 
the valence valleys and so was not discussed. 
From the nature of this new surface it is clear 
that neither helium nor any other inert molecule 
can approach an H atom in efficiency of carrying 
away energy in a linear collision, since so much 
of the region which contributes to the reaction 
of 3 H’s is here shut off by a high hill, and since 
for the available region entrance of the mass 
point into the valley is necessarily at a less 
favorable angle. 

We now consider the case corresponding to 
Fig. 4 in which the inert gas approaches the 
center of gravity of the combining atoms normal 
to the molecular axis. The circular region about 
c now becomes a peak instead. If a particle is 
to leave the plateau and remain in the valley it 
must roll into the valley on the end formed by 
the peak. This is only possible if the relative 
velocities in the 7; and 12 directions lie within 
narrow limits and if the point is directed at the 
peak in just the right way. However, these 
requirements are not quite as stringent as for 
the linear case, i.e., the limits of integration in 
the integral H not quite as narrow. This is 
because the barrier across the valley formed by 
the peak leaves it more accessible from the 
plateau than in the linear case. Given the surface 
it is of course a straightforward problem in 
statistical mechanics to calculate the specific 
rate for this process. This qualitative discussion 
of surfaces shows in an unambiguous way how 
the increased accessibility of the end of the valley 
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to the plateau (which is associated with in- 
creasing reactivity) increases the transfer of 
energy between degrees of freedom. 


COMPARISON WITH EXPERIMENTAL RESULTS 


There is some disagreement as to the relative 
efficiencies of H and H, as third bodies in the 
recombination of hydrogen atoms. The work of 
Steiner indicates that He is at least ten times as 
efficient. The work of Amdur gives contrary 
results. We can examine these relative efficiencies 
from the theoretical point of view which we 
have adopted. In order to have a hydrogen mole- 
cule transfer energy by an exchange of partners, 
it must lie between two atoms on a straight line. 
Thus we have the configuration H—H-H—H 
and reaction will occur when the central 
hydrogen molecule splits and we get two new 
ones formed. There are three factors which 
might make this process less probable than the 
reaction with three hydrogen atoms. The first is 
that there is less probability of finding the four 
atoms on a line than of finding three. The 
difference, however, is not as great as one might 
expect, for the rotational partition function of 
hydrogen is so small that the molecules will have 
the correct orientation about half the time. 
Secondly, one-third of those collisions in which 
the atoms are collinear will have the molecule 
between the two atoms. These combined proba- 
bilities would make the hydrogen molecule 3 as 
effective as the atom. The third factor that could 
reduce the efficiency would be the necessity of 
an activation energy for the approach of two 
hydrogen atoms to a hydrogen molecule. We 
can make no definite statement about this, 
since we have not calculated the necessary 
surface. We thus come to the conclusion that 
Hz is at best 3 as efficient as H. This is in agree- 
ment with the results of Amdur rather than 
those of Steiner. Only if the transfer of energy 
between the rotation of the colliding Hz molecule 
and the associating atoms is much greater than 
our present knowledge of the potential surface 


indicates could the calculated order of efficiencies: 


be reversed. 
The reasoning of the preceding paragraph 


indicates why hydrogen molecules are so much 
more efficient than others for the transfer of 
energy. Other authors have assumed that it is 
because the mechanism is one of transfer of 
rotational energy to vibrational and that since 
hydrogen has such a small moment of inertia it 
can transfer a large amount of energy in a single 
quantum. From the preceding considerations it 
would seem that the efficiency of hydrogen may 
be due rather to the fact that it can so readily 
form a complex with the reacting substances, 
since its moment of inertia is so small that it 
is in a favorable configuration about half of the 
time at room temperature. 

We see that in general, according to the point 
of view which has here been put forth, the 
ability to transfer energy is intimately connected 
with reactivity. This is in agreement with most 
of the available data. Franck and Eucken" have 
reached the same conclusions through more 
qualitative methods. Our calculation for the 3 H 
case shows how any other case may be treated 
quantitatively. The method is valid not only for 
actual chemical reactions but also in cases in 
which no bonds are formed or broken. Thus it 
can be applied to the data obtained in sound 
dispersion experiments, where all the data indi- 
cate that energy is most easily transferred 
among molecules that can react with each other.'4 
The method may also be applied to the study 
of the maintenance of the high pressure rate of 
unimolecular reactions by foreign gases. Here 
too most of the observed facts are fitted by the 
theory. As we might expect, comparatively inert 
substances vary considerably in the efficiency 
with which they maintain the rate. This is 
shown by the work of Volmer and his collabo- 
rators on N2O.'5 Returning to the recombination 
of hydrogen atoms, Bonhoeffer'® has made a 
qualitative determination of the relative eff- 


ass and Eucken, Zeits. f. physik. Chemie B20, 460 
144 Eucken and Becher, Zeits. f. physik. Chemie B20, 467 
(1933); Kneser and Knudsen, Ann. d. Physik 21, 682 
(1935); Richards and Reid, J. Chem. Phys. 2, 206 (1934). 
18 Volmer and Froehlich, Zeits. f. physik. Chemie B19, 
89 (1932). 
16 Bonhoeffer, Zeits. f. physik. Chemie 119, 385, 475 
(1926). 
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ciency of various gases as third bodies. The 
order he finds agrees very well with the relative 
reactivities. The data of Senftleben and Hein,"’ 
however, are in definite disagreement with our 
theory, but the absolute rate for the recombin- 


17 Senftleben and Hein, Ann. d. Physik 22, 1 (1935). 
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ation of hydrogen atoms found by these authors 
is very much lower than that of Amdur or Steiner. 

The method of potential energy surfaces thus 
provides us with a clear picture of homogeneous 
catalytic processes. It forms a convenient basis 
for a satisfactory discussion of the problems of 
energy transfer. 
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A Summary of Experimental Activation Energies of Elementary Reactions Between 
Hydrogen and the Halogens 


J. CARRELL Morris AND Rosert N. PEAse, Department of Chemistry, Princeton University 
(Received July 22, 1935) 


Existing experimental data on reactions between hydrogen and the halogens have been 
analyzed to give values of the activation energies of the intermediate reactions involved. These 
are compared with calculated values due to Eyring and Wheeler. 


VIDENCE has been accumulating over a 
period of years that reactions between 
molecules and free atoms may have appreciable 
activation energies. This evidence is in the main 
not direct, but must be deduced from considera- 
tions as to the most probable mechanism of 
reaction. In this paper the experimental evidence 
for reactions of hydrogen and the halogens is 
reviewed and tentative values of the activation 
energies are suggested. The reactions in question 
are of the types: 


H+X,—-HX+X, X+H:—-HX+H, 
H+HX—X-+H2. 


HYDROGEN-BROMINE REACTIONS 


This reaction system has been quite thoroughly 
covered in both its theoretical and experimental 
aspects by Kassel in his Kinetics of Homogeneous 
Gas Reactions. However, newer numerical data 
have made possible more exact conclusions as to 
the individual reactions. The chain consists of 
the steps: 


ky 
Br2—2Br, (1) 

ke 
Br+H,.—HBr+H, (2) 


k 
H+Br.—HBr+Br, (3) 


k 
H+HBr—H.+Br, (4) 


> 


ks 
Br+Br-—Bry, (5) 


which lead to the theoretical expression : 


d{HBr ] LH: |[Br2}! 
oo = 2ko(ki/ks)? a 
dT 1+k./ks_LHBr ]/[Bre J 








From experimental researches by Bodenstein 
and Lind,! confirmed by Bach, Bonhoeffer and 
Moelwyn-Hughes? 2ke(ki/ks)! and ks/ks have 
been determined. Since k;/ki= Kp, the dissocia- 
tion constant of bromine, which is now well 
known, ke may be evaluated from the knowledge 
of 2ko(ki/k;)*. The temperature coefficient of ke 
will give As, the activation energy of reaction 
(2). If the data of Gordon and Barnes* on Kp 
are used, both sets of results yield A2(Br+He) 
=17.7 Kg.Cals. 


1 Bodenstein and Lind, Zeits. f. physik. Chemie 57, 168 
(1907). 

2 Bach, Bonhoeffer and Moelwyn-Hughes, 
physik. Chemie B27, 71 (1934). 

3 Gordon and Barnes, J. Chem. Phys. 1, 792 (1933). 
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From this value of A» theoretical reaction rates 
may be calculated and are found to agree quite 
well with observed values when reasonable 
molecular diameters (3 10~* cm) are used. 

Reaction (4), H+HBr—H:2+Br, is the reverse 
of the reaction just studied. The difference in 
activation energies between process (2) and- the 
one just given may then be taken as the heat of 
reaction. This may be calculated thermochemi- 
cally from the data of Gordon and Barnes,’ 
who give these values: 


Bro—2Br, 
HBr—3H2+ + Bro, 


AE, =45.23 Kg.Cals. 
AF,=11.99 Kg.Cals. 


and from the recent work of Beutler,‘ who finds 
for hydrogen 


H.—2H, AF, = 102.68 Kg.Cals. 


Proper combination of these leads to the result 


H+HBr—-H,+ Br, AE)= —16.7 Kg.Cals.** 


By direct subtraction a value, As(H+HBr) =1.0 
Kg.Cals., is obtained for reaction (4). This value 
of 1 kilogram calorie is in accord with the results 
of Cremer, Curry, and Polanyi,® who attempted 
to determine the rate of this reaction by passing 
atomic hydrogen into HBr. The reaction was 
too rapid for their measurements, which were 
made by means of the para to ortho-hydrogen 
conversion, so that they could only conclude 
that A, was less than 3.0 Kg.Cals. 

The ratio k,/ks3 was accurately determined by 
Bodenstein and Jung,® who found a value of 
1/8.4 practically independent of temperature 
from 0° to 300°C, indicating that the A/RT 
terms are the same for reactions (3) and (4), for 
otherwise the ratio k,/k3 would vary with temper- 
ature. The constant factor of 1/8.4 may be 
ascribed to an “orientation factor,’”’ but what- 
ever the explanation, it has no bearing on pres- 
ent reasoning. It follows for reaction (3) that 


A;(H+Bre) = 1.0 Kg.Cals. 


4 Beutler, Zeits. f. physik. Chemie B27, 287-302 (1934). 

48 Kassel gives —14.5 Kg.Cals. for the heat of this reac- 
tion. This seems to be an error. 

5 Cremer, Curry and Polanyi, Zeits. f. physik. Chemie 
B23, 445 (1933). 

6 Bodenstein and Jung, Zeits. f. physik. Chemie 121, 
127 (1926). 
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REACTIONS 


HYDROGEN-CHLORINE REACTIONS 


The attempt to evaluate some of the constants 
in the hydrogen-chlorine reaction is a more 
difficult task, since the reaction mechanism has 
not been as definitely settled as in the hydrogen- 
bromine reaction. The Nernst chain mechanism 
is well established, however.’ According to this 
scheme the probable reactions in the photo- 
chemical hydrogen-chlorine combination are 


ky 
Cle+hv—2Cl, (1) 
ke 
Cl+H,—-HCI+H, (2) 
ks 
H+Cl,-HCI+Cl, (3) 


k 
H+HCI—H,+Cl. (4) 


The chain-ending mechanisms may possibly be 
any of the following 


k: 
CI+Cl+M—>Cl,+M (5) 

ke 
CI+X—CIx, (6) 

ky 
H+X—HX, (7) 

kg 
H+H+M—H.+M (8) 

kg 

Rio 
H+0,—HO, (10) 

Ry 


In this scheme X represents any foreign chain- 
breaker other than oxygen, such as the wall, or, 
with Bodenstein,* a volatile compound of silicon. 
M is a third body, either He, Cle or HCI. 

The most significant experiment toward the 
evaluation of activation energies in this series of 


7For a complete discussion see Semenov, Chemical 
Kinetics and Chain Reactions, pp. 98-115. 

8 Bodenstein and Unger, Zeits. f. physik. Chemie B11, 
253 (1931). 
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reactions was a direct determination by Rode- 
bush and Klingelhoeffer® of the rate of reaction 
(2). By passing atomic chlorine into hydrogen 
they determined the probability of a chlorine 
atom’s reacting, using the formula 


HCl formed 


7 Cl introduced 





where Z’ is the number of collisions between Cl 
and H; calculated from kinetic theory, and W is 
the probability of reaction at each collision. 
Rodebush and Klingelhoeffer have calculated A» 
from the temperature coefficient of W. We have 
also calculated A» directly from the formula 


W =e-42/RT, 


The results by these two methods accord very 
well and give for A2,(Cl+H:2) a value of 6 
Kg.Cals. 

The problem may also be indirectly ap- 
proached through a consideration of the photo- 
chemical reaction in the absence of oxygen. 
Bodenstein and Unger,’ using the following set 
of reactions taken from the larger number 
presented before : 


k 
Cle+hy—2Cl, (1) 
ke 
Cl+H.—-HCI+H, (2) 
ks 
H+Cl.—HCI+Cl, (3) 
kg 
H+HCI—H,.+Cl, (4) 
Re 
Cl+X—CIx, (6) 
kz 
H+xX—HX, (7) 


developed this theoretical expression : 
adLHCl] Akoks: Iaps.L He }- [Cle ] 
dt  — kekr[Ho]+kskol Cle ]+kakol HCl] 





On the other hand they found experimentally, 
d(HCl)/dt=kIavs.[He ]. 


® Rodebush and Klingelhoeffer, J. Am. Chem. Soc. 55, 
130 (1933). 
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been 
checked by Potts and Rollefson.'® It could only 
be accounted for if reactions (4) and (7) were 
very slow compared to (3) and (6), respectively. 
When only reactions (1), (2), (3) and (6) are 
considered, the theoretical expression found is 


This empirical equation has_ recently 


d(HC1]/dt = (2k2/ks)[H21/CX Vans. 


which agrees with the theoretical equation if [X ] 
be constant. If we assume that all of the terms 
in this expression are practically independent of 
temperature except ke/ks, the temperature co- 
efficient will depend on ke/kg and gives the 
difference in activation energy between reactions 
(2) and (6). 

Hertel" found the temperature coefficient in 
the absence of oxygen to be 1.37 between 25 
and 35°. This calculates to a difference in 
activation energies of 5.9 Kg.Cals. Potts and 
Rollefson’® by the same method find the differ- 
ence to be 5.8 Kg.Cals. These values agree well 
with the absolute value found by Rodebush and 
Klingelhoeffer for reaction (2) alone, indicating 
that (CI+X) requires no activation. In any case 
this result is a minimum activation energy for 
reaction (2). We, therefore, write 


A2(Cl +He) = 6.0 Kg.Cals. 


The activation for the reverse reaction may 
be obtained from a coupling of the above result 
with thermal data. Best present day values as 
compiled by B. Lewis and von Elbe” give 


HCls$53H2+ 4Cle, 
Cl=22Cl, 


AE,= +21.98 Kg.Cals. 
AFE,=+56.9 Kg.Cals. 


A proper combination of these with the value 
given before for hydrogen dissociation yields 


Cl+H.—-HCI+H, AE,=1.0 Kg.Cals. 


By subtracting this from A2(Cl+He) one obtains 
the activation energy of the reverse reaction, 
A,(H+HCl) =5 Kg.Cals. 

The evaluation of the energy of activation for 
reaction (3) is less certain. The reaction occurs 


10 Potts and Rollefson, J. Am. Chem. Soc. 57, 1027 
(1935). 

' Hertel, Zeits. f. physik. Chemie B15, 325 (1931). 

12 B, Lewis and von Elbe, J. Am. Chem. Soc. 57, 612 
(1935). 
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rapidly at room temperatures": '4 but the meas- 
urements are not of such a type that activation 
energies may be calculated from them. The 
matter must be attacked indirectly through a 
consideration of experiments on the photo- 
chemical combination in the presence of oxygen, 
along with some other data. 

In the presence of oxygen reactions (1), (2) 
and (3) start and continue the chain. Reaction 
(4) may be omitted since it was shown that it is 
negligible compared to (3). The chain-ending 
mechanism is believed to be 


k 
H+0.+M—HO.4+M, (9) 


where M is a third body. The reasons for this 
belief will be brought out later. These reactions 
lead to the following kinetic equation: 


d(HCl] © 2kiks [Cle]? 
ko [M][Oz] 


‘icciintianasican iat [ P 
dt 





(a) 


Since with moderate oxygen content the chain 
length is still of the order of 100 and k,[ Cle } 
involves only the first step, the equation may 
be approximated to 


d(HCl] 2kks [Ch} 
dt ky. [MJ[O2] 





(b) 


The temperature coefficient of the rate in this 
case depends chiefly on k3/ko. Values of the 
temperature coefficient per 10°, found by various 
authors at room temperature are: 1.12 by 
Hertel," 1.19 by Padoa,'® 1.10 by Porter, Bard- 
well and Lind.'* The activation energy calculated 
from these values will represent (A3;—As). 
Evaluation gives for (A;—Ag) 2.1, 2.9 and 1.8 
Kg.Cals., respectively. 

Krauskopf and Rollefson!’ found (A3;—Ag) in 
a different way. The partition of the hydrogen 
reacting in a system containing oxygen and 
chlorine between H:,O and HCl! was directly 


18 Boehm and Bonhoeffer, Zeits. f. physik. Chemie 119, 
385 (1926). 

14 Marshall, J. Phys. Chem. 29, 842 (1925). 

4 Padoa and Buttironi, Atti. Linc. 25, 215 (1916); Gazz. 
47, 6 (1917). 

16 Porter, Bardwell and Lind, J. Am. Chem. Soc. 48, 
2603 (1926). 

17 Krauskopf and Rollefson, J. Am. Chem. Soc. 56, 
325 (1934). 
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determined. Applications of equation (a) above 
and the equation 
d(H,0 |/dt=ki[Cle |, 


assuming one molecule of water formed for each 
hydrogen atom reacting with oxygen leads to 





H. used in H.O Vy 
total He reacting 
2 


~ 342ks/ko-[Cle]/[Or]-M] 


From the observation that when [Oz ]/[Cl. ] is 
20, U=0.57 at 0° and 0.40 at 85°C, they calcu- 
late the difference (A3;—Ag) to be 1.6 Kg.Cals. 

All these values represent a minimum for the 
activation energy A3, since Ay is not likely to 
be negative. This gives A3(H+Cl:) is greater 
than 2 Kg.Cals. 

The upper limit may also be determined. 
Bodenstein and Unger® concluded from their 
data that HCI did not inhibit the reaction in the 
absence of oxygen, and that consequently re- 
action (4) was slow compared to reaction (3). 
This is confirmed by the work of Potts and 
Rollefson.!° Let us therefore make the assump- 
tion that 





ks/k4> 100. 


Since we wish to make sure of having a maximum 
value for A;, an “orientation factor’ of 1/10 
will be allowed for H+HCI as was indicated for 
H+HBr. Then 


ka/ka=(Z18548/27 /0.1Z 240/87) > 100. 


From this is obtained A,—A;>1.4 Kg.Cals. 
Since A, has already been determined to be 
5 Kg.Cals., this gives A3;<3.6 Kg.Cals. 

Another method for obtaining a maximum 
value follows. Bonhoeffer and a co-worker'* tried 
to measure the para- to ortho-hydrogen con- 
version in a mixture of para-hydrogen and 
chlorine when illuminated. Their method was 
accurate enough to determine one percent of 
para- to ortho-conversion, but they were unable 
to detect any change. This means that 


18 Geib and Hartek, Zeits. f. physik. Chemie B15, 116 
(1932). 
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kno 
H +a, ie Heortho tH 


is slow compared with 


k 
H+Cl—>HCI+Cl. 
The assumption may be made then that 


ko /ks<0.01. 
Then we have 


R yo/ka= (Zig 400/B2/Z yo 43/27) <0.01. 


Taking the Z factors in this case to be the 
same, we have 


e~Apo/RT /g—As/RT <).01 


or Aj,o—A2<2.8 Kg.Cals. Since A,» has been 
determined to be in the neighborhood of 7.0 
Kg.Cals.!° a maximum value for A; of 4.2 
Kg.Cals. is obtained. 

Since a minimum value for A; of 1.6 to 2.9 
Kg.Cals. has been obtained and a maximum 
value of 3.6 to 4.2 Kg.Cals., a value, A3(H+Cle) 
= 2-3 Kg.Cals., is indicated. 

It is now necessary to return to a consideration 
of the chain-ending mechanism in the presence 
of oxygen, for which 


H+M-+0,—HO.+M (9) 


has been used. Various experiments have quite 
definitely established that the chain-ending 
mechanism must involve a reaction of hydrogen 
atoms with oxygen rather than chlorine atoms 
with oxygen. Thus, Schwab and Friess”’ passed 
atomic chlorine into oxygen and found no 
measurable reaction whatever. Further, the 
logical product of a reaction of chlorine atoms 
with oxygen, ClOz, is known to be a very strong 
inhibitor for the reaction,”! so that if this product 
were formed, the rate would decrease rapidly. 
No such result has ever been observed. Thirdly, 
a study of yields of HO and HCl in very high 
oxygen-containing mixtures led Krauskopf and 
Rollefson"’ to the conclusion that (Cl+O,) could 
play only a very small part, if any, in stopping 


19 Geib and Hartek, Zeits. f. physik. Chemie, Bodenstein- 
band 849 (1931). 

20 Schwab and Friess, Zeits. f. Electrochim. 39, 586 (1933). 
900)" and MacMahon, J. Chem. Soc. 97, 845 
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the chains. Lastly the types of kinetic expressions 
obtained experimentally indicate that oxygen 
breaks the chains by reaction with hydrogen 
atoms. 

Bodenstein and Schenck” first showed that 
the reaction of hydrogen atoms with oxygen was 
logically a three-body process. They observed 
that although the activation energy of the atomic 
H and O: reaction was 2 Kg.Cals. less than that 
of H and Cl, and consequently should proceed 
50 times as fast as the latter, nevertheless com- 
parative yields of HCl and H,O in H,—O.—Cl. 
mixtures containing equal amounts of Cl, and 
O» showed the rate of H+Cl,. to be twenty times 
faster than H+O.. This meant that only 1 in 
1000 (20 times 50) of the bimolecular collisions 
with sufficient activation energy between hydro- 
gen and oxygen was effective. Since 1 : 1000 is 
approximately the ratio of trimolecular to bi- 
molecular collisions, it appeared that H+O, 
must involve a third body. 

This position has received confirmation from 
Bates** who shows quite conclusively that the 
chain-ending reaction is a three-body process. 
Cook and Bates?* studied the inhibition of the 
photodecomposition of HI in presence of oxygen, 
which involves (H+O,). Considering the re- 
actions 


ka 
H+HI—H-+I, (a) 
ke 


they found that using a bimolecular reaction for 
(8), Rka/kg was about 10 at 100 mm HI and 
100 mm Os, but decreased to 3 when 300—400 
mm of Ne was added. If a trimolecular reaction 
were employed, 


k / 
H+M-+0, — HO.+M (’) 


then the ratio k,/kg,’ remained constant through- 
out. Thus the chain-ending mechanism seems 
quite definitely to be a three-body process. 

The theory on which the values in this 
reaction chain are based is that developed by 
Bodenstein and various co-workers. Some recent 


22 Bodenstein and Schenck, Zeits. f. physik. Chemie 
B20, 420 (1933). 

23 Bates and Lavin, J. Am. Chem. Soc. 55, 81 (1933); 
Cook and Bates, J. Am. Chem. Soc. 57, 1775 (1935). 
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researches, notably the extensive ones of Norrish 
and Ritchie,*4 are at variance with this scheme. 
Thus, while Bodenstein had found a variation 
of the rate of reaction in oxygen-free mixtures 
depending on the first power of the light absorbed, 
Norrish and Ritchie found a dependence on the 
0.6 power. This is explained on the basis that 
the chain-ending mechanism is a recombination 
of chlorine atoms, rather than a reaction of the 
type 
Cl+X—-ClX. 


Thus, either chain-breaking mechanism may 
be correct, depending on the experimental condi- 
tions. This is confirmed by Potts and Rollefson,!® 
who found a first power dependence at room 
temperature and a 0.5 power dependence at 
liquid-air temperatures. Since Cl+He occurs 1 
in 10° collisions at room temperature and only 1 
in 10° at liquid-air temperatures, [Cl] will be 
enormously increased and thus 


Cl+Cl+M-—Clh.+M 


the rate of which depends on the square of the 
chlorine concentration, will be greatly favored at 
liquid-air temperatures as compared to 


Cl+X—CIX 


the rate of which depends only on the first power 
of the chlorine concentration. However, either 
chain-breaking mechanism leads to the same 
type of variation with temperature and allows 
the same conclusions to be drawn concerning 
the activation energy of (Cl+He2). The Boden- 
stein mechanism has been used because Hertel,'! 
and Potts and Rollefson,!° whose values for the 
temperature coefficient were used, both found a 
first-power dependence on light absorbed. 

A second difference is more serious. Bodenstein 
and Unger® found no inhibition by HCl in the 
absence of oxygen, and Bodenstein and Schenck” 
found none in the presence of oxygen. Norrish 
and Ritchie*‘ find serious retardation by that 
substance in both cases. They ascribe this to the 
reaction H+HCl—H:2+Cl which they consider 
to be comparable in rate to H+Cl—HCI+Cl. 
This is not only opposed to Bodenstein’s mechan- 
ism, but also is inconsistent with their own 


*4 Norrish and Ritchie, Proc. Roy. Soc. A140, 112, 713 
(1933), 
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results. For, from thermodynamic considerations 


(Cl+H2) and (H+HCI) must have approxi- 


mately similar rates!® (AE,)=1.0 cals.). On 


Norrish and Ritchie’s assumptions, therefore 


k k 
CI+H.-HCI+H and H+Clh—-HCI+Cl 


have rate constants of the same order of magni- 
tude. Thus, in an equimolecular reaction sys- 
tem, since these two reactions must be pro- 
ceeding at the same rate in order to maintain 
steady state atomic concentrations 


kel Cl JL He ]=kslH JLCle], 
(Cl]j=[H ]. 


If this be true, then CI+CI+M--Cl.4+M, 
H+Cl+M-—HCI+M, and H+H+M-—H.2+M 
should all be effective in breaking chains. Yet 
neither of the latter two reactions enter into 
Norrish and Ritchie’s kinetic expression and they 
state that a mechanism which considers these as 
chain-breaking reactions does not fit the experi- 
mental facts. Either both 


H+HCI—H:2+Cl 
H+H+M—H:+M 
H+Cl+M—HCI+M 


whence 


and 


must enter the kinetic expression or neither can 
enter. The evidence indicates the latter. 

So, although Norrish and Ritchie have un- 
doubtedly obtained excellent experimental re- 
sults, their explanation of them does not hold. 
On the other hand, their results in the presence 
of oxygen are easily explained by adopting the 
three-body process (H+M-+O,) as the chain- 
breaking mechanism, where HCl is an effective 
third body. Their results in the absence of oxygen 
are only partly explained by the trimolecular 
recombination of chlorine atoms. However, some 
light has recently been thrown on this point by 
the experiments of Rodebush and Spealman,?® 
which seem to indicate that H+HCI may occur 
rapidly as a catalytic reaction on the surface of 
Pyrex. If the walls of Norrish and Ritchie’s 
vessels were active in this way, their results are 
easily accounted for. In view of these facts the 


*5 Rodebush and Spealman, J. Am. Chem. Soc. 57, 1040 
(1935). 
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Bodenstein mechanism seems to be the correct 
one to use to obtain the rates of the homogeneous 
reactions. 


HYDROGEN-IODINE REACTION 


Work on the atomic reactions in this case is 
quite scanty since the ordinary combination is a 
molecular process. However, there is sufficient 
evidence to obtain approximate values. 

From his studies of the photodecomposition of 
hydrogen iodide in presence of oxygen, Bates?® 
concludes that for the reactions 


k 
H+HI—H.+I1, (1) 

ky 
H+0.+M—HO.+M, (9) 


k; is about ten times ky when ky is regarded as a 
bimolecular constant. Later work indicated that 
reaction (9) is in fact trimolecular, as it is written 
above. Since the ratio of bimolecular to tri- 
molecular collisions is of the order of 1000, k; 
would be 1000 times ky if every collision were 
effective in both reactions. Allowing an orienta- 
tion factor of 10 for reaction (1), as was found 
for the analogous reaction H+HBr—H:+Br, 
this is reduced to a factor of 100, which is to be 
compared to the observed ratio of 10. One 
concludes that reaction (1) may require an 
activation energy, Ai(H+HI) ~1 Kg.Cal. 

The value for the reverse of this reaction may 
be found from thermal data. Beutler’s‘ value 
for hydrogen was 

H.—2H, AE, = 102.68 Kg.Cals. 
Landolt-Bérnstein”® give 


H.+1L<22HI, AK)=—2.8 Kg.Cals. 


The value for the dissociation of iodine is taken 


26 Landolt-Bérnstein, Tabellen, Erster Erganzungsband, 
p. 809. 
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from Brown?? 


I,—2I, 


AF) = 35.4 Kg.Cals. 
These lead to a value 


H+HI@H2+1, AE, = —32.3 Kg.Cals. 


The activation energy of the reverse reaction is 
obtained by addition and gives A3(I+He2) =33 
Kg.Cals. 

From the observation of Farkas and Bon- 
hoeffer?® that 15 percent of iodine completely 
inhibits the photodecomposition of HI they 
conclude that 


k 
H+I.—HI+1 (2) 


occurs 100 times as frequently as reaction (1) or 
therefore at practically every collision. This 
means that A,(H+I,) ~0. 


COMPARISON OF VALUES 


A comparison of the results obtained by this 
study with those obtained by theoretical calcu- 
lation?® by Eyring’s method by using a 14 per- 
cent Coulombic term is shown in Table I. The 


7 


TABLE I. 








H+Cle H+Bre H+I2 H+HCl H+HBr H+HI Cl+He Br+He I+He 





Experi- 
mental 3 1.0 0 5 1 1 6.0 17.7 33 
Eyring 2.4 2.2 14 11.8 6.6 7.6 12.8 24.0 40 








results agree quite well as to order, but the 
calculated values tend to run higher. It may 
be noted that a larger Coulombic term would lead 
to better agreement in the majority of cases. 


27 Brown, Phys. Rev. 38, 709 (1931). 
28 Farkas and Bonhoeffer, Zeits. f. physik. Chemie 132, 
235 (1928). 
29 Private communication, H. Eyring and A Wheeler. 
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The Group Relation Between the Mulliken and Slater-Pauling Theories 
of Valence 


J. H. VAN VLECK, Harvard University 
(Received October 7, 1935) 


By means of the group theory of characters, it is shown that there is an intimate relation 
between Mulliken’s molecular orbitals and the Slater-Pauling directed wave functions. One 
can pass from the former to the latter by making a simple transformation from an irreducible 
to a reducible representation. Consequently the same formal valence rules are usually given 
by either method, and one can understand generally why wave functions of the central atom 
which are nonbonding in Mulliken’s procedure are likewise never employed in constructing 


Pauling’s “hybridized”’ linear combinations. 
ing’s “hy 





WO distinct viewpoints have been particu- 
larly developed in applying quantum 
mechanics to problems of valence,' namely, the 
Heitler-London-Pauling-Slater method of the 
electron pair bond, and the method of molecular 
orbitals used by Hund, Lennard-Jones, Mulliken, 
and others. The two procedures represent 
different approximations to the solution of a 
complicated secular equation. The method of 
molecular orbitals permits factorization into 
one-electron problems, but at the expense of 
adequate cognizance of the terms due to electron 
repulsion, which are too fully recognized in the 
H-L-P-S procedure. A characteristic feature of 
the latter is the “hybridization,”’ whereby linear 
combinations of states of different azimuthal 
quantum number for the central atom are neces- 
sary in fields of, for example, tetrahedral sym- 
metry. It was shown by the writer that in the 
case of carbon compounds the two theories, 
though superficially different, predicted similar 
results on geometrical arrangement.’ It is the 
purpose of the present paper to show that the 
equivalence is general in the sense that one 
formulation will give the same formal stereo- 
chemical v4lence principles as the other. Thus 
it is futile to discuss whether the Mulliken or 


1L. Pauling, J. Am. Chem. Soc. 53, 1367 (1931); J. C. 
Slater, Phys. Rev. 38, 1109 (1931); R. S. Mulliken, Phys. 
Rev. 40, 55; 41, 49, 751; 43, 279 (1932-3); J. Chem. Phys. 
1, 492 (1933); 3, 375, 506 (1935). For other references, or 
more detailed introduction on the methods which we com- 
pare, see J. H. Van Vleck and A. Sherman, Rev. Mod. 
Phys. 7, 167 (1935). Other workers besides Mulliken have 
contributed to the molecular orbital procedure, but we 
sometimes refer to the latter as Mulliken’s method, since 
we are concerned with the application to polyatomic mole- 
cules in the light of symmetry groups, an aspect considered 
primarily by Mulliken. 

2 J. H. Van Vleck, J. Chem. Phys. 1, 219 (1933). 
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Pauling theory will give better working rules in 
compounds formally amenable to electron pair 
treatment. Both, for instance, suggest the 
now classic Pauling square configuration for 
Ni(CN),- ~ in view of the diamagnetism of this 
ion. (A tetrahedral model would give para- 
magnetism with either method.*) One can only 
inquire which procedure involves the more 
reasonable hypotheses. It seems to us that in the 
case of the transition elements, one must prob- 
ably decide in favor of the Mulliken formation 
as a simple qualitative description, though 
perhaps a poor quantitative approximation. It 
is difficult to believe, for instance, that the 
Fe(CN),* radical has the Pauling structure 
Fe*-(CN)., since the Fe ion certainly is unwilling 
to swallow four extra electrons. The conventional 
ionic model Fet+(CN-)., on the other hand, 
probably goes too far in the other direction. 
The Mulliken viewpoint has here the advantage 
of allowing an arbitrary distribution of charge 
between Fe and (CN)¢, depending on how one 
weights the various atomic orbitals in forming a 
molecular orbital as a linear combination of them. 
Only a limited significance should, however, be 
given to any purported preference between the 
two methods, as each represents a solution of the 
secular equation only under certain extreme 
conditions. The true wave function is in reality 
a combination of H-LJ-M and H-L-P-S func- 
tions, along with many ingredients intermediate 
between these two extremes, and so either theory 
is bound to have some semblance of truth. The 
latter functions, for instance, can be amplified 


Cf. Pauling, reference 1, and J. H. Van Vleck and A. 
Sherman, Rev. Mod. Phys. 7, 206, 221 (1935). 
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by taking linear combinations with terms repre- 
senting different stages of charge transfer until 
finally just the right polarity is obtained. 


MOLECULAR ORBITALS 


We shall confine our discussion to the case 
where a central atom attaches m atoms arranged 
in some symmetrical fashion (tetrahedron, 
square, etc.) characteristic of a crystallographic 
point group.‘ Let ¥(T) be a wave function of the 
central atom which has the proper symmetry, 
i.e., whose transformation scheme under the 
covering operations of the group is that charac- 
teristic of some irreducible representation I. 
Let y; be a wave function of attached atom 7. 
We shall assume that only one orbital state need 
be considered for each attached atom, and that 
this state is either an s state or else is symmetric 
(as in a 20 bond) about the line joining the 
attached to the central atom. The method of 
molecular orbitals in its simplest form® seeks to 
construct solutions of the form 


V=V(T) +2 ai. (1) 


The coefficients a; must be so chosen that 2a; 
transforms in the fashion appropriate to the 
irreducible representation I’. Now the important 
point is that bases for only certain irreducible 
representations can be constructed out of linear 
combinations of the y;. To determine which, one 
ascertains the group characters associated with 
the transformation scheme, usually reducible, 
of the original attached wave functions y; before 
linear combinations are taken. This step is easy, 
as the character xp for a covering operation D 
is simply equal to g, where qg is the number of 
atoms left invariant by D. This result is true 
inasmuch as D leaves g of the atoms alone, and 
completely rearranges the others, so that the 
diagonal sum involved in the character will 
contain unity g times, and will have zeros for the 
other entries. The scheme for evaluating the 
characters is reminiscent of that in the group 


4See R. S. Mulliken, Phys. Rev. 43, 279 (1933) or refer- 
ences 6 and 7 if further background is desired on the aspects 
of group theory and crystallographic symmetry which we 


use. 

5 Called by Mulliken the LCAO (‘‘linear combination of 
atomic orbitals”) form. For a critique of this type of ap- 
proximation see R. S. Mulliken, J. Chem. Phys. 3, 375 


(1935). It appears to have been first suggested by Len- 
nard-Jones. 
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theory of molecular vibrations employed by 
Wigner and by E. B. Wilson, Jr.,° but is simpler 
since we are not interested in displacements of 
atoms from equilibrium, and in consequence all 
nonvanishing entries are unity rather than some 
root of unity. After the characters have been 
found, the determination of the constituent 
irreducible representations proceeds in the usual 
way by means of the theorem that the unre- 
solved characters must equal the sum of the 
primitive characters contained therein. 

As an illustration, we may consider a complex 
containing six atoms octahedrally arranged, i.e., 
located at the centers of the six cube faces. Then 
the symmetry group is the cubic one O,, and the 
characters associated with the arrangement of 
attached atoms are 


E CG. Cy CG! C3 I 1G, IC, IC,’ ICs 
xm™6220004 0 2 6 £=(@) 


Here x(C,), for instance, means the character for 
the covering operation consisting of rotation 
about one of the fourfold or principal cubic axes 
(normals to cube faces) by 27/4. Any rotation 
about such an axis leaves two atoms invariant, 
and hence x(C2) = x(C4) =2. On the other hand, 
x(C2’) =x(C3) =0 since no atoms are left in- 
variant under rotations about the twofold or 
secondary cubic axes (surface diagonals) or 
about the threefold axes (body diagonals). In- 
version in the center of symmetry is denoted by 
I. By using tables of characters for the group O,, 
one finds that the irreducible representations 
contained in the character scheme (2) are, in 
Mulliken’s notation,‘ 


Aig, E,, Ite (3) 


The irreducible representations corresponding 
to various kinds of central orbitals are shown 
below: 


orbit s pb dy de f 
Ai, Tw E, Tou es Tu, “a (4) 





rep. 


The notation for the various kinds of d wave 
functions is that of Bethe,’ viz., 


6 E. Wigner, Gott. Nachr., p. 133 (1930); E. B. Wilson, 
J 3a)) Chem. Phys. 2, 432 (1934); Phys. Rev. 45, 706 
1934). 
7H. Bethe, Ann. d. Phvsik 3, 165 (1929). 
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TABLE I. 














SYMMETRY CENTRAL ORBITALS 
d 


ATTACHED ORBITALS 
f No. REPRESENTATIONS 





Tetrahedral (74) y, E(dy), T2(de) 
Trigonal (D3,) y, Ay te 


Tetragonal (Dn) 4 Aig, Buy, Bag, E, 


A, Ti, T2 
A,’, Az’, A2”, EB’, BE” 


Agu, Biu, Bou, 2E 3 > oF 
Quy luy Quy u Ay, Aou, Big, Bax, | E, 








W(dy1) = (1/12)4f(r)(322—P’), (5) 
¥(dy2) = 2f(r)(?—y"), 


Y(de1)=f(r)xy, (de) =f(r)xz, (6) 
¥(des) =f(r) yz. 


By comparison of (3) and (4) we see that/is it 
impossible for de orbitals to form any partner- 
ships with attached orbitals, in agreement with 
Mulliken’s conclusion® that dy rather than de 
is particularly adapted to forming octahedral 
bonds. Mulliken’s arguments were mainly of a 
qualitative nature. The preceding considera- 
tions enable us to formulate the situation more 
succinctly, as they show that the de orbitals are 
entirely nonbonding. 

In Table I we give the irreducible representa- 
tions, in Mulliken’s notation, contained in the 
central and attached orbitals for compounds of 
other types of symmetry. When 3 or 4 atoms are 
trigonally or tetragonally attached, we have 
supposed that the plane of these atoms is a 
plane of symmetry, as in (NO;)~ or Ni(CN),~~. 
When there is no such symmetry plane, as in 
NHs, the distinctions between u and g, or be- 
tween primes and double primes, are to be 
abolished,? and the symmetries degenerate to 
Cx, Cy instead of D3,, Dan. When 6 atoms are 
attached in the scheme D3,, or 8 in D4, they are 
arranged respectively at the corners of a trigonal 
and a square prism. 

The results given in Table I are obtained by 
the same method as in the octahedral example. 
The explicit forms of the linear combinations of 
the attached orbitals which transform _ir- 
reducibly, or in other words the values of the 
coefficients a; in (1) have been tabulated by Van 


8 R. S. Mulliken, Phys. Rev. 40, 55 (1932). 

®In adapting Table I to the case C3y, the following ir- 
regularity, however, is to be noted: one must replace Aj’, 
A,”, Ao’, Ae”, respectively, by A1, Az, Ao, Ai rather than 
by A1, Ai, Az, Az as one would guess. 


Vleck and Sherman” in many instances, and so 
need not be repeated here. The a; for the octa- 
hedral case are also given in Eqs. (2)—(7) of the 
following paper. 

Note particularly that in the tetrahedral 
complexes, the de orbitals of the central atom are 
bonding, as there are attached orbitals of similar 
group properties with which they can combine, 
while the dy orbitals are nonbonding. The re- 
verse was true of octahedral compounds—a re- 
sult at first a little surprising in view of the 
isomorphism of the groups Ta and O,. This re- 
versal was also deduced by Mulliken*® from the 
geometrical study of the way the central wave 
functions “‘overlap.” 

It will be observed that if eight atoms are 
attached, their full bonding power is not utilized 
unless one includes f wave functions for the 
central atom, since the representation Bo, is not 
included in s, p, or d. Now f wave functions 
usually have too high energy to be normally 
available, or else are so sequestered in the in- 
terior of the atom as to be of no value for bonding 
because of small overlapping. Even if eight atoms 
are attached at the corners of a cube, central f 
wave functions must be included in order to 
realize all possible bonding partnerships, for 
results always true of tetragonal symmetry 
surely apply to cubic symmetry, which is a 
special case of the latter. On the other hand, no 
f functions are needed for six atoms attached 
either octahedrally or at the corners of a trigonal 
prism. We thus have an indication of why it is 
that coordination numbers of six are common in 
nature, while those of eight are rare. 


METHOD OF DIRECTED ELECTRON PAIRS 


We now turn to the method of Pauling and 
Slater. Here the procedure is to use hybridized 


10 J. H. Van Vleck and A. Sherman, Rev. Mod. Phys. 7, 
219 (1935). 
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central orbitals, i.e., linear combinations of 
orbitals of different azimuthal quantum number, 
in such a way that the resulting central wave 
function projects out especially in some one 
direction in space, and so is adapted to form an 
electron pair with one particular attached atom. 
Thus for tetrahedral compounds Pauling and 
Slater! use wave functions which are linear 
combinations of s and #, or alternatively as 
Pauling! shows, of s and de wave functions. 
For octahedral compounds, Pauling finds that 
sp'dy? combinations are appropriate, and sp*dy 
for tetragonal. Hultgren! proves that eight 
atoms cannot be attached (at least symmetri- 
cally) by means of unidirectional electron pair 
bonds formed from s, p, and d wave functions. 
Incidentally, the present paper shows that 
sp*d*f functions are needed to hold eight atoms.” 
It will be noted that the wave functions in- 
volved in the Pauling unidirectional linear 
combinations are precisely those which are 
bonding in the method of molecular orbitals. 
For example, Pauling, like Mulliken, makes no 
use of dy orbitals for tetrahedral compounds, or 
of de for octahedral. Such coincidences have 
hitherto appeared something of a mystery, but 
asimmediate explanation, as follows, is furnished 
by group theory. 

In the Pauling-Slater theory, one desires the 
central wave functions to possess unilateral 
directional properties so as to be correlated with 
one particular attached atom. Hence the P-S 
central functions must have the same transforma- 
tion properties as do those y; of the attached 
atoms before linear combinations of the latter 
are taken. Thus the problem of finding the 
linear combinations of the central orbitals which 
exhibit the proper directional properties is 
simply the reverse of finding the proper linear 
combinations of the attached orbitals in the 
Mulliken procedure. The difference is only that 
in the P-S theory, the linear combinations are in 
the central rather than attached portion, and 
their construction corresponds to transforma- 
tion from an irreducible representation to a 


1 R, Hultgren, Phys. Rev. 40, 891 (1932). 

12 Similar conclusions on the type of bonds necessary to 
attach eight atoms have also been obtained in unpublished 
work of R. S. Mulliken. 
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reducible one, of structure similar to that be- 
longing to the original y,’s, rather than to the 
inverse transformation. Clearly, the same irre- 
ducible representations are needed in the con- 
struction of a given reducible representation as 
those contained in the resolution of the latter 
into its irreducible parts. Pauling has obviously 
shown considerable ingenuity in constructing his 
wave functions without using the status in terms 
of group theory. 


GENERALITY OF THE RESULTS 


The argument underlying Table I, etc., os- 
tensibly assumed that the molecular orbital be 
expressible as a linear combination of atomic 
orbitals, but is readily seen to be still applicable 
provided only that the charge cloud of any at- 
tached orbital be symmetric about the line 
joining the given attached atom to the central 
one. Hence the atomic orbitals can be of what 
James calls the flexible type, i.e., contain parame- 
ters which can be varied in the Ritz method, 
and which allow for the fact that chemical com- 
bination distorts the atomic orbitals from what 
they would be in the free condition. This admis- 
sion of flexibility is fortunate, for it is well 
known that it, is a bad quantitative approxi- 
mation® to express a molecular orbital as a 
linear combination of undistorted atomic or- 
bitals. The Ritz variational problem is, of course, 
to be of the 1 rather than electron type, so that 
the generality in our analysis by means of 
molecular orbitals is roughly comparable with 
that in the Hartree method. 

One thing which the preceding analysis does 
not do is to tell us what is the best arrangement 
of atoms in case the symmetry group does not 
uniquely determine this arrangement. For in- 
stance, by examining the overlapping of wave 
functions, Hultgren" finds that when six atoms 
are attached at the corners of a trigonal prism, 
the binding is firmest if the sides of the prism are 
square. This fact cannot, however, be inferred 
from our group theory considerations, as there 
are no additional elements of symmetry when 
the sides are square rather than rectangular. 


13 This transformation has been explicitly given by the 
writer in the case of methane (J. Chem. Phys. 1, 177 (1933)), 
but he did not discuss its group-theoretical significance. 
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Certain complex salts, notably ferro- and ferricyanides, 
have susceptibilities much lower than those predicted by 
the Bose-Stoner “‘spin only” formula. The first interpreta- 
tion was that given by Pauling on the basis of (I) directed 
wave functions. In the present paper it is shown that 
alternative explanations are possible with (II) the crystal- 
line potential model of Schlapp and Penney, or with (IIT) 
Mulliken’s method of molecular orbitals. In any of the 
theories, the interatomic forces, if sufficiently large, will 
disrupt the Russell-Saunders coupling, and make the 
deepest state have a smaller spin, and hence smaller 
susceptibility, than that given by the Hund rule. This 


situation is not to be confused with that in normal para- 
magnetic salts, such as sulphates or fluorides, where only 
the spin-orbit coupling is destroyed. The similarity of the 
predictions with all three theories is comforting, since any 
one method in valence usually involves rather questionable 
approximations. Because of this similarity, a preference 
between the theories cannot be established merely from 
ability to interpret the anomalously low magnetism of the 
cyanides. Covalent bonds, as in cyanides, seem to be more 
effective in suppressing magnetism than are ionic ones, 
as in fluorides, but so far the evidence to this effect is 
empirical rather than theoretical. 





N most salts of the iron group, the suscepti- 
bility has approximately the value 


x=4NBS(S+1)/3kT, (1) 


where 8 is the Bohr magneton he/4rmc, and 
where the spin S has the value given by the 
Hund rule that the ground state is that of 
maximum multiplicity compatible with the Pauli 
principle. Formula (1) was suggested by Bose 
and Stoner. The theoretical explanation has been 
known for some time,' and is that the crystalline 
field is so strong as to destroy the coupling of 
the spin angular momentum S and orbital 
angular momentum L to a quantized resultant J. 
It can be shown that with such an “electric 
Paschen-Back effect,”’ the orbital magnetism is 
largely quenched, and the susceptibility has 
nearly the value (1). 

There are, however, certain compounds of 
elements of the iron group, notably iron cyanides 
and the various cobaltammines, for which the 
susceptibility has a value very much lower than 
that given by (1). In fact, these compounds are 
diamagnetic if they involve a complex ion con- 
taining an even number of electrons (e.g., 
Fe(CN).*-), or have a susceptibility of an order 
of magnitude corresponding to one free spin, i.e., 
to S=4 in (1), if this number is odd (e.g., 
Fe(CN),*-). The first specific model accounting 
for this behavior was given by Pauling,” but the 


1J. H. Van Vleck, The Theory of Electric and Magnetic 
Susceptibilities, Chap. XI. 

?L. Pauling, J. Am. Chem. Soc. 53, 1367 (1931); L. 
Pauling and M. L. Huggins, Zeits. f. Krist. 87, 205 (1934). 


explanation may be couched in more general 
language as follows: the interatomic forces are so 
very large as to destroy not merely the spin-orbit 
but also the Russell-Saunders coupling, and make 
the deepest state that of lowest possible spin 
rather than of maximum spin as given by the 
Hund rule.* Large spin would be an advantage 
as far as a free atom is concerned, but the point 
which we wish to make is that the interatomic 
energy is decreased by lowering the total spin 
regardless of whether one makes the calculation 
in any of three ways: 


(I) the directed electron pair bond, some- 
times called the Heitler-London-Pauling-Slater 
approximation, 

(II) the method of crystalline fields, particu- 
larly adapted to ideally ionic compounds, and 

(III) Mulliken’s method of molecular orbit- 
als.4 


There is thus a formal similarity between the 
results of various approaches, in many respects 
paralleling that on the subject of directional 
valence stressed in the preceding paper. The 
cyanide anomalies have been previously inter- 
preted only on the basis of I, although Mulliken‘ 
did state that magnetic behavior might often be 
equally intelligible with the mechanism III. 

The particular examples which we shall con- 
sider are compounds in which the central or 


3 For this rule see, for instance, Pauling and Goudsmit, 
The Structure of Line Spectra, p. 165. 
*R.S. Mulliken, Phys. Rev. 40, 55 (1932). 
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paramagnetic atom has a coordination number of 
six resulting from six neighboring atoms or ions 
octahedrally arranged. The clusters Fe(CN).*-, 
Fe(CN),*-, FeF.*-, in KyFe(CN)¢, KsFe(CN)<, 
(NH,)3FeF¢, respectively, are structures of this 
character. It is not necessary that the crystal as 
a whole have cubic symmetry. For instance, 
K;Fe(CN). is monoclinic. From these examples 
we trust it will be sufficiently clear how the 
similarity between the various methods is demon- 
strated in other cases. 

The deepest state of the Fet** ion is a ®S 
term, and indeed the observed susceptibility® of 
(NH,)3FeF, approximately equals that given by 
(1) with S=5/2. The susceptibility® of K;Fe- 
(CN)., however, is much nearer the value given 
by (1) with S=4. The interatomic forces should 
thus be relatively more important in the cyanides 
than in the fluorides, and it is, in fact, known 
that the former are firmer compounds than the 
latter. 


I. THE METHOD OF DIRECTED ELECTRON PAIR 
Bonpbs 


The magnetic behavior of Fe(CN),.*~ is ex- 
plained by Pauling on the basis of the structure 
Fe*-(CN),. This model need not be literally true, 
but is taken to be sufficiently typical, ie., a 
sufficiently common phase in the resonance 
through various stages of polarity, to serve for 
purposes of discussion. According to Pauling,? 
the Fe*~ ion attaches its six atoms of coordination 
by means of six sp*dy? electron pair bonds. 
These six bonds consume the 4s, 4p, and dy 
states,’ and house six electrons of Fe*-. (Only 
one electron can be assigned to a state of the 
central atom which is used for bonding purposes ; 
for if a state houses two electrons, their spins 
compensate each other, and they are not free to 
form electron pair bonds with outside atoms.) 
The Fe* ion has 11 electrons in all to house 
(apart from completed inner shells). The five 
electrons not absorbed by the electron pair bonds 
can only be accommodated in the three 3de 
states. Thus two of these states will be filled 






5 E. Cotton-Feytis, Ann. Chim. 4, 9 (1925). 
6 Cf. for instance, L. C. Jackson, Proc. Roy. Soc. A140, 
695 (1933). 

7 For explanation of the notation dy, de see Eqs. (5)—(6) 
of the preceding paper in this issue. 
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twice. The electrons which are thus “doubled 
up” must have their spins antiparallel because 
of the Pauli principle, and hence make no 
contribution to the susceptibility. The fifth 
electron, however, has a free spin, so that the 
susceptibility should be given approximately by 
(1) with S=} provided one can overlook orbital 
contributions to the susceptibility. Actually they 
are not negligible, as Howard shows in the next 
paper. In Fe(CN)¢.4, there is an additional 
electron to house, so that the fifth electron loses 
its private state, and there can only be dia- 
magnetism, in agreement with experiment. 


II. THE METHOD OF THE CRYSTALLINE 
POTENTIAL 


The simplest model for explaining the varying 
magnetic properties is that used with success by 
Penney and Schlapp® primarily on hydrated 
sulphates of the rare earth and iron groups. 
The basic assumption is that the effect of 
neighboring atoms can be represented by means 
of a static potential. One may distinguish be- 
tween three cases depending on the strength of 
the crystalline field: (a) the field is so weak that 
the inner quantum number J has a meaning 
(b) it is so largé as to prevent J, but not L or S 
from being a good quantum number, and (c) it 
is still more powerful, and able to destroy 
Russell-Saunders coupling, i.e., the compounding 
of the individual /’s to a resultant L, so that L 
loses its validity as a quantum number. Cases 
(a) and (b) are those studied by Penney and 
Schlapp, (a) being characteristic of the rare 
earths, and (b) of the hydrated sulphates of the 
iron group. We shall now show that (c) furnishes 
an explanation of the behavior of the cyanides 
and other salts having abnormally low suscepti- 
bilities. 

In the compounds in which we are at present 
interested, the crystalline field is of dominantly 
cubic symmetry. In a field of this symmetry 
type, it can be shown by means of group theory,’ 
or otherwise, that a d electron has the Stark 
pattern shown in Fig. 1. The separation of the 


8W. G. Penney and R. Schlapp, Phys. Rev. 41, 194 
(1932); R. Schlapp and W. G. Penney, ibid. 42, 666 (1932); 
O. Jordahl, ibid. 45, 87 (1934); R. Janes, ibid. 48, 78 (1935). 

9 See, for instance, H. Bethe, Ann. d. Physik 3, 143 (1929). 
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cubic nearly cubic 


Splitting in method of 
molecular orbitals 


Fic. 2. 


Splitting in method of 
crystalline potentials 


Fic. 1. 


levels which branch out from a common origin 
on the left is, by hypothesis, comparatively 
small, and owes its existence to the deviations of 
the field from perfect cubic symmetry. Do not 
confuse Fig. 1, which is for one electron, and 
the closely related crystalline Stark diagrams 
given by Penney, Schlapp, and Van Vleck for 
whole atoms in which the electrons are space 
quantized collectively rather than individually. 
Fig. 1 is right side up, as Gorter!® shows, so long 
as the negative ions surrounding the para- 
magnetic cation are octahedrally arranged. 
Fig. 1 should be inverted in case the coordin- 
ation number is four rather than six, i.e., the sym- 
metry tetrahedral rather than octahedral. Now 
if we assume the conventional polar structure 
Fe+++(CN~-)., the Fet+++ ion has five electrons 
outside closed shells, and there are five electrons 
to be housed in the various states of Fig. 1. 
Clearly thedeepest energy isobtained by assigning 
the five electrons to the three components of de. 
Then all the electrons but one double up, and 
there is only one free spin, while in Fe++(CN-), 
there would be no free spin. Thus with a suffi- 
ciently large field to make individual space 
quantization a good approximation, the deepest 
state is quite different from that °S given by the 
Hund rule. 

It is interesting to note the significance of the 
state d5®S in terms of Fig. 1. Usually Russell- 
Saunders coupling involves collective quanti- 
zation in such a way as not to permit any simple 
interpretation in terms of a one-electron diagram 


10 C, J. Gorter, Phys. Rev. 42, 437 (1932). 
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such as Fig. 1. However, when the spin quantum 
number equals half the number of electrons, as 
in d°*S, the spins are all mutually parallel. 
Such an alignment is possible only if each elec- 
tron has its own private orbital. Hence a state 
of the configuration d° with S=5/2 involves one 
electron in each of the five levels of Fig. 1. This 
distribution has an invariant significance (one to 
a state) regardless of how the axis of space 
quantization is chosen. Hence this state must be 
an S state, as is also proved by other means and 
as is reflected in the high magnetic isotropy of 
manganous and ferric salts of the sulphate 
variety. When, however, one has the cyanide 
case, there are three different ways in which the 
five electrons may be distributed among de, 
since the private orbital may be the lower, 
middle, or upper component of de. True cubic 
symmetry is achieved only if the three compo- 
nents are populated equally. Hence K;Fe(CN)¢ 
can exhibit a high degree of magnetic anisotropy. 
This point is studied quantitatively by J. 
Howard in the next paper. 

As we have seen, the energy due to the 
crystalline field is lowest when the Hund rule is 
broken down. However, the internal energy of 
the Fe+++ ion is lower with S=5/2 than with 

=}. (This is the meaning of the Hund rule.) 
Thus whether the total internal plus external 
energy is lower or higher for S=5/2 than for 
S=} is a question purely of the size of the 
crystalline potential relative to the strength of 
the Russell-Saunders coupling. The fact that 
both cases are actually realized suggests that 
the two energies are of the same order of magni- 
tude, and such indeed appears to be the case. 
The crystalline fields deduced by Schlapp, 
Penney, Jordahl, and Janes*® for the hydrated 
sulphates usually amount to about 3 volts." 


In Crt+** an over-all splitting of about 8 volts is ob- 
tained by Schlapp and Penney, a value which seems unduly 
high and out of line with their other results. As they inti- 
mate, the explanation of this discrepancy is probably that 
in Crt+*, the computation of the splitting is unusually 
sensitive to small experimental errors in the determination 
of the absolute value of the susceptibility because the 
latter has very nearly the ‘‘spin-only”’ value in this par- 
ticular ion. Schlapp and Penney employed the Leiden data 
on chrome alum. It is interesting to note that a splitting 
(3 volts) of about the usual size is yielded by Janes’ recent 
measurements on K;Cr(SCN).4H.O. Incidentally, the 
separation 17,200 cm™ which Janes obtains in cupric salts 
relates to the over-all splitting, rather to the constant D as 
stated in his article, and so is not unreasonable. 
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The strength of the Russell-Saunders coupling 
is measured by the separation of the various 
levels arising from the configuration d* in the 
free ion. Spectroscopic data indicate that this 
separation amounts to only about 4 or 5 volts.” 
Hence the crystalline field would not have to be 
much larger than in the examples studied by 
Schlapp and Penney in order to overpower the 
Russell-Saunders coupling,'* and so it is plausible 
that this situation is realized in the cyanides. 


Even in the case studied by Schlapp and Penney, 
labeled (b) above, the separation between the various 
states of the configuration 3d* is not larger than the 
crystalline potential. Consequently one naturally wonders 
whether they were justified in neglecting matrix elements 
between states of different Z but similar S belonging to 
the same configuration. Neglect of matrix elements be- 
tween states differing both in L and S is legitimate if the 
intervals between such states are large compared with the 
spin-orbit interaction, i.e., with multiplet widths. This 
condition is practically always fulfilled. In Fe**+* or Mnt+ 
there can be no question of rigor, since the configuration 
d® has only one sextet state. Also in Cr*++ or Fe** (d‘ or d®) 
there is only one quintet, while in Cu** (d°) there is only 
one state of any character, a doublet. In Crt++* or Cott, 
however, the states of highest multiplicity are *F and ‘P, 
while in Nit*, they are *F and *P. The separation between 
these two states is only about 2 volts! in either Co** or 


12 Spectroscopic data are not available on the Fe*** ion. 
However, the separation of the various terms belonging to 
the configuration d> in Fe*** should be somewhat greater 
than (probably a little less than double) the corresponding 
separations in the homologous ion Cr*, and the interval 
d’ 6S —d 4G, for instance, is known to be 2.5 volts in Cr*. 

‘3 Unfortunately it does not appear possible to estimate 
directly from theoretical considerations the crystalline 
splitting to be expected even with an ideal ionic structure 
and assumed interatomic separations obtained from Paul- 
ing’s atomic radii. So one can only deduce the splittings 
empirically from the magnetic data. The difficulty is that 
one does not know well enough the effective charge Z to be 
used in computing the d wave functions. Clearly Z should 
be somewhat greater than 4, the value corresponding to 
perfect screening, and somewhat less than the effective 
charge Zi, deduced from ionization potentials, as the value 
of Z to be used in computing 7‘ etc. is less than that in- 


volved in 1/r. Extrapolated spectroscopic data indicate 
that Zi, is about 6. In unpublished calculations, Howard 
finds that the interval dy —de amounts to 4 volts if Z=4 
and to 1 volt if Z=6, provided the further assumptions are 
made that the distance Fe—F is 1.91A, and that the F~- 
ions act like point charges. 

14 The intervals d’ 4F—d7 4P of Cot+ and d*3F—d8*P of 
Ni** have not been observed directly, but probably do not 
differ greatly from the intervals d?(4Fyas? — 74P)4s? of Col, 
and d°(3F)4s?—d8(8P)4s? of Ni I, respectively. These latter 
intervals are known and are both 1.9 volts. That the addi- 
tion of the 4s electrons does not change too materially the 
separations of the core states is indicated, for instance, by 
the fact that the frequency difference d*(3F)4s 4F—d4(8G)4s 4G 
in Cr II deviates only 4 percent from the difference 
d**F—d*3G in Cr III (see Bacher and Goudsmit’s tables). 
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Nit*. Spedding" criticizes the work of Schlapp and Penney 
on the ground that they do not consider the perturbing 
effect of *P in their calculations of the *F state of Nit**, 
or of *P on the *F of Crt+*, Cot+. However, Spedding 
has overlooked the fact that in Cr*++* or Nit* the Stark 
component of the F level which is important for the 
magnetic calculations is that which belongs to the repre- 
sentation A2!® of the cubic group. On the other hand, *P 
or ‘P belongs to the representation 7), and so is incapable 
of perturbing Ag, if the field is really cubic. This statement 
is not quite true when one allows for the fact that the 
spin-orbit distortion makes the symmetry different from 
that characteristic of orbit alone. In other words, one 
ought rigorously to consider the irreducible transformation 
properties of the spin-orbit rather than just orbital wave 
functions. However, the spin-orbit distortion is subordi- 
nate, and so the perturbation due to this cause is negligible. 
Furthermore the Stark component Az of °F or 4F involves 
different rhombic representations than do *P or ‘*P, and 
so cannot be perturbed by the latter even when one 
considers the deviations from cubic symmetry, as long as 
rhombic symmetry is preserved. Hence Schlapp and 
Penney’s conclusions on the nearly perfect magnetic 
isotropy of Nit+, Cr**+ and their close conformity to (1), 
are unaffected. In the case of Co**+, the ground level is 
the Stark component 7; of d’*F, and can be perturbed 
by d’74P. We hope to study this effect later more fully. 
In Co**, anyway, the agreement which Schlapp and 
Penney obtained with experiment was qualitative rather 
than quantitative. Rough preliminary examination sug- 
gests that it will probably be improved by considering the 
perturbation by *P. Hence, we conclude that all the vital 
features of Schlapp and Penney’s calculations are un- 
affected by considering the incipient breakdown of Russell- 
Saunders coupling, and that in Cot* the quantitative 
agreement may actually be bettered. 


III. THe METHOD oF MOLECULAR ORBITALS 


This method differs from the crystalline field 
procedure II in that the structural unit for the 
wave function is the whole complex ion (e.g., 
Fe(CN),°-) rather than the single central atom. 
Both II and III utilize one-electron wave func- 
tions, and so group theory based on symmetry 
properties can still be used to obtain information 
about the character of the levels. However, 
there is the difference that if a central orbital is 
bonding, a given representation appears more 
often than in method II because of the fact that, 
in the language of the preceding paper, there are 
both central and attached orbitals belonging to 
the same irreducible representation. The proper 


1° F, H. Spedding and G. C. Nutting, J. Chem. Phys. 2, 
421 (1935). 

16 For explanation of the notation for the group repre- 
sentations see R. S. Mulliken, Phys. Rev. 43, 279 (1933). 
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linear combination of such central and attached 
orbitals is found by solving a secular equation. 
In the case of the octahedral group, the wave 
functions belonging to the irreducible represen- 
sentation E, are 


¥(E,) = ay(3dy1) + (1 —a*)*(1/12)*(2ys 
+ 2Y6—¥i-Ya— 2 — Ys), 


VE) = (1 — a) b(3dy1) — a(1/12)*(2ps 
+2f6—vi- va 2 - Vs), 


¥' (Ey) =ayp(3dy2) 
i (1 ses a®)§3 (1 +yYi-— yo —Ys), 


v'"" (E,) = (1 —a*) ty (3dy2) 
—az(ityi-—ve—Ys), 


where the wave function of attached atom 7 is 
denoted by yi, and the explicit forms of the 
central wave functions ~(dy) are as given in 
Eqs. (5)—-(6) of the preceding paper. Attached 
atoms 1, 4 are supposed located on the x axis; 
2, 5 on the y; and 3, 6 on the z. The value of a 
is determined, at least in principle by solving a 
secular equation. Similarly, the wave functions 
corresponding to A,, and 7}, are 


¥(A1,) =By(4s) 
+(1- 8) 1G) i tvetystyitystye), 


¥(Tin) acs vy (4 poz) + (1 —7)(5)?(fi— Ys), 
W (Tin) = (Abo y) + (1 — 77)4(3)! (2 — Ys), 
W'" (Tin) = (Apo 2) +(1—7)*(2)*(Ws— Ye), (7a) 


together with four other wave functions (4b)— 
(7b) which are similar except that 8, (1— £7)! 
are replaced, respectively, by (1—§*)!, —8, etc. 
The notation po, means a p state with m,=0 
when the quantization is relative to the x axis. 
The basic theoretical principles underlying the 
construction of wave functions such as (2)-(7) 
are outlined in the preceding article." 

The lower roots of the quadratic secular 
equations associated with (2)—(7) represent lower 
energy than for the free atom, i.e., are bonding 
levels, while the upper roots give higher energy 
than for the free atoms, i.e., are antibonding. 


(2a) 


(2b) 


(3a) 


(3b) 


(4a) 
(Sa) 


'7 Cf. also J. H. Van Vleck and A. Sherman, Rev. Mod. 
Phys. 7, 218-222 (1935). 


(6a) 
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The de wave functions are nonbonding, as shown 
in the previous paper. The situation is illustrated 
in Fig. 2. All these statements on the positions 
of the levels are surely true as long as the im- 
portant term for bonding is the so-called Hund 
resonance integral, which is the off-diagonal 
matrix element connecting the central and 
attached portions. Since diagonal sums are 
unaffected by the addition of off-diagonal ele- 
ments, it follows that the sum of the two roots 
of the quadratic must equal an energy level of 
the free central atom, plus one of the attached 
atom. Hence if one root is lowered in virtue of 
the resonance integral, the other must be raised. 
Let us now discuss the case of Fe(CN),*. 
We can consider each CN radical as contributing 
one electron (viz., the electron which corresponds 
to a 20 wave function directed towards the 
central atom). There are thus 17 electrons in all 
to house (apart from completed inner shells). 
Twelve electrons can be housed in the lower 
roots of the secular equations associated with 
(2)-(7), since the Pauli principle allows two 
electrons to an orbital state. The remaining five 
electrons are then to be accommodated in the 
three nonbonding, i.e., purely central de orbitals, 
which are preferable to antibonding states. One, 
and only one, of these five electrons can have a 
private orbital, and so only one spin contributes 
to the susceptibility. With Fe(CN),.‘, the last 
six electrons are snugly accommodated in the 
three de levels, and so there is diamagnetism. 
Of course, if the bonding action is weak enough, 
several electrons may prefer to occupy private 
orbitals, even though antibonding, in order to 
secure large spin and low internal energy for the 
central atom. This is a refinement not included 
in the crude forms of the method of molecular 
orbitals, for it is the interelectronic exchange 
energy, not incorporated in the one-electron 
problem, which is lowered by making the spins 
parallel. Thus we see that, as in the other 
methods, the susceptibility will have the low 
values characteristic of the cyanides only if the 
bonding action is larger than the Russell- 
Saunders structure of the central atom. 


COMPARISON OF THE VARIOUS METHODS 


We have now seen that all three descriptions 
are capable of explaining the anomalously low 
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magnetism of the cyanides. The question natu- 
rally arises as to which corresponds the closest 
to reality. It seems to us that in the cyanides 
one must decide in favor of III or possibly I. 
One reason is that method III is the most 
general, since it is noncommittal on the amount 
of polarity, which can be anything between 
Fe*+(CN-), and Fe®-(CN+)., depending on the 
size of the coefficients a, 8, y. In fact, II can be 
regarded as a special case of III, wherein these 
coefficients are zero and the CN orbitals have 
lower energies than those of the central atom, 
so that the twelve electrons digested by (2-7) are 
located on the CN radicals. Method I is also to 
a certain extent a special case of III, correspond- 
ing to 

a=B=y=V}, (8) 


but we have seen in the preceding paper that in 
any event I differs from III in assuming localized 
or unidirectional bonds. Eq. (8) may be regarded 
as the condition for ideal covalency, as when it 
is satisfied, the twelve bonding electrons spend 
equal amounts of time on the (CN). group and 
on the central atom. The actual behavior is not 
as extreme as in (8), for the structure Fe*-(CN)< 
given by (8) requires too much negative charge 
on the Fe ion, as mentioned at the beginning of 
the preceding paper. 

A decision between II and III is furnished by 
the fact that (NH,)3FeF, obeys the ordinary 
formula (1), andyso has a much higher magnetism 
than KgFe(CN)e. Now it is known from other 
evidence!® that (NH,4)3FeF, is much more polar in 
its structure, i.e., more nearly (NH,*);Fe**(F-)., 
than is K;Fe(CN).. With the method of the 
crystalline potential, the fields responsible for 
the Stark splitting are those due to the neg- 
ative charges surrounding the paramagnetic 
cation. Hence if method II were always appli- 
cable, the quenching of magnetism would be 
more complete in the fluoride than in the 
cyanide, contrary to experiment. On the other 
hand, in the method III of molecular orbitals, 
the strength of the bond or separation of the 
levels stands in no immediate relation to the 
amount of polarity, and it is perfectly possible 
to have a covalent structure quench magnetism 


18 L. Pauling, J. Am. Chem. Soc. 54, 988 (1932); Pauling 
and Huggins, reference 2. 
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more completely than an ionic one. In fact, if 
the resonance integrals are the determining 
factor in chemical bonds, one comes nearer to 
(8) the greater the interatomic forces available 
to suppress the magnetism. 

Pauling has argued that the suppression of 
magnetism is evidence in favor of the electron 
pair mechanism I in Fe(CN),*°-. However, we 
have seen that II or III can also be effective in 
quenching magnetism, so that it does not seem 
possible directly to infer the bond type from 
magnetic behavior. We agree with Pauling that 
covalent bonds are more effective in destroying 
magnetism than ionic, but the reasoning by 
which we reach this conclusion is somewhat 
different and is based on the empirical com- 
parison of (NH,4)3FeFs and K;Fe(CN), made 
above. Once it has been established empirically 
in a few cases that covalent bonds destroy 
magnetism more than ionic, one is probably 
justified in extrapolating to other cases, and 
supposing the bonds to be covalent when most 
of the magnetism has been destroyed. 

If method III is to be preferred to II, the 
calculations of Schlapp and Penney,* based 
ostensibly on II, must be given a somewhat 
different interpretation than previously. They 
by no means lose applicability, for they are 
based largely on the symmetry group properties, 
and so retain practically as much significance 
with III as with II, except that the size of the 
crystalline potential is not to be taken too 
literally. The splittings of a few volts may now 
relate to the magnitude of the Hund resonance 
integrals rather than of the crystalline potential, 
but remain comparable to the Russell-Saunders 
structure, making it still reasonable that one 
should sometimes have the cyanide behavior and 
sometimes the sulphate. 

One feature which is common to methods I, 
II, III is that nickel salts should have a suscepti- 
bility corresponding to S=1 rather than to S=0 
in (1) even if the cubic splitting is large compared 
with the Russell-Saunders structure. This state- 
ment is true if the grouping about the nickel ion 
is either tetrahedral or octahedral. The demon- 
stration has been given by Pauling? with me- 
chanism I and by Van Vleck and Sherman!’ 
with III. If II is used, one has only to note 
that Ni** has two more electrons than Fett, so 
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that there are eight rather than six electrons to 
be housed in the five levels of Fig. 1. Hence 
regardless of whether Fig. 1 is upright or in- 
verted, two electrons can have private orbitals 
without materially diminishing the energy, pro- 
vided the rhombic splitting in Fig. 1 is small. 
To explain the observed diamagnetism of K,Ni- 
(CN)s, one is thus led to predict a square 
(tetragonal) rather than tetrahedral configura- 
tion for the Ni(CN),-~ ion. This prediction is 
confirmed by x-ray measurements!® on BaNi- 
(CN),.4H2O. On the other hand, Ni(NH3),SO, is 
paramagnetic, and this fact suggests that there 
may be a tetrahedral structure for the group 
(NHs)«. 

Our arguments to explain the diamagnetism of 
K,Fe(CN). apply equally well to diamagnetic 
cobaltic compounds such as the cobaltammines, 
since Cot++ is isoelectronic with Fet*+. The 
writer has shown elsewhere” that if the coordi- 
nation number is six, and if (1) is valid (the 
sulphate rather than cyanide case), cobaltous 
compounds should show marked anisotropy and 
deviations from Curie’s law, while the reverse 
should be true of nickel ones. If the coordination 


number is four, the roles of cobalt and nickel 


19H. Brasseur, A. de Rassenfoss and J. Piérard, Zeits. f. 
Krist. 88, 210 (1934). 

20J. H. Van Vleck, Phys. Rev. 41, 208 (1932). In this 
connection see reference 10 regarding the relation between 
coordination number and the sign of the crystalline field. 
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should be interchanged. Blue and pink cobaltous 
salts are supposed to involve coordination num- 
bers of four and six, respectively.”! It would 
therefore be interesting to see if one could 
prepare a blue cobaltous salt conforming closely 
to Curie’s law. Attempts of R. B. Janes to do 
this have so far proved unsuccessful. The ex- 
planation may be either that the magnetic 
dilution is insufficient, so that the exchange forces 
between paramagnetic ions complicate the tem- 
perature dependence, or else that one is in the 
critical transition region wherein the mechanism 
for suppressing the magnetism below (1) just 
comes into play. In the latter event both states 
with S=} and S=3/2 are populated, and the 
temperature variation will not be simple. 

In conclusion, we may say that it should be 
regarded as reassuring and significant that all 
three mechanisms I, II, III permit a rational 
understanding of magnetism lower than that 
given by (1) and the Hund rule. Since none of 
the methods of approximation can be regarded 
as satisfactory from a quantitative standpoint, 
a property which is common to all three must 
be regarded as on a much firmer basis than if it 
is characteristic only of one approach. 

The writer wishes to thank Mr. John Howard 
for interesting and valuable discussions. 


21 R. Hill and O. R. Howell, Phil. Mag. 48, 833 (1924). 
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The magnitudes, anisotropy and temperature depend- 
ence of the principal magnetic susceptibilities of K3Fe(CN)« 
are calculated by the method of crystalline potentials, and 
prove to be in good agreement with experiment. The 
assumption is made that the interatomic forces in the 
Fe(CN)s~—~ complex destroy Russell-Saunders coupling. 
This results in a susceptibility corresponding to one free 
electron spin with a contribution, except at very low 
temperatures, from the unquenched orbital angular mo- 
mentum. The effect of the latter is sufficient to make the 
susceptibility at room temperature almost twice as large 
as the value which would be obtained if the spin alone 
were considered. The superposition of a small rhombic 


crystalline field suffices to produce the large observed 
magnetic anisotropy. The behavior of the susceptibility 
of K;Fe(CN). is in sharp contrast to the higher 
effective magneton numbers and magnetic isotropy of 
many ferric compounds, such as the sulphates. In these 
the ground state of Fe**+ is the normal state of the free 
ion, namely a °S state, so that five spins contribute to the 
susceptibility. Although the calculations ostensibly use the 
method of crystalline potentials, they are really based on 
general group properties and so apply even if the bonds in 
the cyanide are covalent, as they probably are, rather than 
ionic. 
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INTRODUCTION 


HE magnetic measurements made by Jack- 
son! on single crystals of K3;Fe(CN), in the 
temperature range 75°—290°K reveal, in addition 
to the abnormally low susceptibility reported by 
other investigators, a remarkable anisotropy 
amounting to almost 100 percent at low tempera- 
tures. Compounds of ferric iron in general ex- 
hibit a high degree of magnetic isotropy and a 
mean susceptibility corresponding to a spin of 
5/2. This latter behavior is to be expected if the 
Fe*** ion in the crystal is in the ®S state, which 
is the ground state of the free ion.? In the second 
of the foregoing articles’ Van Vleck has demon- 
strated how the low susceptibility and magnetic 
anisotropy of K;Fe(CN). may be explained 
upon the assumption that the interatomic 
forces are sufficiently powerful to destroy 
Russell-Saunders coupling.‘ A quantitative study 
of the susceptibility of this compound, based on 
this assumption, is presented in the present 
paper. 
GENERAL PROCEDURE 
The principal molar susceptibilities may be 
calculated from the expression® 


Ly nL(W,,?/kT) — 2W,,2 Je—W v0 /kT 


a eT 
ne Wake 


(1) 
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Here L is Avogadro’s number. The other con- 
stants are the coefficients of an expansion of the 
energy W, in powers of the applied magnetic 
field 3 


Wr=Wo+WMR+W PIE: ++. (2) 


The summation is over all states of the system. 
In the so-called method of crystalline potentials® 


1L. C. Jackson, Proc. Roy. Soc. A140, 695 (1933). 

2See J. H. Van Vleck, Phys. Rev. 41, 208 (1932). For 
experimental data on the magnetic isotropy of ferric and 
manganous salts see Jackson, reference 1; Krishnan, 
Chakravorty and Banerjee, Phil. Trans. Roy. Soc. 232, 
99 (1933); Krishnan and Banerjee, Nature 135, 873 (1935). 

3 We shall refer to Van Vleck’s first and second papers 
in this issue as V.V. I and V.V. II. 

4 This explanation was first advanced by Van Vleck and 
W. G. Penney, Phil. Mag. 17, 961 (1934), footnote 6. 
See Linus Pauling, J. Am. Chem. Soc. 53, 1367 (1931); 
54, 988 (1932) for an explanation of the low susceptibility 
based on the concept of localized orbitals. 

5See Van Vleck, The Theory of Electric and Magnetic 
Susceptibilities, Ch. VII. 

° H. Bethe, Ann. d. Physik 3, 133 (1929); J. H. Van Vleck, 
reference 5, Ch. XI; W. G. Penney and R. Schlapp, Phys. 


this system is taken to be the magnetic ion sub- 
ject to the constant ionic field of its neighbors in 
the crystal. The problem of finding the sus- 
ceptibility is thus reduced simply to that of 
determining the relative energy levels W,,° of the 
magnetic ion in a crystalline field of a certain 
symmetry and finally the energies in the form 
(2) when a magnetic field is superposed. This is 
accomplished by means of a perturbation calcu- 
lation in which the unperturbed system has the 
type of coupling most appropriate to the size 
of the crystalline field. 

For the time being we shall use the language of 
the method of crystalline potentials and regard 
the interatomic forces of the Fe(CN),.~ ~~ com- 
plex as electrostatic. But we shall later see at the 
conclusion of the paper, that our calculations 
are independent of the nature of the bonding 
mechanism assumed. 

In the present problem we assume that the 
coupling of the individual /’s to a resultant 
orbital angular momentum L is destroyed by a 
large cubic field, which is the dominant portion 
of the crystalline potential. In this respect our 
calculations differ from those of Schlapp and 
Penney,*® on Cr, Co and Ni compounds, wherein 
only the spin-orbit interaction is broken down. 
The unperturbed state of the Fe*** ion in our 
problem is consequently taken to consist of five 
non-interacting 3d electrons. Under the action of 
the cubic potential the single d electron level 
splits into a triply and a doubly degenerate level, 
which we shall designate de and dy respectively.’ 
Gorter® has shown that in an octahedral complex 
such as Fe(CN),~~~ the triplet will lie lower. We 
now make the following approximation. We 
assume that the splitting of the d level is so large 
that we may neglect contributions to the sus- 
ceptibility of perturbation matrix elements 
connecting the. triplet and singlet levels. In 
order to uncouple the individual /’s, this splitting 
must, of course, exceed 3104 cm, the ap- 
proximate separation® of different multiplets 
arising from the d*° configuration of Fet**. 


Rev. 41, 194 (1932); 42, 666 (1932); O. M. Jordahl, ibid. 
45, 87 (1934). 

7For the full significance of the terms de and dy see 
Eqs. (5) and (6) of V.V. I. 

8 C. J. Gorter, Phys. Rev. 42, 437 (1932). 

9 Cf. footnote 12 in V.V. II. 
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Introducing the electron spin we obtain from 
the triply degenerate de level six single-electron 
states, among which the five electrons are to be 
distributed. The energy matrix for this de’ con- 
figuration is the same as that for a single electron 
in a de state except for a change in the signs of the 
proportionality factors in the perturbation. The 
reason for this is that five electrons represent a 
shortage of one from a complete and hence non- 
magnetic group of six electrons in the de® con- 
figuration. Consequently, by virtue of our dis- 
carding the effect of the dy level, we may carry 
through the calculation of the energy levels as 
we would for a one-electron system except for 
this change in signs. The perturbation may then 
be written 


H=[(Ax?+By?—(A+B)2] 
+A’l-s+35B(m,+2m,), (3) 


where A and B are the coefficients of the quad- 
ratic term of a rhombic crystalline potential, 
which is introduced in order to account for the 
observed magnetic anisotropy; A’ is the coeffi- 
cient of the spin-orbit coupling; 3C is the applied 
magnetic field; and 8 is the value of the Bohr 


magneton eh/4xmc. The pattern of energy levels 
produced by the perturbation (3) for the dé 
configuration is shown in Fig. 1. 

Apparently we have specialized things some- 
what in assuming that the rhombic potential is a 
quadratic function of x, y, z. But this is not 
really true, as the energy of our system in a 
rhombic field has three eigenvalues, of which one 
serves as a trivial additive constant since we are 
interested only in the relative energies. Hence 
the full generality of the secular problem arising 
from the lifting of the cubic degeneracy by a 
rhombic field is realized by (3), since (3) in- 
volves two independent constants, A and B, 
which determine the relative eigenvalues com- 
pletely. 


EXPRESSIONS FOR THE PRINCIPAL 
SUSCEPTIBILITIES 


To calculate the energy values of the system 
one must diagonalize the secular matrix cor- 
responding to the perturbation (3). The matrix 
elements are given by Schlapp and Penney and 
Jordahl.* The matrix for the de configuration is 


Cusic 
FIELO 


SPIN-ORBIT 


« , RHOMBIC 
INTERACTION 


 MAONETIC 
FIELO 


FIELO 


/ 
/ 
ee 
de 

\ 

\ 

\ 

\ 
\ | 


on ~ 





Fic. 1. Pattern of energy levels for the de® configuration. 


the same as the sextic corresponding to the 
irreducible representation I; given on p. 90 of 
Jordahl’s paper.’® In the absence of a magnetic 
field this sextic factors into two identical cubics. 

The similar behavior of two of the principal 
susceptibilities (see Fig. 2) suggests that the 
rhombic field, which is responsible for magnetic 
anisotropy, is very nearly axial in symmetry. 
For convenience the rhombic field in (3) may be 


written 
(E+A)x*+ (E—A)y*—2E2’, (4) 


i.e., we set A=E+A and B=E—A. We assume 
that A is small compared with E£. In the limit 
when A=O and the crystalline field has axial 
symmetry it is possible to solve rigorously the 
cubic secular equation involving the elements of 
the axial field and spin-orbit coupling. The com- 
plete matrix corresponding to (3) corrected by 
(4) is then diagonal in all but the elements con- 
taining A, the deviation from axial symmetry, 
and the magnetic field 3. A transformation is 
next made which leaves in off-diagonal positions 
only terms containing square or higher powers of 
A, which are negligible. Lastly, the effect of the 
magnetic field on the energy values is computed 
in the form (2). The final expressions for the 
relative energy levels are given below. The first 
set of three energy levels is for the magnetic 
field parallel to the a and b crystallographic 
axes.!! The upper sign corresponds to the a axis, 


10 Jordahl’s matrix as a whole is inapplicable to our prob- 
lem. His system is a d® configuration which is one electron 
short of a completed configuration d'°. Our system is a d® 
configuration and only because we have chosen to neglect 
the effect of the dy level, does it become one electron short 
of the closed configuration de®, 

11 Although K;Fe(CN). has monoclinic symmetry the 
deviation the orthorhombic symmetry is very slight 
(8 =90° 6’). Consequently the crystallographic axes may be 
considered to coincide with the Cartesian axes along which 
the magnetic field is taken in the calculation. The suscepti- 
bilities which we calculate are for the magnetic field in the 
direction of these axes, whereas the experimental data with 
which the susceptibilities are compared are given for the 
magnetic field along the principal magnetic axes of the 
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the lower sign to the } axis. The second set of energy levels is for the field along the c axis. 
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The meaning of the symbols is as follows: 


K=X8; 6=3\—30—-u; 
u=—2Er/7; m=u—v—X; 
w= —2Ar/7; n=u+v—n; 
v= —A’/2; p=u—v+tr; 
A=(w?+2uv+v")!; g=utvt+r; 
§6=3\+30+4u; t= p/m. 


The coefficients A and B of the rhombic field are 
given by the expressions u-+-w=—2Ar?/7 and 
u—w=—2Br/7. 

When these expressions are incorporated in 
(1) we obtain the principal magnetic suscepti- 
bilities. It is found convenient to express the 
susceptibilities in terms of the square of the em- 
pirical effective Bohr magneton number yz. This 
is defined by the relation x,=(L6?/3kT) (uz?) q, 
where g=a, 6b or c, corresponding to the 
crystallographic axis along which the magnetic 
field is directed. 


COMPARISON WITH EXPERIMENT 


The formulas for the principal susceptibilities 
involve three parameters, A’, A and B, which 
can be evaluated from the observed values of the 
three susceptibilities at a given temperature. 
Actually values of the parameters are chosen to 
reproduce as nearly as possible both the correct 
magnitudes and correct temperature dependence 


crystal. Jackson’s measurements show, however, that the 
latter system of axes, in the temperature range 75° to 
290°K, is never out of coincidence with the crystallo- 
graphic axes by more than a rotation by 3° about the b 
axis. For our purposes this small rotation may be ignored 
and the magnetic and crystallographic axes are taken to 
coincide. 





of the susceptibilities (see Fig. 2). We thus ob- 
tain A’=—510 cm-, Ar?=108 cm- and Br 
=179 cm. The value obtained for A’, the 
coefficient for spin-orbit coupling, is in excellent 
agreement with the value 500 cm~ estimated by 
Van Vleck and Penney‘ from spectral data. No 
appreciable difference would result if the value 
500 cm were used rather than 510 cm™. It is 
interesting that the spin-orbit coefficient can 
thus be indirectly calculated, say to within ten 
percent perhaps, by magnetic calculations. 

In Fig. 2 it may be seen that the magnitude, 
anisotropy and temperature dependence of the 
magnetic susceptibilities are fairly well dupli- 
cated by the theoretical results. The slight 
divergence with increasing temperature which 
Jackson observed for the two similar. suscepti- 
bilities along the a and b axes may well be due to 
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variations in the force field in the crystal caused 
by the increasing energy of vibration. With a 
constant crystalline potential one would expect 
any magnetic anisotropy to diminish with 
temperature, if it changes at all, as the differ- 
ence in the coefficients of the rhombic field be- 
comes less important with respect to kT. The 
kink in the experimental curve for the suscepti- 
bility along the ¢ axis probably also arises from 
a change in the crystalline field. 

Some of the discrepancy with experiment may 
be due to the fact that we have neglected the 
matrix elements connecting dy and de which arise 
from the electrostatic energy ) e?/r:; coupling 
the five electrons together. Hence our calcula- 
tions are not as accurate as those of Schlapp and 
Penney (except for their work on cobalt; cf. 
portion of V.V. II in fine print).” 


GENERALITY OF OuR CALCULATIONS— 
CONCLUSION 


We have so far used the language of the method 
of crystalline potentials for the sake of conven- 
ience in presenting the procedure of calculation. 
However, our calculations still apply with either 
the method of molecular orbitals or with the use 
of the Slater-Pauling localized bonds. 

We have spoken of the Fe*** ion subject to an 
ionic field whereas the ion is probably much less 
polar in K;Fe(CN).>. Consequently, an addi- 
tional perturbation not explicitly included in (3) 
is the spilling over of some of the CN~ electrons 
onto the Fe+++ ion, which represents a tendency 
toward covalent bond formation. A convenient 


12 In Schlapp and Penney’s article, reference 6, the inter- 
electronic energy Le?/r;;, which we shall call (a), is diagon- 
alized in their unperturbed problem, but at the expense of 
considering the complete effect of (b) the cubic potential. 
They retained in Co only the part of (b) diagonal in L. 
We diagonalize (b) rather then (a) in the unperturbed 
system, since by hypothesis (b) is larger than (a) in the 
cyanides (unlike the sulphates). (Even at that we retain 
only matrix elements of (b) internal to the configuration d'. 
Otherwise, the problem involves an infinite number of 
states and is of hopeless complexity.) Fortunately, except in 
Co, there are no matrix elements of (b) nondiagonal in L 
in the states of high spin in which Schlapp and Penney were 
interested (cf. fine print in V.V. II.); otherwise their 
accuracy would be no higher than ours. 
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way of treating covalent bonds is provided by 
the method of molecular orbitals, and is dis- 
cussed more fully in V.V. II. In this method 
there are two types of dy levels, corresponding 
to the upper and lower roots of a secular equa- 
tion which determines the proper linear combi- 
nations of the orbitals of the iron atom with the 
orbitals of the CN~- groups. (Cf., for instance, 
Fig. 2 in V.V. II.) The lower roots are bonding 
and, as they are each inhabited twice, do not 
contribute to the paramagnetic susceptibility. 
The upper roots are antibonding, and are so 
high that we consider them to be virtually 
uninhabited and neglect their effect. We are 
thus concerned only with the de orbitals. 
Fortunately they are nonbonding and do not 
combine with CN- orbitals. Thus the covalent 
“spilling-over”’ effects do not vitiate our calcu- 
lations of the magnetic susceptibility. This result 
is also true if we use Pauling’s localized bonds 
instead of the method of molecular orbitals, for 
his sp'dy? directed bonds represent completed 
electron pairs which do not contribute to the 
susceptibility. Thus in his theory the magnetism 
results from the dé orbitals, which are nonbond- 
ing and so amenable to our type of calculation. 

Pauling assumed‘ that the magnetism of the 
de group arises entirely from the spin. Our calcu- 
lation has shown that actually the orbital angular 
momentum of the de group contributes appre- 
ciably to the susceptibility if the field is nearly 
cubic. This explains why the observed suscepti- 
bility is almost twice as high as that corre- 
sponding to one free spin." 

In conclusion I would like to thank Professor 
J. H. Van Vleck for suggesting this problem and 
for his invaluable advice during the course of the 
work, 


13 The spin-only value of uz is 1.73 for one free spin. The 
reason that Pauling did not seem to need orbital contribu- 
tions to the susceptibility is that he apparently used the 
value ug=2.00 quoted by Bose, Zeits. f. Physik 65, 677 
(1930). This value is out of line with the higher values 
2.3-2.4 reported by Jackson! and by Ishiwara, Sci. Rep. 
Tohoku Univ. 3, 303 (1914); W. Biltz, Zeits. f. anorg. 
Chemie 170, 161 (1928); and L. A. Welo, Phil. Mag. 6, 
481 (1928). 
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Symmetry Considerations Concerning the Splitting of Vibration-Rotation Levels in 
Polyatomic Molecules 


E. BRIGHT WILSON, JR.,* Mallinckrodt Chemical Laboratory, Harvard University 
(Received September 30, 1935) 


The interaction of rotation and vibration and perhaps other effects may split degenerate 
vibration-rotation energy levels of symmetrical polyatomic molecules into a number of com- 
ponents. The permutation symmetry of molecules containing several identical atoms provides 
certain restrictions on this splitting. This paper discusses the maximum number of fine- 
structure components, their quantum weights when nuclear spins are taken into account, and 
the selection rule for transitions. All arguments are based solely on symmetry considerations 
so that no estimate of the magnitude of the splitting is given. 





DEFINITION OF COORDINATES 


N discussing the motions of the atoms in a 

polyatomic molecule (as distinct from the 
motion of the electrons) it is useful to use the 
Eulerian angles 6, y, x to describe the orientation 
of the molecule in space and the normal coordi- 
nates Q;, Qe, Q3, «++, Qw to describe the mutual 
positions of the atoms in the molecule. It is 
important to define these coordinates more 
definitely. Let the angles 0, y, x define a set of 
rotating Cartesian axes X, Y, Z, with origin at 
the center of gravity of the molecules. The 
position of any atom 7 can be specified by the 
Eulerian angles and the coordinates x, y, 2 
relative to the rotating system. However, there 
are three too many coordinates so that three 
relations between the coordinates must be 
written down in order to define the rotating sys- 
tem uniquely. These relations are conveniently 
chosen to be 


=m;(aiyi—bix;) =0, 
Ym ;(bi2i—Ciyi) =0, (1) 
©m;(ci.xi— ;2;) =(, 


where a;, ;, c; are the coordinates (in terms of 
the rotating system) of the equilibrium position 
of the ith atom and m; is the mass of the atom. 
A transformation from the coordinates x;, yi, 2; 
to normal coordinates Q:, Qe, ---, Qw can then 
be carried out in the usual manner and it will 
be found that the relations (1) are equivalent to 
those usually made in normal coordinate treat- 
ments; namely, that to the first approximation 


* Junior Fellow of the Society of Fellows. 
la C, Eckart, Phys. Rev. 47, 552 (1935). 


there is no rotational angular momentum with 
respect to the rotating system of coordinates. 


AN APPROXIMATE WAVE EQUATION 


It is natural to set up the wave equation in 
terms of these coordinates. This involves first 
expressing the kinetic energy in terms of 8, ¢, x, 
Q:, Qe, -:*, Qw and then using the result to 
construct the wave equation in the usual man- 
ner.! This is a possible, but complicated, pro- 
cedure. However, if we assume that we can 
neglect in the classical Lagrangian expression for 
the kinetic enefgy those terms involving the 
normal coordinates to the first or higher powers 
(which means classically that we are assuming 
small vibrations), then the wave equation re- 
sulting from this approximate Lagrangian will 
be simply the sum of two parts; one the equation 
for the rigid rotator (coordinates 6, y, x), the 
other the equation for harmonic vibrations (co- 
ordinates Qi, Qe, ---, Qy).* 

The solutions of this equation are products of 
the type 


y=R(8, Y) x) V(Qi, ict Qvy), (2) 


where R(@, ¢, x) is the rotational wave function 
and V(Q:, ---, Qw) is the vibrational wave 
function. 

For molecules whose equilibrium configura- 
tions possess a certain minimum symmetry; 
namely, at least a threefold axis of. symmetry, 


1B. Podolsky, Phys. Rev. 32, 812 (1928). 

* We assume here that the potential function is harmonic 
but this assumption does not change the final conclusions. 
This method of justifying the usual separation into rota- 
tional and vibrational parts does, however, depend on the 
choice of the rotating coordinate system defined by (1). 
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both the rotational and vibrational energy levels 
may be degenerate when the approximation 
considered here is used. Thus for CH:2F the 
functions R are the symmetrical top solutions 
with the energy levels 


Wa, x=(h?/82*) { J(J+1)/ 


A+K*(1/C—1/A)}, (3) 


so that levels with | K| >0 are doubly degenerate- 
For this molecule certain of the fundamental 
vibrational levels are also doubly degenerate so 
that, to this approximation, CH;F possesses 
some fourfold degenerate energy levels because 
of its threefold symmetry axis. 

It will be shown that this high degeneracy is 
not really required by the symmetry so that 
these energy levels may be split by terms which 
have been omitted in obtaining this simple 
approximation. Some degeneracy, however, will 
remain on account of the symmetry. 





SYMMETRY CONSIDERATIONS? 


The true energy operator, H, is unchanged by 
any permutation of identical particles; it is 
therefore invariant under the group of permuta- 
tions which are equivalent to rotations of the 
molecule. The wave functions RV are not in 
general invariant under these permutations since 
the coordinates @, y, x, Qi, -**, Qy, are not. 
Thus the permutation (123) of the three H atoms 
in CH;F changes x into x+2z/3 and therefore 
R(x) into e?***/8R(x). 

The functions RV belonging to a given ap- 
proximate energy level form a representation’ of 
the group of permutations equivalent to rota- 
tions. This representation is in general reducible 
so that by forming the correct linear combina- 
tions of these functions RV the representation 
can be reduced and each combination will corre- 
spond to one of the irreducible representations of 
the group. The technique of finding the number 


? This paper is an extension of a previous publication 
(J. Chem. Phys. 3, 276 (1935)), which considered the same 
symmetry questions but which obtained the total statistical 
weights of the unsplit energy levels. Here exactly the same 
methods are used, only the split levels will be considered 
and their separate statistical weights determined. 

3 For the meaning and use of the group theory terms em- 
ployed here see E. Wigner, Gruppentheorie (Vieweg and 
Sohn); B. L. van der Waerden, Die gruppentheoretische 
Methode in der Quantenmechanik (Springer, Berlin). 






VIBRATION-ROTATION LEVELS 


of wave functions with each type of symmetry 
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IN POLYATOMIC 


has been treated in detail in a previous paper? 
and will be assumed here. 

If wa, wz, etc., are the linear combinations of 
the functions RV which reduce the representa- 
tion then the perturbed energy levels are deter- 
mined by integrals of the type fYsHwadr. It is 
a consequence of group theory that if y4 and Wz 
belong to different irreducible representations 
the integral is identically zero. Therefore the 
energy matrix, if set up in terms of Wa, Wz, etc., 
will factor into diagonal blocks, the elements of 
each block involving only functions of one 
symmetry type or irreducible representation. 
All other elements of the matrix will vanish. 
Furthermore, blocks corresponding to doubly 
degenerate representations of the group will 
occur twice, so that the energy levels which are 
the roots of such a block of the secular equation 
will be doubly degenerate. Similar remarks apply 
to triply degenerate representations. 





CH;F as AN EXAMPLE 


These statements may be better understood 
if illustrated by an example, CH;F. From 
reference 2, Table X, it is found that if the 
rotational quantum number K has a value not 
divisible by three the rotational functions R for 
that level have the symmetry E with respect to 
the group of permutations equivalent to rota- 
tions. If the vibrational energy level under 
consideration is also doubly degenerate with the 
symmetry £, there will be four wave functions 
of the type RV, having the symmetry EXE 
=2A+E; that is, there can be formed four 
independent linear combinations of the four 
functions RV, such that two of them are sym- 
metrical with respect to permutations equivalent 
to rotations whereas the other two form a 
degenerate pair. Applying the theorems of the 
previous paragraph we obtain the energy matrix 
(or secular equation) shown in Fig. 1, in which 
only the elements in the heavy squares are 
different from zero and the elements of the two 
squares marked E are equal. As far as the 
symmetry is concerned, therefore, the originally 
fourfold level may split, under the influence of 
coupling terms, etc., into three levels, two of 
which are the roots of the quadratic equation 
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Fic. 1. Fic. 2. 


Fic. 1. Energy matrix for a sample approximate level of 
CH;F. Heavy lines enclose non-zero elements. 
Fic. 2. Energy matrix for a sample approximate level 
of shee 


obtained from the upper block of Fig. 1, while 
the remaining level is doubly degenerate and 
comes from either of the two equal lower blocks 
of the secular equation. 


METHANE AS AN EXAMPLE 


Methane may be similarly treated. As a 
sample let us consider the level with J equal to 2 
and with one quantum of one of the triply 
degenerate vibrations excited. From Table III, 
reference 2 we find that the symmetry of this 
level is (E+T)XT=A+E+4T. The secular 
equation for this level therefore has the form 
shown in Fig. 2, in which the lower three blocks 
are identical and the two blocks above these are 
identical. This energy level may therefore split 
into one single level (A), one double level (£), 
and four triple levels (7). 

In connection with methane it may be pointed 
out that it is possible, but not proved, that even 
the nonvibrating rotational levels of methane 
may be split. For example, the level with J 
equal to 5 has the symmetry E+3T so that in 
the nonvibrating state it is possible for a per- 
turbation with tetrahedral symmetry (such as 
conceivably the stretching terms) to split this 
eleven-fold level into one double level and three 
triple levels. It requires further study of a more 
quantitative nature in order to decide whether 
or not such a splitting actually occurs to an 
appreciable degree or not. This treatment merely 
shows that the symmetry of the molecule does 
not require such a high degree of degeneracy. 


JR. 


SELECTION RULES 


The electric moment of a molecule is a function 
of the coordinates of the nuclei and is unaffected 
by interchange of coordinates of identical nuclei. 
It is thus symmetrical. The integral J Yauyadr, 
determining the transition probability, is there- 
fore zero if y4 and yz are in different irreducible 
representations. The selection rule is that tran- 
sitions occur only between levels of the same 
symmetry. This rule does not replace the ordi- 
nary selection rules but is an additional re- 
striction. 


STATISTICAL WEIGHTS OF COMPONENT LEVELS? 


The calculation of the statistical weights of 
the components of a level (important for the 
computation of relative intensities) is carried out 
in exactly the same manner as in reference 2. 
Thus for methane the nuclear-spin functions 
have the symmetry 5A+E£+3T. The total wave 
function must have the symmetry A. Therefore 
the components of the sample vibration-rotation 
level previously discussed will have the weights 
shown in Fig. 3 (multiplied by 2/+1). 


, 
MOLECULES WITH SEVERAL FRAMEWORKS 


Both of the examples above .are molecules 
possessing two ‘‘frameworks” in the sense used 
in reference 2. As discussed in that reference, 
however, this fact does not need to be taken 
into account in calculating the statistical weights 
of the levels since the total wave function must 
have one definite symmetry with respect to 
interchanges of identical particles and there will 


A {5A+E+3T } 5 
E :2A+6E+6T; 2 
T i 3A+3E+BT! 3 
7 
T 
r 








A+E+4T} 





A is ¥e 














3 
: it 3 
I nm: m iw 

Fic. 3. Number of components, their symmetries and 
statistical weights (factor 2/+1 omitted) for a sample 
level of methane. Regions I and II give the symmetries 
before and after splitting, spins not included. Region III 
gives the symmetries when thespin function ys(5A +E£+3T) 
is included. Region IV gives the statistical weights of the 
components when the exclusion principle is applied (A 
only). The vertical order and relative spacing of the com- 
ponents has no significance. 





STRUCTURES OF SO:z, 


be only one linear combination of the functions 
based on the various frameworks which will have 
the correct symmetry. 


CONCLUSION 


It is to be emphasized that only symmetry 
considerations have been used in this paper so 
that the results are necessarily incomplete since 
they do not yield any information regarding the 
magnitudes of the splittings to be expected. 
Actual energy levels may not be split to an 
observable extent or they may show the splitting 
indicated here only incompletely. It is possible 
also that certain of the perturbations, especially 
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the rotation-vibration coupling, may split the 
approximate energy levels to such an extent that 
a component of one approximate level may 
nearly coincide with a component of another 
approximate level. In spite of these complica- 
tions, however, the symmetry does restrict the 
maximum number of components which are 
possible and gives the selection rules and sta- 
tistical weights involved. Experimental results‘ 
exist which may be capable of interpretation in 
terms of these effects but further work is nec- 


essary. 


4W. B. Steward and H. H. Nielsen, Phys. Rev. 47, 828 
(1935). 





DECEMBER, 1935 


JOURNAL OF CHEMICAL 


PHYSICS VOLUME 3 


The Molecular Structures of Sulfur Dioxide, Carbon Disulfide, and Carbonyl Sulfide 


Pau C. Cross* anp L. O. Brockway, Gates Chemical Laboratory, California Institute of Technology 
(Received October 7, 1935) 


The structures of the molecules SO2, CS2, and COS have 
been investigated by the electron diffraction method with 
the following results; SO,, S-O=1.46+0.02A; CS, 
C-—S=1.5440.03A; COS, C—O=1.16+0.02A, C-—S 
= 1.56+0.03A. The types of bond arrangement compatible 
with these interatomic distances are discussed. In SO, 
the molecule resonates between the structures having 


HE molecular structures of sulfur dioxide, 

carbon disulfide, and carbonyl] sulfide have 
been studied by the method of electron diffrac- 
tion of gases. A jet of vapor is made to intersect 
a beam of forty-kilovolt electrons whose diffrac- 
tion pattern is registered-on x-ray film. Visual 
measurements on the maxima and minima are 
correlated with those on the theoretical curves 
based on various molecular models. In these 
curves the scattering power of an atom for 
electrons is assumed to be constant with angle 
and the inelastic scattering is neglected. A 
detailed discussion of the procedure used in 
interpreting electron diffraction photographs 
has already been published.’ 


* National Research Fellow in Chemistry. 


t No. 501. 
1L. Pauling and L. O. Brockway, J. Chem. Phys. 2, 867 


(1934). 


single-double and double-single bonds between the sulfur 
and the two oxygen atoms, with a bond angle of 122°+.5°. 
CS, is a linear molecule with the structure having the 
two double bonds predominating over those having a single 
and a triple bond. In COS the double-double bond arrange- 
ment and the structure having the triple carbon-oxygen 
bond predominate. 


SULFUR DIOXIDE 


The sample of sulfur dioxide was redistilled 
from the purest commercial product. 

The heavier photographs showed four maxima. 
The first was very strong; the second, strong 
and rather broad; the third and fourth were 
weaker but very sharp. The theoretical intensity 
curve, Fig. 1, is based upon a triangular model 
with an angle of 120° between the sulfur-oxygen 
bonds. No other curves are shown since the 
theoretical diffraction pattern is very insensitive 
to changes in bond angle. This effect in molecules 
with heavy central atoms and light outer atoms 
has been discussed for chlorine dioxide.” 

The values obtained for the sulfur-oxygen 
distance are given in Table I. The third column 


*L. O. Brockway, Proc. Nat. Acad. Sci. 19, 303, 874 
(1933). 
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Fic. 1. Theoretical intensity curves for sulfur dioxide and 
carbon disulfide. 


shows the observed positions of the maxima 
and minima, where @ is the scattering angle and 
dis the wave-length associated with the electrons. 
The fourth column gives the positions of the 
maxima and minima on the theoretical curve 
and in the last column are listed the values of 
the sulfur-oxygen distance. The final value is 
S—O=1.46+0.02A. 

An earlier electron diffraction investigation of 
sulfur dioxide® leads to the value 1.37A. This was 
obtained with the aid of photographs showing 
only one maximum and is accordingly less 
reliable than the present result based upon 
measurements of four maxima and two minima. 

The bond angle in SO, has been evaluated from 
the studies of the vibrational spectrum in the 
infrared region and by the Raman technique.‘ 
With the aid of the fundamental frequencies 
(ve =1361, v-=1152, vs=525 cm~')® correspond- 
ing to the three modes of vibration, the treatment 
of the molecule, with the two-constant ‘‘valence 


3R. Wierl, Ann. d. Physik 8, 544 (1931). 

4C. R. Bailey and A. B. D. Cassie, Proc. Roy. Soc. A137, 
605 (1933); A138, 531 (1932). See also A. Jonescu, Comptes 
rendus 196, 1476 (1933); 197, 35 (1933); and J. H. Clements, 
Phys. Rev. 47, 224 (1935). We have not considered Jon- 
escu’s structure because of the incomplete nature of the 
rotational analysis upon which it is based, and the de- 
pendence upon Wierl’s erroneous value of theS— Odistance. 

& 3) M. Badger and L. G. Bonner, Phys. Rev. 43, 305 
(1933). 
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TABLE I. Sulfur dioxide. 











FOUR PHOTOGRAPHS ; A =0.0606 
Max. Min. ssn of? x 
1 5.625 7.86 (i 
2 9.656 14.30 
3 11.99 17.27 
F 13.76 20.° 1. 
4 16.19 «ae 
4 18.15 26. 9 ‘ 
Average 


S—O=1.46+0.02A 








force” potential function and the assun: 
that the vibration is harmonic, leads to a 
of 122°+5°. If this angle and the vibrat... 
frequencies are used in conjunction with 
Badger relation between the force constant 
interatomic distance,® the sulfur-oxygen separ 
tion is calculated to be 1.48A, in satisfacte 
agreement with our results. 

The electronic structures which seem likely J 
contribute to the ground state of sulfur dioxid 
are as follows: 


A structure can be written having single bonds 
and only three electron pairs on one of the oxygen 
atoms; but this is less likely than the third of the 
structures above because of the greater electron 
affinity of oxygen. In view of the theoretical 
possibility of placing five electron pairs on the 
sulfur atom, structures may be considered 
having two double bonds or a single and a triple 
bond, but these undoubtedly make very small 
contributions. 

The observed bond distance is 0.05A less 
than the sum of the covalent radii’? for double 
bonds. The double bond distance required by 
resonance between I and II above would be 
shortened by the polar character of the bonds. 
The inclusion of III in the resonating structures 
would not affect the distances but this structure 
is probably involved because of the tendency of 
sulfur to assume the double plus charge as in the 
sulfates and thionates. 


®*R. M. Badger, J. Chem. Phys. 2, 128 (1934). 
7L. Pauling, Proc. Nat. Acad. Sci. 18, 293 (1932). 
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Fic. 2. Theoretical intensity curves for carbonyl] sulfide. 


The angle is determined by the separate elec- 
tronic structures and the value determined 
spectroscopically is just that to be expected 
between a single and a double bond on a tetra- 
hedral atom. 


TABLE II. Carbon disulfide. 











EIGHT PHOTOGRAPHS © A =0.0607A 

Max. Min. frst x a 
4.713 7.25 1.538 
la 6.312 9.85 1.560 
Z 7.623 11.52 1.512 
2 8.698 13.55 1.559 
2a 10.63 16.10 1.515 
3 11.63 17.84 1.535 
3 12.65 19.84 1.568 
3a 14.58 22.40 1.536 
4 15.54 24.12 1.553 
4 16.81 26.12 1.554 


Average 1.543 
C—S=1.54+0.03A 











CARBON DISULFIDE 


The carbon disulfide photographs were char- 
acterized by four strong maxima each followed 
by a weak maximum and a deep minimum. 
This appearance, which is also observed in the 
theoretical curve corresponding to a linear 
model (Fig. 1), is due to the fact that one of the 
two different interatomic distances in the mole- 
cule is just twice the other. The results given in 
Table II lead to a carbon-sulfur distance of 
1.54+0.03A in essential agreement with Wierl’s 
value’ of 1.58A. 

The electronic structures to be considered are 
the three having the double-double, the single- 
triple and the triple-single arrangement of 
bonds. If the three were of nearly equal impor- 
tance the triple bond distance should be observed 
as in the analogous case of carbon dioxide. The 
table of radii gives 1.63A for the double bond and 
1.43A for the triple bond. Observations on the 
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TABLE III. Carbonyl sulfide. 








» =0.0602A 


@1,35 


FIVE PHOTOGRAPHS 


‘ 4x sin 0/2 
Max. Min. Ad 1.35 1.30 





1.156 
1.171 
1.147 
1.144 
1.135 
1.153 
1.172 


1.151 


6.18 
140 
8.90 
10.10 
11.45 
12.90 
14.60 
15.90 
16.60 
18.20 1.162 
20.10 Lisi 
Average 1.157 
(Ave. Dev. 0.011 
1.16+0.02A 
1.56+0.03A 


6.05 
7.51 
8.80 
9.95 
11.15 
12.60 
14.30 


16.15 
17.80 
19.60 


1 5.23 
6.41 
7.67 
8.70 
9.82 

10.93 
12.20 
13.33 
14.04 
15.31 
16.65 


C-O= 
C-S= 








spectroscopic molecules P2 and PN§® indicate that 
the ratio of triple to single bond radius for ele- 
ments in the second row of the periodic table 
should be larger than the value 0.79 used in the 
radius table. From 79>=1.86A in P=P the triple 
bond radius for phosphorus is 0.93A. Adding to 
this the triple bond radius for nitrogen (0.55A) 
we get 1.48A in agreement with the observed value 
for P=N of 1.49A. Accordingly with the value 
0.93A for the triple bond radius of phosphorus 
the triple-single bond radius ratio becomes 0.846 
for second row elements. Applying this ratio to 
sulfur and adding the triple bond radius of car- 
bon we obtain 1.49A for the predicted value of 
the carbon-sulfur triple bond. 

The observed value 1.54A is 0.05A larger than 
the triple bond distance; therefore the contribu- 
tion of the second and third structures above is 
considerably less important than that of the 
first. This conclusion is supported by the very 


8G. Herzberg, Ann. d. Physik 15, 677 (1932); J. Curry, 
‘ie and G. Herzberg, Zeits. f. Physik 86, 348 
(1933). 
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small tendency of second row elements to form 
triple bonds; in fact, the two compounds of 
phosphorus mentioned in the preceding para- 
graph are the only known examples, and they 
exist only in an electrical discharge. 


CARBONYL SULFIDE 


Theoretical curves were calculated for four 
linear models of carbonyl sulfide corresponding 
to different ratios (a values in Fig. 2) of the 
carbon-sulfur to carbon-oxygen distances. That 
the second ring was observed to be very nearly 
as strong as the third and the maximum beyond 
the fourth was scarcely more than a shelf fixes 
the ratio of the two distances at about 1.34. 
Table III includes the calculations for the two 
models with ratios of 1.35 and 1.30. The final 
values are C—O=1.16+0.02A and C—S=1.56 
+0.03A. 

Dornte® obtained electron diffraction photo- 
graphs of carbonyl sulfide and reported calcula- 
tions of the theoretical scattering for only one 
model (corresponding to a distance ratio of 1.40). 
This ratio and the distance which he reports 
(C—O=1.13+0.05A, C—S=1.58+0.08A) differ 
from ours by less than the uncertainty which he 
ascribes to his measurements. 

The value observed for the carbon-sulfur 
separation is essentially the same as that in the 
CS. molecule while the carbon-oxygen distance 
agrees with the value 1.15A in carbon dioxide 
and with the triple bond radius sum 1.13A. The 
foregoing discussion of the character of the 
carbon-sulfur bond applies also to the case of 
carbonyl sulfide. 

We are grateful to Professor Linus Pauling 
for suggestions and criticisms. 


*R. Dornte, J. Am. Chem. Soc, 55, 4126 (1933). 
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The thermodynamic properties of sulfur dioxide, carbon disulfide, and carbonyl sulfide are 
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calculated from the molecular constants obtained by electron diffraction and Raman and 
infrared spectra. — (F°— Ep) /T 298.1 = 50.95, 48.28, 47.39; $298.1 = 59.40, 56.84, 55.34 for SO, 
CS.(g), and COS, respectively. The results are applied to several reactions involving these 
compounds. The free energies of formation of CS,(l) and COS at 298.1°K are +15.24 and 


— 40.48 kg.cal. 








HE results of the preceding paper! give the 
following values for terms which appear 

in the approximate equations for the entropy and 

free energy functions such as given by Kassel :? 


SO2, $R1n Lal eI c= —262.93, R In p= —1.377, 
CSz, Rin J=—172.03, Rin p= —1.377, 
COS, Rin J=—173.23, Rin wp=0. 


The vibration frequencies were taken as follows: 


SOz, 1361, 1152, 525,3 
CS2, 1523, 655.5, 396.8 (2), 
COS, 2079, 859, 527 (2). 


Because of the lack of adequate data to determine 
the anharmonic factors, the molecules were 
treated as harmonic oscillators and the vibra- 
tional contributions determined from the tables 
of Gordon and Barnes.* This also ignores the 
Fermi interaction, but the errors introduced are 
not serious for work of the accuracy intended 
herein. 

The results shown in Table I were obtained 
from the equations; 


SO», (F°—E,°)/T= —4R In T-—5.49 
| +((F°—Eo°)/T Jviv., 
S°=4R In T+13.44+ Sy in., 


* National Research Fellow in Chemistry, 1933-35. 
Now at Stanford. 
+ No. 502. 


1 P. C. Cross and L. O. Brockway, J. Chem. Phys. 3, 821 
(1935). 

2L. Kassel, J. Am. Chem. Soc. 55, 1351 (1933). 

3R. M. Badger and L. G. Bonner, Phys. Rev. 43, 305 
(1933). 

4D. M. Dennison and N. Wright, Phys. Rev. 38, 2077 
(1931). 


5C. R. Bailey and A. B. D. Cassie, Proc. Roy. Soc. 
A140, 605 (1933). 

6 A. R. Gordon and C. Barnes, J. Chem. Phys. 1, 297 
(1933). 
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C,°=7.946+ Cviv.; 
(F°— Eo) /T = —(7/2)R In T—7.93 
+[(F°—Eo°)/T lviv., 
5° = (7/2)R In T+14.88+ Svin., 
C,°=6.954+ Cvin.; 
COS, (F°—£E,°)/T= —(7/2)R In T—7.40 
+((F°—Eo°)/T lviv., 
S°=(7/2)R In T+14.354+ Sviv., 
C,°=6.954+ Cvin.. 


CS», 


The values for SO, are probably accurate to 
within 0.15 unit, while those for CS, and COS 
should be within 0.1 unit of the true values. 


SULFUR: DIOXIDE 


A similar calculation for SOz has recently been 
made by Gordon,’ who applied his results to the 
reactions 


1So+O2=SOn, (1) 
SO+402=SOz, (2) 
SO.+3H2*=H.S°+2H20," (3) 
2CO24+4S,2=SO.+2C0." (4) 


(References are to the sources of (F°—E,°)/T 
values.) Since the new values of —(F°—E,°)/T 
are 0.38 larger than Gordon’s, the values of log 
K, and log Kz in Gordon’s Table I should be in- 


7A. R. Gordon, J. Chem. Phys. 3, 336 (1935). 

8 W. F. Giauque, J. Am. Chem. Soc. 52, 4816 (1930). 

9P. C. Cross, J. Chem. Phys. 3, 168 (1935). 

10 A. R. Gordon, J. Chem. Phys. 1, 308 (1933). 

11L. Kassel, J. Am. Chem. Soc. 56, 1838 (1934). 

12C, W. Montgomery and L. Kassel, J. Chem. Phys. 2, 
417 (1934). Corrected by Cross, reference 9. : 

13 J. O. Clayton and W. F. Giauque, J. Am. Chem. Soc. 
55, 5071 (1933). 
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TABLE I. Free energy, entropy, and heat capacity of SOz, CS2, and COS. 
































—(F°—Eo)/T So Co 
T SO2 CS2(g) COs SO2 CS2(g) COs SO2 CS2(g) COS 
298.1 50.95 48.28 47.39 59.40 56.84 55.34 9.51 10.91 9.92 
400 53.49 50.90 49.81 62.32 60.19 58.41 10.35 11.85 10.95 
500 55.50 53.04 51.79 64.72 62.91 60.93 11.08 12.52 11.66 
600 57.21 54.89 53.50 66.79 65.24 63.11 11.67 13.00 12.21 
700 58.71 56.51 55.01 68.62 67.27 65.03 12.11 13.37 12.64 
800 60.05 57.97 56.37 70.26 69.08 66.74 12.44 13.65 12.98 
900 61.27 59.30 57.61 71.74 70.70 68.28 12.70 13.86 13.26 
1000 62.39 60.51 58.75 73.09 72.16 69.69 12.90 14.02 13.49 
1100 63.42 61.63 59.81 74.33 73.50 70.98 13.05 14.15 13.68 
1200 64.37 62.67 60.79 75.47 74.74 72.18 13.17 14.26 13.84 
1300 65.27 63.65 61.71 76.53 75.89 73.30 13.27 14.35 13.97 
1400 66.11 64.57 62.58 tae 76.96 74.34 13.35 14.42 14.08 
1500 66.91 65.44 63.40 78.44 77.97 75.31 13.42 14.48 14.17 
1600 67.66 66.25 64.17 79.30 78.91 76.22 13.47 14.53 14.25 
1700 68.37 67.02 64.90 80.12 79.79 77.09 13.52 14.57 14.32 
1800 69.04 67.75 65.60 80.90 80.61 77.91 13.56 14.61 14.38 
TABLE II. 
Eq. Reactions AE AF°/T =—Rin K 
T =298.1 400 600 900 1500 
(3) 3H.2+SO.2.=H2:S+2H:0 —47.98 —152.2 — 109.7 —67.5 — 38.5 —14.5 
(4) 2CO.+4S2=2CO0+S0, 47.65 139.0 97.2 56.5 29.4 8.1 
(5) CO.+H2S = COS+H.20 7.0 23.7 17.6 11.5 7.4 4.3 
(6) COS+H.2S=CS.+H20 8.4 31.3 23.9 16.7 11.8 8.0 
(17) CO+43S2.=COS —22.2 —57.9 — 38.3 —19.1 —6.4 —3.6 
(18) 2COS=C0O2+CS:2 1.4 7.6 6.4 Ie 4.4 3.8 
(19) CS,=C(graph) +S. 3.6 12.6 9.7 6.9 5.2 3.8 
> 
creased by 0.08. By using the values of AE»° From the reactions 
suggested by Gordon, several values of AF°/T 
“ J ' H.+3S:=H»S, AE=—19.6, (7) 
=—RInK have been calculated for reactions 
(3) and (4) and given in Table II. Comparisons H.+302.=H.0O, AE »°=—57.1,! (8) 
with experimental results are found in Table ITI. 
2S=Sz, AE,°= —_= 102.6," (9) 
CARBON DISULFIDE AND CARBONYL SULFIDE 20 =O», AE, = —117.3," (10) 
Since the thermal measurements on these one obtains 
nds are known to be subject to some 
mcneliggenir E,(H.0) — Ey°(H2S) = —44.9. (11) 


errors, values of AEo° were obtained from equi- 
librium measurements. The results of Terres and 
Weseman" on the reactions 


CO2.+H2S=COS+H:20, AE o®=7.0kg. cal. (5) 
and 
COS+H.S=CS, +H.0, AEo®=8.4 (6) 


gave exceptionally consistent results having no 
pronounced trend with temperature. 


14E,. Terres and H. Weseman, Zeits. f. angew. Chemie 
45, 795 (1932). 


Combining with the results on reactions (5) and 
(6) gives 


E,°(COS) — ECO) =51.9, (12) 
E,(CS2) —E,(COS) =53.3 (13) 
Since EF (CO) —E,(COz) =125.4," (14) 
E.(CO) —£E,(CS2) =20.2, (15) 


1 A. Christy and S. M. Naudé, Phys. Rev. 37, 903 


(1931). 
16 B. Lewis and G. von Elbe, J. Am. Chem. Soc. 57, 612 
(1935). 
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TABLE III. 




















—-Rink Calc.- 
Eq. T Calc. Obs. Ref. Obs. 
(3) 1160 — 25.1 —27.1 18 +2.0 
1362 —18.2 —19.8 18 +1.6 
1473 —15.2 — 16.2 18 +1.0 
1645 —11.2 —11.7 18 +0.5 
(4) 1275 13.7 16.3 19 —2.6 
1460 8.9 10.9 19 —2.0 
(5) 623 2 1123 14 —0.1 
673 10.3 10.3 14 0 
723 9.5 9.5 14 0 
773 8.8 8.8 14 0 
823 8.2 8.2 14 0 
873 ta ia 14 +0.2 
(6) 973 11.1 11.3 14 —0.2 
1023 10.7 11.0 14 —0.3 
1073 10.3 10.3 14 0 
1123 9.9 9.7 14 +0.2 
1173 9.6 9.0 14 +0.6 
(17) 533 23.9 22.5 20 +1.4 
575 20.7 19.1 20 +1.6 
1073 ye FB 21 +0.1 
1223 —0.2 —1.4 21 +1.2 
(18) 533 2.5 8.2 22 —2.7 
575 Ee 9.4 22 —4.1 
573 5.3 5.3 21 0 
673 4.9 4.9 21 0 
753 4.7 aa 21 +1.0 
E,.(CO) —E,(COS) =73.5. (16) 


These results are expressed as differences inas- 
much as there still seems to be some doubt as to 
the absolute values of Eo° for carbon compounds. 

The £,° differences given here have been ap- 
plied to the reactions 


CO+ 3S2= COS, AE, = — 22.2, (17) 
2COS =CO2+CSo, AE,°=1.4, (18) 
CS.=C(graphite)!7+Se, AEo®= 3.6. (19) 

The reaction 
CO =C (graphite) +302, AEo®= 27.4!" (20) 


was also employed to obtain A£,° for reaction 
(19). The results are summarized in Tables II 
and III. 


17 J. O. Clayton and W. F. Giauque, J. Am. Chem. Soc. 
54, 2610 (1932). 








bo 
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The comparisons on reactions’ (17) and (18) 
are satisfactory considering that the values of 
Lewis and Lacey on (18) are admittedly only 
approximate. The large discrepancies between 
observed and calculated results for (19) have 
been discussed by Randall.** The experimental 
measurements are in such poor agreement with 
each other that it seems justified to consider the 
indirectly calculated values as the most reliable 
so far available. The experimental results are 
not given here because this reaction will be con- 
sidered more fully in a forthcoming paper by 
Professor Randall. 

The free energies of formation of CS, and COS 
may be obtained from these results as follows: 


C(graphite)!7+S,"=CSo(g), AFoos.1= —3.75, 


2S(rh)=Se, AF 98.1 = +19.44,? 


CS2(g) =CS2(1), AFoos.1.= —0.45,"4 
C(graphite) + 2S(rh) = CSo(1), AFeos.1= +15.24. 


The value obtained from a rather uncertain 
value of the heat of combustion of CSz is 17.15.*4 


C(graphite) +302” 
+4S. = COS: A Fogs.1 = — 50.20, 
S(rh) = 3S, A Foo3.1 = 9.72, 


C(graphite) +302 
+S(rh) =COS, A Foog.1 = — 40.48. 


This value is in essential agreement with the 
thermal value of —39.80.*4 


18M. Randall and F. R. Bichowsky, J. Am. Chem. Soc. 
40, 368 (1918). 

19 J. B. Ferguson, J. Chem. Soc. 40, 1626, 1900 (1918). 

20 Lewis and Randall, Thermodynamics (1923), p. 581. 
Recalculated with Gordon’s value for the free energy of 
formation of'Ss, reference 7. 

*1 A. Stock, W. Siecke and E. Pohland, Ber. 57, 719 
(1924). 

22 G. N. Lewis and W. N. Lacey, J. Am. Chem. Soc. 37, 
1976 (1915). 

23 Randall, Am. Chem. Soc. San Francisco, August, 1935. 

*4 Parks and Huffman, Free Energies of Some Organic 
Compounds (1932), pp. 204-206. 

25H. L.. Johnston and M. K. Walker, J. Am. Chem. Soc. 
55, 172 (1933). 
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The Molecular Structure of Nickel Carbonyl 


L. O. Brockway AND Paut C. Cross,* Gates Chemical Laboratory,** California Institute of Technology 
(Received October 7, 1935) 


Electron diffraction by the vapor of Ni(CO), leads to a 
molecular model in which the four carbonyl groups have 
a tetrahedral arrangement about the nickel atom with the 
distances Ni—C =1.82+0.03A and C—O=1.15A. These 
distances are compatible with resonance between two 
electronic structures in which the C—O bond resonates 
between triple and double electron pair bonds and the 


Ni—C bond between single and double electron pair bonds. 
Nickel carbonyl is the first quadricovalent compound of 
neutral nickel whose structure has been determined, and 
its tetrahedral configuration is contrasted with the square 
arrangement of bonds in the quadricovalent compounds of 
divalent nickel ion. 





HE molecular structure of nickel carbony] 

and its relation to the more general problem 
of the orientation of the four bonds in quadri- 
covalent nickel compounds has interested chem- 
ists for more than ten years. In 1926 Blanchard 
and Gilliland! gave a detailed discussion of the 
chemical properties of nickel carbonyl and 
deduced an electronic structure which was in- 
tended to be compatible with these properties 
and at the same time would allow for a complete 
shell of eighteen electrons around the nickel 
atom. It is interesting to note that their sug- 
gestion that four electrons are shared between 
each carbon and the nickel atom is equivalent 
in the modern concept of valence to double 
electron pair bonds between the carbon and 
nickel atoms, to which we shall resort in order 
to account for our observations. 

In 1930 Sugden? proposed a cyclic structure 
for nickel carbonyl in which the nickel and the 
four carbon atoms form a five-membered ring 
with an oxygen atom attached to each carbon. 

Pauling*® predicted from a theoretical consider- 
ation of the relative strengths of bonds formed 
by s, p, and d eigenfunctions that in the com- 
pounds of divalent nickel ion in which four 
covalent bonds are attached to the nickel atom 
the four bonds should be directed toward the 
corners of a square. In support of this prediction 
he cited the isomorphism of K2Ni(CN),4:H2O 


* National Research Fellow in Chemistry. 

** Contribution No. 503. 

1A. A. Blanchard and W. L. Gilliland, J. Am. Chem. 
Soc. 48, 872 (1926). 

2S. Sugden, The Parachor and Valency (Albert Knopf, 
New York, 1930), p. 189. 
932). J. Am. Chem. Soc. 53, 1367 (1931); 54, 988 
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with the corresponding palladium compound and 
the known‘ square arrangement of the bonds in 
the chloropalladites and chloroplatinites. Because 
of later misinterpretations it should be pointed 
out that the same theoretical considerations 
suggest a tetrahedral arrangement when four 
bonds are formed with neutral nickel, as in the 
case of nickel carbonyl. 

Sugden® found further evidence in favor of 
the square arrangement in quadricovalent nickel- 
ous compounds in the preparation of intercon- 
vertible isomers of the nickel derivatives of 
methyl benzyl-, methyl u-propyl-, and methyl 
n-butyl glyoxime. These were explained as cis- 
and trans-isomers, which could not exist with a 
tetrahedral arrangement of the four bonds. A 
similar relation was found much earlier by 
Tschugaeff*® for the nickel derivatives of mono- 
ethylglyoxime but he made no attempt to ex- 
plain this isomerism. Additional crystallographic 
evidence is found in the isomorphism’ of 
BaNi(CN),-4H:O and BaPt(CN),-4H.0. Even 
more conclusive is the recent work of Cox, Pink- 
ard, Wardlaw and Webster® on the nickelous deri- 
vative of salicylaldoxime in which two molecules 
of the aldoxime are joined by four bonds to the 
nickel atom. Oscillation and rotation x-ray 
photographs show that the crystal belongs to the 
space group P2,/m—Cz2,*°. With two molecules 
in the unit cell (as determined from the dimen- 

4R. G. Dickinson, J. Am. Chem. Soc. 44, 2404 (1922). 

5S. Sugden, J. Chem. Soc. 1932, 246. H. Covell and S. 
Sugden, ibid. 1935, 621. 

® Tschugaeff, J. Russ. Phys. Chem. Soc. 42, 1466 (1910). 

7H. Brasseur, A. de Rassenfosse and J. Pierard, Zeits. f. 
Krist. 88, 210 (1934).H. Brasseur and A. de Rassenfosse, 
Mem. Soc. Sci. Liége. [3] 18, 1 (1933). 


8E. G. Cox, F. W. Pinkard, W. Wardlaw and K. C. 
Webster, J. Chem. Soc. 1935, 459. 














MOLECULAR STRUCTURE OF NICKEL 


sions of the unit and the density) the molecule 
is required to have a center of symmetry. This 
is possible only if the molecule is planar with a 
trans-configuration of the two salicylaldoxime 
groups. Robertson® has similarly investigated 
the nickel derivative of phthalocyanine and 
found that the molecule has a center of sym- 
metry, which requires a planar configuration of 
the four bonds to the nickel atom. 

For the case of nickel carbonyl the zero elec- 
tric moment observed by Sutton and Bentley’ 
indicates that the molecule must have either the 
square or tetrahedral configuration (with the 
carbonyl groups collinear with the nickel-carbon 
bonds) as opposed to a cyclic arrangement. 
These authors favored the tetrahedral configura- 
tion on the basis of Pauling’s prediction. 

The Raman and ultraviolet absorption spectra 
of nickel carbonyl have been studied by Dadieu 
and Schneider,"' Anderson,” Thompson and 
Garratt,!* and Duncan and Murray." The first 
two investigators pointed out that the almost 
identical frequency of vibration of the carbonyl 
groups in nickel carbonyl and the carbon 
monoxide molecule is strong evidence for an 
identical structure of the carbon-oxygen bond in 
the two cases. Duncan and Murray investigated 
both the Raman and ultraviolet absorption 
spectra as well as the vibrations of a mechanical 
model of a square planar molecule. They ascribe 
incorrectly to Pauling a prediction of the square 
configuration for nickel carbonyl. From the fair 
agreement between their assignment of the fre- 
quencies in the spectra and the frequencies 
observed in the model and from a calculation of 
the symmetry number they conclude that the 
plane square model is correct or at least more 
probable than the tetrahedron. Wilson’ points 
out that the failure of Duncan and Murray to 
observe the selection rules in assigning the 
observed frequencies to modes of vibration of 


9 J. M. Robertson, J. Chem. Soc. 1935, 615. 

10 L. E.Sutton and J. B. Bentley, Nature 130, 314 (1932). 
L. E. Sutton, R. G. A. New and B. Bentley, J. Chem. Soc. 
1933, 652. 

1 Dadieu and Schneider, Anz. Acad. Wiss. Wien, Math. 
Naturwiss. Klasse 68, 191 (1931). 

12 J. S. Anderson, Nature 130, 1002 (1932). 

13H, W. Thompson and A. P. Garratt, J. Chem. Soc. 
1934, 524. 

‘4B. F. Duncan and J. W. Murray, J. Chem. Phys. 2, 
636 (1934). 
1 E. B. Wilson, Jr., J. Chem. Phys. 3, 59 (1935). 
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the molecule has led to an incorrect assignment 
of the frequencies; hence, their conclusions are 
invalidated. 

Laird and Smith'* have recently redetermined 
the parachor value of nickel carbonyl but reached 
no conclusion concerning the configuration of 
the four bonds. 

The present investigation was undertaken in 
the hope of obtaining decisive evidence on the 
structure of nickel carbonyl and determining 
the values of the interatomic distances, which 
would be useful in considering the types of bonds 
involved. 

The sample of nickel carbonyl was obtained 
from the Eastman Kodak Company. This was 
vacuum distilled and pumped off while frozen 
with liquid air in order to remove carbon 
monoxide. The small residue of nickel hydroxide 
produced by hydrolysis was left behind in the 
distillation. The compound was photographed 
with the sample at room temperature with a 
diffraction point to film distance of 10.43 cm 
and with the use of electrons having a wave- 
length of 0.0606A. 

The photographs show eight maxima with a 
characteristic appearance of the diffraction pat- 
tern. The first maximum is weak while the second 
and third are both very strong and very sharp. 
The fourth maximum is weak and lies close to 
the third while the minimum following is broad 
and well defined. The fifth maximum is strong 
and sharp with the sixth, seventh, and eighth 
rings falling off in intensity in regular order. 

Quantitative estimates were made of the 
apparent intensities of the successive rings and 
these data were combined with the observed 
values of the function (47 sin 6/2)/X to obtain 
an atomic distribution curve in the manner 
described by Pauling and Brockway.!’ This 
curve, reproduced in Fig. 1, represents the 
average distribution of scattering power around 
the atoms in the molecule of nickel carbonyl. 
Since the electrons in the beam are scattered 
chiefly by their interaction with the atomic nuclei 
the positions of the maxima in the distribution 
curve represent internuclear distances: in the 


16 F, W. Laird and M. A. Smith, J. Am. Chem. Soc. 57, 
266 (1935). 

17. Pauling and L. O. Brockway, J. Am. Chem. Soc., in 
proof. 
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Fic. 1. Curve showing the observed radial distribution of scattering power in Ni(CO),. The two rows of 
arrows show the distances in the tetrahedral (upper) and square (lower) models based on the observed Ni—C 


and C—O distances. 


molecule. Accordingly without any assumptions 
about the structure of the molecule we know 
that it must contain atomic separations near the 
values 1.18A, 1.82A, 2.97A, and 3.96A. This 
result may be combined with the usual pro- 
cedure’ for interpreting electron diffraction 
photographs in the consideration of square and 
tetrahedral arrangements of the atoms in nickel 
carbonyl. 

The molecular models for which calculations 
of the theoretical intensity of electron scattering 
were made are alike in having four carbonyl 
groups joined to a central nickel atom with the 
line joining the carbon and oxygen atoms in the 
respective carbonyl groups passing through the 
central atom. This restriction is reasonable on 
the basis of a triple bonded carbon monoxide 
joined by a single bond to the nickel atom. 
Models were included in which the oxygen atoms 
are attached to the nickel as well as those 
having carbon atoms adjacent to the nickel atom. 
In the square configuration the carbonyl groups 
are directed toward the corners of a square; in 
the tetrahedral arrangement they point to the 


'S L. Pauling and L. O. Brockway, J. Chem. Phys. 2, 867 
(1934). 


vertices of a regular tetrahedron. Calculations 
were made for each of these arrangements with 
the ratio of carbon-oxygen to nickel-carbon 
distances varying from 0.55 to 0.65. Five of 
these curves are shown in Fig. 2 in which the 
intensity of elastically scattered electrons (calcu- 
lated with atomic scattering powers which are 
constant over the scattering angle) is plotted 
against (47 sin 6/2)/X with an assumed Ni—C 
separation of 2.00A to fix the horizontal scale. 

The general qualitative appearance of the 
photographs is well reproduced by the tetra- 
hedral curves in Fig. 2. The square models with 
a C—O/Ni-—C ratio of 0.625 or greater are 
unsatisfactory because they would require the 
seventh maximum to be appreciably stronger 
than the sixth and they fail to show the subsidi- 
ary fourth maximum following the strong third 
maximum. For values of the ratio below 0.625 
the fourth maximum appears at a position about 
halfway between the third and fifth whereas 
the fourth is observed on the photographs to be 
much closer to the third maximum. 

The quantitative comparison is given in Table 
I. The first column of data is calculated from the 
measured ring diameters and the wave-length 
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Fic. 2. Theoretical intensity of electron scattering by Ni(CO), with constant atomic 
scattering factors. The upper two curves correspond to square models and the lower 


three to tetrahedral models. 


used. The second represents the visually esti- 
mated intensities used for the atomic distribution 
curve. The next five columns give the values of 
the indicated points on the theoretical curves, 
while the last five columns show the correspond- 
ing values for the size of the molecule in terms 
of the assumed distances. The assumed nickel- 
carbon separation of 2.00A is multiplied by the 
a values to give the experimentally determined 
values of the Ni—C distance. The C—O distance 
is then determined by means of the ratios used 
in the respective models. 

The square models in the table are definitely 
inferior to the tetrahedral ones since they show 
an average deviation in thc values for the size 





of the molecule more than three times as great 
as that for the tetrahedral models. An even 
more serious discrepancy is found in the devia- 
tion of the value calculated from the fourth 
maximum. This maximum is subject to the St. 
John effect,'® in which the apparent position of 
an unsymmetric ring is shifted in the direction 
of the greater contrast. In this case the effect 
would lead to too low a value of the molecular 
size. It will be observed in the tetrahedral 
models that the value from the fourth maximum 
is around five percent below the average from 
other rings whereas in the only square model 


19 L. O. Brockway and F. T. Wall, J. Am. Chem. Soc. 56, 
2374 (1934). 
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TABLE I. b 
ae t 
Max. Min 2™8i0. 9/2 INTEN- TETRAHEDRAL SQUARE TETRAHEDRAL SQUARE n 
r SITY X0.650 X0.625 X0.600 X0.625 X6.600 40.650 20.625 40.600 40.625 20.690 t 
1 2.63 5 
2 486 10 430 4.32 4.34 4.05 4.18 0.885 0.889 0.893 0.833 0.860 
35.92 5.35 5.38 5.43 5.50 5.45 004 909.917 929.921 
3 6.92 10 623 630 6.38 6.87 6.80 001 911 — 922 093.983 c 
4 1.74 6.96 7.05 7.15 8.00 000 (9111924 1,034 
4 8.70 6 745 7.55 7.60 8.45 (856) (868) (874) 972 s 
5 9.96 8.97 905 9.20 9.03 9.27 000 = 908-924 906 931 
5 11.09 7 10.12 10.25 10.35 10.10 10.20 913 924 1934 ‘O11 920 
6 12.14 11.20 11.30 11.40 11.35 11.38 923 931 939 .935 938 
6 13.17 3 11.86 12.08 12.25 12.10 12.32 001 918 ~— 1930 019 936 
7 14.08 12.58 12.77 12.92 12.74 13.15 394 907 1918 004.934 
7 15.01 2 13.45 13.58 13.64 13.65 13.83 .896 -904 -908 -909 921 
8 16.25 14.85 15.05 15.20 14.85 15.00 014-926 1936 0141923 
8 17.38 2 15.90 16.12 16.25 16.25 16.22 015 9281935 935.934 
Average -904 -914 924 -917 .939 
Ave. Dev. 007 = 007. S007 023.027 
Final Value Ni —C = 1.82 +0.03A Ni-C 1.808 1.828 1.848 
C—O =1.15+0.02A CcC-O 1.175 1.142 1.109 
Tetrahedral arrangement of atoms 
c 
which shows the fourth ring at all the value is hedral and a square molecule, respectively, of 1 
four percent too high. Since the deviation is in nickel carbonyl which has the Ni—C and C—O ( 
the wrong direction it can be ascribed only to distances given above. The height of the arrows c 
the unsatisfactory character of the square model. is proportional to the products of the scattering t 
Among the tetrahedral models we eliminate powers of the respective pairs of atoms. The 
those having a ratio below 0.625 since they give large, sharp maximum at 2.97A corresponds to € 
values for the C—O distance lower than that the Ni—Oand the C—C distances which coincide \ 
observed in carbon monoxide. Since carbon in the tetrahedral model. The hump at 1.82A is c 
monoxide is observed to have the triple bond the Ni—C separation. The less reliable hump ( 
distance, there is no form of covalent bond which for C—O appears at 1.18A but it should be ‘ 
could exist between the carbon and oxygen observed that the difference between Ni—O and I 
atoms in nickel carbonyl which would have a_ Ni-—C is just 1.15A. t 
smaller interatomic distance. The effect on the The atomic distribution curve also corrobo- 
triple bond itself of forming the extra bond to rates our choice of the tetrahedral arrangement. t 
the nickel atom could scarcely be that of In the experience of Pauling and Brockway"’ the ( 
strengthening the bond and shortening the dis- three distances in the square model at 2.57A, é 
( 


tance. Accordingly, we take the observed dis- 
tance in carbon monoxide, 1.14A, as the lower 
limit for the carbon-oxygen separation in nickel 
carbonyl. 

An upper limit on the C—O/Ni—C ratio may 
be placed at 0.65 since in models above that 
value of the ratio, some of the qualitative 
intensity comparisons occur in the wrong direc- 
tion (such as seventh maximum greater than 
sixth). Accordingly we choose as most probable 
the tetrahedral model with a ratio of 0.63 and 
having Ni—C=1.82+0.03A and C—O=1.15 
+0.02A. 

Inspection of the atomic distribution curve 
now leads to its interpretation in the light of the 
foregoing result. The vertical arrows in Fig. 1 
show all of the interatomic distances in a tetra- 


2.97A, and 3.33A would not combine to give 
the sharp maximum observed at 2.97A. 

Two practically independent methods of treat- 
ing the electron diffraction photographs thus 
lead to the conclusion that nickel carbonyl is a 
tetrahedral molecule with the dimensions given 
above. 


DISCUSSION 


The observed carbon-oxygen distance, 1.15A, 
agrees with the value for carbon monoxide 
molecule and thus supports the previous evidence 
for a similarity of bonding in the carbonyl group 
in this compound and carbon monoxide. Since 
the latter compound is known?’ to resonate 


20. Pauling, J. Am. Chem. Soc. 54, 1000 (1932). 











MOLECULAR STRUCTURE OF NICKEL 


between the double bond and triple bond struc- 
tures the question of finding this same state in 
nickel carbonyl is raised. With a triple bond 
between carbon and oxygen only a single electron 
pair bond could be formed between nickel and 
carbon as indicated in the first of the following 
structures. The formal charges are obtained by 


C 0 
c C 
“4- + *e 
O:::C : MZ: 6:::0: 0:7: Ofc nNa:: ¢ 0 
c c 
o* :0 
I zz 


dividing the electrons in the bond and they 
indicate the directions of the semipolar bonds 
(sometimes represented by arrows). The nickel- 
carbon distance is observed to be 0.18A smaller 
than the sum of the covalent single bond radii. 
It is improbable that the nickel radius is in 
error by this amount or that the single bond 
would be shortened that much by its semipolar 
character (in spite of the rather large negative 
charge on the nickel atom); accordingly, it is 
suggested that this shortening is produced by 
resonance with the structure II in which all of 
the bonds are represented as double electron pair 
bonds. This represents an extreme state since 
the four bonds would perhaps not assume this 
character all at the same time. The double bonds 
are formed by placing two unshared pairs on 
oxygen and bringing one of the unshared pairs 
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of 3d electrons on nickel into the nickel-carbon 
bond. Since both the electrons and eigenfunctions 
are available, nickel does not share with the 
elements of the first row in the periodic table the 
restriction to four electron pair bonds. Structure 
II gives each atom a zero formal charge. Reso- 
nance between I and II explains the observed 
shortening of the nickel-carbon bond and allows 
the carbon-oxygen bond to exhibit a similar 
resonance between the double and triple states 
which it shows in carbon monoxide. For each of 
the bonds in nickel carbonyl the stronger of the 
resonating components determines the distance. 

The directional properties* of structure I are 
determined by the fact that all five of the 3d 
eigenfunctions of neutral nickel are probably 
occupied by unshared pairs and the 4s and three 
4p functions used in forming the four single bonds 
have the usual tetrahedral arrangement. For 
structure II no prediction can be made since the 
problem of four double bonds formed from d‘sp* 
has never been studied. Since structure I prob- 
ably predominates in the resonance it is likely 
that its directional properties determine those 
of the normal state of the molecule; but in fact 
the observed tetrahedral configuration cannot 
be used as a criterion for the proposed resonance. 
It should again be noted that although this is the 
first definitely tetrahedral configuration observed 
in quadricovalent nickel compounds all of the 
other substances studied have been compounds 
of divalent nickel ion. 

We are deeply grateful to Professor Linus 
Pauling for inany suggestions and criticisms. 
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The Conductance of Salt Crystals 


Recently Jost! has derived an expression for the con- 
ductance of salt crystals by postulating in detail the exact 
mechanism involved in the transport of ions through the 
crystal. The expression obtained by Jost does not appear 
adequate to account for the conductance behavior of 
various salts throughout the whole range of temperature, 
and if it can be shown that the same form of expression 
could be obtained by very general considerations, then it 
would appear that the mechanism postulated by Jost is 
not necessarily the correct one. 

We find it possible to obtain the Jost equation by 
assuming that when an ion becomes activated it moves a 
certain average distance x before it loses its mobility. 
The number of ions activated per cubic centimeter per 
second, is given by the expression originally used by one 
of the authors? in the theory of unimolecular reactions 
vnoe~@/RT, Here v is the frequency of the ion vibration in 
the lattice, 2» is the number of ions per cubic centimeter, 
and Q may be identified with E/2+U in Jost’s formula. 
The diffusion coefficient D is then given by the expression 
3x2ve—@2/RT, and when reasonable numerical values are 
introduced the value obtained for D is within the same 
order of magnitude as Jost’s value. 

It is not assumed that the conductance of salt crystals 
can be described completely in terms of expressions as 
simple as the one given above, but we do believe that any 
reasonable mechanism that may be postulated will lead to 
essentially the same type of formula. 

W. H. RopesusH 
T. G. CooKE 
1 Jost, J. Chem. Phys. 1, 466 (1933). 


2 Rodebush, J. Am. Chem. Soc. 45, 606 (1923); J. Chem. Phys. 1, 
440 (1933). 


Vibrations of Benzene and Raman Spectra of Benzene-d 
and Benzene-d, 


In Table I we give the preliminary results of measure- 
ments, which we were performing on benzene-d and p- 
benzene-d2, while Wood! published Raman frequencies of 
benzene-ds. Experimental error cannot account for the 
increase from 2270 (CsH;D) to 2279 (CsH.D2). Several 
lines of p-benzene-d, are more diffuse than the correspond- 
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TABLE I. 








p-CsHsD2 CeHe CsHsD p-C6HsDe2 
3053 cm! | — cm=! 1007 cm™ 1008 cm” 
2279 991 980 979 
1587 849 851 848 
1565 = 638? 
1176 606 602 602 


CsHsD 


3059 cm=! 
2270 
1593 
1571 
1176 


CeHe 


3059 cm=! 
3047 
1606 
1584 
1174 











TABLE II. 








II vj (D) /Tv;(CeHe) 
i t 


exp. exp. (with 
(according v9 = 1178, 
to Wilson) v10 = 


0.723 
1.002 
0.715 
0.541 
0.741 





p-C6eHsDe2 
Ce6De 








ing lines of benzene-d, indicating perhaps an equilibrium 
between the three isomers of benzene-d2. 

Applying a relation between the frequencies of isotope 
molecules, previously derived and proved in the case of 
deuterochloroform,? to Wood’s results and to our own, 
we can test Wilson’s’ correlations of the benzene frequencies 
to the vibrations. The discrepancies between experimental 
and calculated values indicate, that not all correlations can 
be correct. (See Table II.) Indeed, the exchange of the 
correlations of frequencies 9 and 10 yields satisfying 
agreement. The agreement yields a confirmation of Wood’s 
correlations of frequencies 1, 2, 6, 7, 8 and 8. 

We have begun with measurements on other isotopes of 
benzene. Dr. Ph. Gross, Vienna, to whom the authors are 
indebted for directing their attention to the preparation 
of the two benzenes, will shortly report on measurements 
of the vapor pressures. 

O. REDLICH 
W. Srricks 
Technische Hochschule, 
Institut fiir physikalische Chemie, 
Vienna, 
October 25, 1935. 
1R. W. Wood, J. Chem. Phys. 3, 444 (1935). 
20. Redlich, Zeits. f. physik. Chemie B28, 371 (1935). 


3 E. B. Wilson, Jr., Phys. Rev. 46, 146 (1934); cf. also ibid. 45, 706 
(1934). 
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